Numerical storm surge model with higher order finite difference method of lines for the coast of Bangladesh

In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.

RBFs Meshless Method of Lines for the Numerical Solution of Time-Dependent Nonlinear Coupled Partial Differential Equations

In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.

A New Approach Based on Iterative Method for the Characterization of a Micro-Strip Line with Thick Copper Conductor

In this work, we applied two electromagnetic models for the characterization of a planar structure including a flat, thick copper conductor. Indeed the first model is consisted by modeling two metal ribbons without bulkiness, placed one above the other at a distance of h2 equal to the thickness of the thick conductor. This approach has been implemented and tested by the iterative method. The results of simulations have been compared with those calculated by the Ansoft HFSS software, and they are in good concordance, validating the method of analysis used. The second model is based on the calculation of the effective permittivity of the medium containing the thick conductor. This medium consists of a metallic region of complex relative permittivity , the rest of this medium is filled with air er2 = 1. The effective permittivity eeff calculated from these two relative permittivity er2 and . Comparing the simulation results of this new formulation of the iterative method with those calculated by the software Ansoft HFSS shows that they are in good matching which validates the second model.

A CALCULATING METHOD OF THE KILL PROBABILITY ATTACK AREA FOR AAM

This paper provides a calculating method which can be used in calculation of the kill probability attack area for every AAM. At first, attack area of AAM and kill probability of every characteristic point are obtained by combining trajectory calculation with kill probability calculation. Then, coordinates of a fire point relative to standard kill probability value in terms of standardization method are found. At last, equivalent kill probability curve equations are formulated by means of curve fitting method.

Simulation on Calculation Accuracy of Three Methods for Live Line Measuring the Parameters of Transmission Lines with Mutual Inductance

Live line measurement methods can reduce the loss of power outages and eliminate interference. There are three live line measurement methods including integral method, differential method and algebraic method. A simulation model of?two coupled parallel transmission lines spanning on the same towers is built in PSCAD and the calculation errors of these three methods are compared with different sampling frequencies by using of Matlab. The effect of harmonic on calculation is also involved. The simulation results indicate that harmonic has the least effect on the algebraic method which provides stable result and small error.

A New Modification of the Method of Lines for First Order Hyperbolic PDEs

A new modification of the Method of Lines is proposed for the solution of first order partial differential equations. The accuracy of the method is shown with the matrix analysis. The method is applied to a number of test problems, on uniform grids, to compare the accuracy and computational efficiency with the standard method.

Two Methods for Calculating the Size Distribution of Ferritin’s Outer Diameter

Ferritin is a kind of iron-storage protein widely found in animals and plants. The dynamic light scattering (Dynamic Light Scattering) method is used in the laboratory to determine the ferritin size. This paper presents two methods for calculating the outer diameter size distribution of ferritin, both of which assume that ferritin is approximately spherical. The ferritin data file was obtained from the PDB website and was calculated using the coordinate data of the amino acids to which the amino acids belong. The first method is based on the calculation of the sphere center;the second method is based on the method of the farthest distance atom pair. The outer diameter size distribution curves obtained by the two methods are basically consistent with the experimental methods. The paper also compares the calculation results and performance of the two methods. Both methods are versatile and can be used to calculate the size distribution of the globular proteins.

Method of Lines for the Chiral Nonlinear Schrödinger Equation

In this paper, we solve chiral nonlinear Schrodinger equation (CNSE) numerically. Two numerical methods are derived using the explicit Runge-Kutta method of order four and the linear multistep method (Predictor-Corrector method of fourth order). The resulting schemes of fourth order accuracy in spatial and temporal directions. The CNSE is non-integrable and has two kinds of soliton solutions: bright and dark soliton. The exact solutions and the conserved quantities of CNSE are used to display the efficiency and robustness of the numerical methods we derived. Interaction of two bright solitons for different parameters is also displayed.

A Slip-Line Method for Calculating Extrusion Force of Steel Helmet with Cold Extrusion Moulding

This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution of slip-line is established based on the plastic streamline. Theextrusion force of plastic moulding of the steel helmet is calculated by using of slip-line method.Furthermore, an applied example is given.

THEORY AND METHOD FOR WETLAND BOUNDARY DELINEATION

Based on the analysis of the subjectivity of wetland boundary criteria and their causes at present, this paper suggested that, under the condition that the mechanism of wetland formation process has not been understood, "black box" method of System Theory can be used to delineate wetland boundaries scientifically. After analyzing the difference of system construction among aquatic habitats, wetlands and uplands, the lower limit of rooted plants was chosen as the lower boundary criterion of wetlands. Because soil diagnostic horizon is the result of the long-term interaction among all environments, and it is less responsive than vegetation to short-term change, soil diagnostic horizon was chosen as the indicator to delineate wetland upper boundary, which lies at the thinning-out point of soil diagnostic horizon. Case study indicated that it was feasible using the lower limit of rooted plants and the thinning-out point of soil diagnostic horizon as criteria to delineate the lower and upper boundaries of wetland. In the study area, the thinning-out line of albic horizon was coincident with the 55.74m contour line, the maximum horizon error was less than 1m, and the maximum vertical error less than 0.04m. The problem on wetland definition always arises on the boundaries. Having delineated wetland boundaries, wetlands can be defined as follows: wetlands are the transitional zones between uplands and deepwater habitats, they are a kind of azonal complex that are inundated or saturated by surface or ground water, with the lower boundary lying at the lower limit of rooted plants, and the upper boundary at the thinning-out line of upland soil diagnostic horizon.

Erratum to:“An Approximate Method for Calculating the Fluid Force and Response of A Circular Cylinder at Lock-in” [China Ocean Engineering,22(3),2008,371-384]

Because of my carelessness,Eq.（1）in the paper "An approximate method for calculating the fluid force and response of a circular cylinder at lock-in"(China Ocean Engineering,22（3）,2008,pp.373)should be f’-1.0/U’-5.0=f’;-1.0/5.75f’;-5.0,not f’=U’/5.75. My apology is hereby given.

ON THE LIMIT CASE OF THE STEP-REDUCTION METHOD FOR CALCULATING NON-UNIFORM BEAM WITH VARIOUS SECTIONS

In this paper, the step reduction method is discussed, which was advanced by Prof. Yeh Kai-yuan for calculating a non-uniform beam with various sections. The following result is proved. The approximate solution by this method approaches the true solution if the number of the steps approaches the infinity. However, the measure of the error between the limit solution and the ture solution is not in the pure mathematics sense but in the mechanics sense.

Analysis and Application of the Synthetic Relative Measuring Method in On-Line Monitoring for Capacitive Equipment in Power Systems

On-line measurement for dielectric loss angle can effectively monitor the insulation condition of capacitive equipment in power systems. Synthetic relative measuring methods not only markedly overcome the shortcomings of traditional absolute measuring methods but also greatly improve the accuracy of dielectric loss angle measurement. However, synthetic relative measuring methods based on two or three pieces of capacitive equipment do not have the characteristic of generality. In this paper, a principle of synthetic relative measuring method is presented. The example of application for synthetic relative methods based on three and four pieces of capacitive equipment running in the same phase is taken to present the failure judgment matrices for N pieces of equipment. According to these matrices, the fault condition of N pieces of capacitive equipment can be watched, which is more general. Then some problems needing to be concerned along with two diagnostic methods used in diagnostic system are introduced. Finally, two programmable flow charts for the two methods are given and corresponding examples demonstrate their feasibility in practice.

Type Synthesis of 4-DOF Parallel Kinematic Mechanisms Based on Grassmann Line Geometry and Atlas Method

Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms（PKMs）, but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom（DOFs） and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R（T： translational DOF; R： rotational DOF）, are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs＇ type synthesis. Moreover, this method can provide theoretical guidance for other PKMs＇ type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods＇ problems in terms of complicated, not intuitive and unsuitable for practical application.

Calculating Method for Influence of Material Flow on Energy Consumption in Steel Manufacturing Process

From the viewpoint of systems energy conservation, the influences of material flow on its energy consumption in a steel manufacturing process is an important subject. The quantitative analysis of the relationship between material flow and the energy intensity is useful to save energy in steel industry. Based on the concept of standard material flow diagram, all possible situations of ferric material flow in steel manufacturing process are analyzed. The expressions of the influence of material flow deviated from standard material flow diagram on energy consumption are put forward.

SMOOTHING NEWTON ALGORITHM FOR THE CIRCULAR CONE PROGRAMMING WITH A NONMONOTONE LINE SEARCH

In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space with the circular cone. Based on the relationship between the circular cone and the second-order cone(SOC), we reformulate the CCP problem as the second-order cone problem(SOCP). By extending the nonmonotone line search for unconstrained optimization to the CCP, a nonmonotone smoothing Newton method is proposed for solving the CCP. Under suitable assumptions, the proposed algorithm is shown to be globally and locally quadratically convergent. Some preliminary numerical results indicate the effectiveness of the proposed algorithm for solving the CCP.

GLOBAL COVERGENCE OF THE NON-QUASI-NEWTON METHOD FOR UNCONSTRAINED OPTIMIZATION PROBLEMS

In this paper, the non-quasi-Newton＇s family with inexact line search applied to unconstrained optimization problems is studied. A new update formula for non-quasi-Newton＇s family is proposed. It is proved that the constituted algorithm with either Wolfe-type or Armijotype line search converges globally and Q-superlinearly if the function to be minimized has Lipschitz continuous gradient.

Analysis and management strategies for the low voltage problems of rural power grid based on the distribution network flow calculation method

For a long time, because of the lack of investment capital and enough attentions, the overall constructions of rural power grid were far behind than the urban power grid in Chongqing Jiangbei Power Company. The low voltage problems were highlighted in the rural power grid due to the characteristics of rural power grid. Using the distribution network flow calculation method, we evaluated the low voltage problems of the rural power grid which belongs to Chongqing Jiangbei Power Company. In addition, we collected the data of distribution transformers in electricity consumption peak period. Some practical management strategies were proposed by the analysis and evaluation of potential and appeared low voltage problems.

Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Study on the automatic process of line heating for pillow shape plate

This paper focuses on the process for pillow shape plate by line heating technique, which is widely applied in the production of ship hull. Based on the analysis of primary parameters and experimental data in line heating process, the amount of local contraction generated by line heating has been illustrated. Then, combining with the computational result of local deformation determined by shell plate development, an optimization method for line heating parameters has been studied. This prediction system may provide rational arrangements of heating lines and technical parameters of process. By integrating the prediction system into the line heating robot for pillow shape plate, the automatic process of line heating for pillow shape plate can be achieved.