Numerical Modeling the Flood and Pollutant Transport Processes in Residential Areas with Different Land Use Types

Large-scale flooding causes widespread disaster, and harmful pollutant concentration in water following flood affects public safety and the environment. In this study, a numerical model for solving the 2D shallow water equations and the solute transport equation is proposed to simulate overland flood and pollutant transport caused by floods. The present model is verified by comparing the predictions with the analytical solutions and simulation results; sufficiently high computational accuracy is achieved. The model is also used to simulate flood inundation and pollution spread in the area of Hun and Taizi Lane (HTL) in China due to river dike breaches; the results show that the coupling model has excellent performance for simulating the flooding process and the temporal and spatial distribution of pollutants in urban or rural areas. We use remote sensing techniques to acquire the land coverage in the area of HTL based on Landsat TM satellites. The impacts of changed land use on mitigation of flooding waves and pollutant spread are investigated; the results indicate that the land cover changes have an obvious influence on the evolution process of flood waves and pollutant transport in the study areas, where the transport of pollutants is very dynamic during flood inundation in HTL area. Furthermore, the motion of pollutants considering anisotropic diffusion is more reasonable than that due to isotropic dispersion in simulating pollutant transport associated with the flood in urban or farmland environments.

Investigation of Finite-Difference Schemes for the Numerical Solution of a Fractional Nonlinear Equation

The article discusses different schemes for the numerical solution of the fractional Riccati equation with variable coefficients and variable memory, where the fractional derivative is understood in the sense of Gerasimov-Caputo. For a nonlinear fractional equation, in the general case, theorems of approximation, stability, and convergence of a nonlocal implicit finite difference scheme (IFDS) are proved. For IFDS, it is shown that the scheme converges with the order corresponding to the estimate for approximating the Gerasimov-Caputo fractional operator. The IFDS scheme is solved by the modified Newton’s method (MNM), for which it is shown that the method is locally stable and converges with the first order of accuracy. In the case of the fractional Riccati equation, approximation, stability, and convergence theorems are proved for a nonlocal explicit finite difference scheme (EFDS). It is shown that EFDS conditionally converges with the first order of accuracy. On specific test examples, the computational accuracy of numerical methods was estimated according to Runge’s rule and compared with the exact solution. It is shown that the order of computational accuracy of numerical methods tends to the theoretical order of accuracy with increasing nodes of the computational grid.

Parametric and V&V study in a fundamental CFD process: revisiting the lid-driven cavity flow

Purpose The openings on aircraft structures can be modeled from an aerodynamical point of view as lid-driven cavities (LDC). This paper aims to show the primary verification and validation (V&V) process in computational fluid dynamics (CFD, and to investigate the influences of numerical settings on the efficiency and accuracy for solving the LDC problem. Design/methodology/approach To dig into the details of CFD approaches, this paper outlines the design, implementation, V&V and results of an efficient explicit algorithm. The parametric study is performed thoroughly focusing on various iteration methods, grid density discretization terms and Reynolds number effects. Findings This study parameterized the numerical implementation which provides empirical insights into how computational accuracy and efficiency are affected by changing numerical settings. At a low Reynolds number (not over 1,000), the time-derivative preconditioning is necessary, and k = 0.1 can be the optimal value to guarantee the efficiency, as well as the stability. A larger artificial viscosity (c = 1/16) would relieve the calculating oscillation issue but proportionally increase the discretization error. Furthermore, the iteration method and the mesh quality are two key factors that affect the convergence efficiency, thus need to be selected “wisely”. Practical implications The study shows how numerical implementation can enhance an accurate and efficient solution. This workflow can be used to determine the best parameter settings whenever CFD researchers applying this LDC problem as a complementary design tool for testing newly developed solvers. Originality/value The studied LDC problem is representative of CFD analysis in real aircraft structures. These numerical simulations provide a cost-effective and convenient tool to understand the parameter sensitivity, solution receptivity and physics of the CFD process.

The Improved Element-Free Galerkin Method for 3D Helmholtz Equations

The improved element-free Galerkin (IEFG) method is proposed in this paper for solving 3D Helmholtz equations. The improved moving least-squares (IMLS) approximation is used to establish the trial function, and the penalty technique is used to enforce the essential boundary conditions. Thus, the final discretized equations of the IEFG method for 3D Helmholtz equations can be derived by using the corresponding Galerkin weak form. The influences of the node distribution, the weight functions, the scale parameters of the influence domain, and the penalty factors on the computational accuracy of the solutions are analyzed, and the numerical results of three examples show that the proposed method in this paper can not only enhance the computational speed of the element-free Galerkin (EFG) method but also eliminate the phenomenon of the singular matrix.

The Improved Interpolating Dimension Splitting Element-Free Galerkin Method for 3D Potential Problems

This paper proposes the improved interpolating dimension splitting element-free Galerkin (IIDSEFG) method based on the nonsingular weight function for three-dimensional (3D) potential problems. The core of the IIDSEFG method is to transform the 3D problem domain into a series of two-dimensional (2D) problem subdomains along the splitting direction. For the 2D problems on these 2D subdomains, the shape function is constructed by the improved interpolating moving least-squares (IIMLS) method based on the nonsingular weight function, and the finite difference method (FDM) is used to couple the discretized equations in the direction of splitting. Finally, the calculation formula of the IIDSEFG method for a 3D potential problem is derived. Compared with the improved element-free Galerkin (IEFG) method, the advantages of the IIDSEFG method are that the shape function has few undetermined coefficients and the essential boundary conditions can be executed directly. The results of the selected numerical examples are compared by the IIDSEFG method, IEFG method and analytical solution. These numerical examples illustrate that the IIDSEFG method is effective to solve 3D potential problems. The computational accuracy and efficiency of the IIDSEFG method are better than the IEFG method.

Probability Density Evolution Algorithm for Stochastic Dynamical Systems Based on Fractional Calculus

With the increasing load and speed of trains, the problems caused by various random excitations (such as safety and passenger comfort) have become more prominent and thus arises the necessity to analyze stochastic dynamical systems, which is important in both academic and engineering circles. The existing analysis methods are inadequate in terms of computational accuracy, computational efficiency, and applicability in solving complex problems. For that, a new efficient and accurate method is used in this paper, suitable for linear and nonlinear random vibration analysis of large structures as well as static and dynamic reliability assessment. It is the direct probability integration method, which is extended and applied to the random vibration reliability analysis of dynamical systems. Dynamical models of the dynamic system and coupled system “three-car vehicle-rail-bridge” are established, the time-varying differential equations of motion are derived in detail, and the dynamic response of the system is calculated using the explicit Newmark algorithm. The simulation results show the influence of the number of representative points on the smoothness of the image of the probability density function and the accuracy of the calculation results.

Numerical Modeling on Crack Propagation Based on a Multi-Grid Bond-Based Dual-Horizon Peridynamics

Peridynamics (PD) is a novel nonlocal theory of continuum mechanics capable of describing crack formation and propagation without defining any fracture rules in advance. In this study, a multi-grid bond-based dual-horizon peridynamics (DH-PD) model is presented, which includes varying horizon sizes and can avoid spurious wave reflections. This model incorporates the volume correction, surface correction, and a technique of nonuniformity discretization to improve calculation accuracy and efficiency. Two benchmark problems are simulated to verify the reliability of the proposed model with the effect of the volume correction and surface correction on the computational accuracy confirmed. Two numerical examples, the fracture of an L-shaped concrete specimen and the mixed damage of a double-edged notched specimen, are simulated and analyzed. The simulation results are compared against experimental data, the numerical solution of a traditional PD model, and the output from a finite element model. The comparisons verify the calculation accuracy of the corrected DH-PD model and its advantages over some other models like the traditional PD model.

Determination of Local Stresses and Strains within the Notch Strain Approach: Efficient Implementation of Notch Root Approximations

An estimation of the elastic-plastic stress state using elasticity-theoretical input data is an essential part of the service life estimation with the local strain approach in general and a German guideline based on it, in particular. This guideline uses two different notch root approximations (an extended version of Neuber’s rule and an approach according to Seeger and Beste) for this estimation. Both require the implementation of Newton’s method to be iteratively solved. However, many options are left open to the user concerning implementation in program code. This paper discusses ways in which notch root approximation methods can be implemented efficiently for use in software systems and elaborates an application recommendation. The following aspects and their influence on the computational accuracy and performance of Newton’s method are considered in detail: influence of the formulation of the root finding problem, determination of the derivative required for Newton’s method and influence of the termination criterion. The investigation shows that the advice given in the abovementioned guideline indeed leads to a conservative implementation. By carefully considering the investigated aspects, however, the computational performance can be increased by approximately a factor of 2–3 without influencing the accuracy of the service life estimation.

Toward Computational Accuracy in Realistic Systems to Aid Understanding of Field-Level Water Quality Issues

Contemplating what will unfold in this new decade and those after, it is not difficult to imagine the increasing importance of conservation and protection of clean water supplies. A worrying but predictable offshoot of humanity’s technological advances is the seemingly ever-increasing chemical load burdening our waterways. In this perspective are presented a few modest areas where computational chemistry modelling could provide benefit to these efforts by harnessing the continually improving computational power available to the field. In the acute event of a chemical spill incident, true quantum-chemistry-based predictions of physicochemical properties and surface-binding behaviors can be used to help decision making in remediating the spill threat. The chronic burdens of microplastics and perfluorinated “forever chemicals” can also be addressed with computational modelling to fill the gap between feasible laboratory experiment timescales and the much-longer material lifetimes. For all of these systems, field-level accuracy models will avail themselves as the model computational systems are able to incorporate more realistic features that are relevant to water quality issues.

Properties Analysis of Disk Spring with Effects of Asymmetric Variable Friction

The primary objective of this fundamental research is to investigate the mechanical properties of the disk spring when the friction at the contact edges is asymmetric and varies with the load. The contact mechanics study shows that the static friction and static friction coefficient on fractal surfaces change depending on the normal load. In this paper, a fractal contact model based on the W-M function is used to explore the connection between the static friction and the normal load. Subsequently, taking into account the asymmetry of the contact surface at the edge, the variable static friction coefficient is brought into the existing model to obtain an improved static model of the disk spring. Different fractal dimensions, frictional states and free heights are considered under quasi-static loading condition, the relative errors between this paper and the method using Coulomb friction are also calculated, and experimental validation was performed. The static stiffness and force hysteresis of the disk spring for different forms of asymmetric variable friction are discussed. It is shown that using the variable friction model can improve the computational accuracy of the disk spring model under small loads and help to improve the design and control accuracy of preload and vibration isolation equipment using the disk spring as a component.