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.

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.

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.

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.

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.

RESEARCH ON THE EVALUATION METHOD OF FRICTION AND LUBRICATION IN DEEP DRAWING

The deformation characteristic of bland in deep drawing is discussed. It is pointed out that the friction and lubrication conditions in for drawing are different from that in mechanical motion or machine work or other plastic process. The common test methods in laboratories are analyzed. It shows that though all those test methods can test the friction coefficient, the probe test method is most suitable for the research of friction and lubrication and the process in deep drawing, for this method is identical with the actual work condition either from the test principle or deformation status of the blank. Last the successful application in the deep drawing simulator newly developed the the probe method are intro- duced in detail.

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.

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.

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.

Method of Lines for Third Order Partial Differential Equations

The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.

Recent improvements of actuator line-large-eddy simulation method for wind turbine wakes

In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.

An Improved Method for Doppler Wind and Thermo dynamic Retrievals

A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: the beginning and ending times of the data assimilation period, simultaneously from two successive volume scans by using the weak form constraints provided by the mass continuity and vorticity equations. As the retrieved wind fields are expressed by Legendre polynomial expansions at the beginning and ending times, the time tendency term in the vorticity equation can be conveniently formulated, and the retrieved winds can be compared with the radar observed radial winds in the cost function at the precise time and position of each radar beam. In the second step, the perturbation pressure and temperature fields at the middle time are then derived from the retrieved wind fields and the velocity time tendency by using the weak form constraints provided by the three momentum equations. The merits of the new method are demonstrated by numerical experiments with simulated radar observations and compared with the traditional least squares methods which consider neither the precise observation times and positions nor the velocity time tendency. The new method is also applied to real radar data for a heavy rainfall event during the 2001 Meiyu season in China.

An Efficient Approach for Simulation of Water Levels due to the Nonlinear Interaction of Tide and Surge Along the Coast of Bangladesh

The ultimate goal and highlight of this paper are to explore water levels along the coast of Bangladesh efficiently due to the nonlinear interaction of tide and surge by employing the method of lines(MOLs)with the aid of newly proposed RKAHeM(4,4)technique.In this regard,the spatial derivatives of shallow water equations(SWEs)were discretized by means of a finite difference method to obtain a system of ordinary differential equations(ODEs)of initial valued with time as an independent variable.The obtained system of ODEs was solved by the RKAHeM(4,4)technique.One-way nested grid technique was exercised to incorporate coastal complexities closely with minimum computational cost.A stable tidal oscillation was produced over the region of interest by applying the most influential tidal constituent M2 along the southern open boundary of the outer scheme.The newly developed model was applied to estimate water levels due to the non-linear interaction of tide and surge associated with the catastrophic cyclone April 1991 along the coast of Bangladesh.The approach employed in the study was found to perform well and ensure conformity with real-time observations.

STUDY ON THREE-DIMENSIONAL FINITE BODIES CONTAINING CRACKS USING THE FINITE ELEMENT METHOD OF LINES

The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.

Investigation on the Dynamic Responses of a Truss Spar Platform for Different Mooring Line Groups

The dynamic responses of any floating platform arc dependent on the mass, stiffness and damping characteristics of the body as well as mooring system. Therefore, it is very essential to study the effect of individual contributions to the system that can finally help to economise their cost. This paper focuses on the effect of mooring stiffness on the responses of a truss spar platform, obtained by different grouping of lines. The study is part of our present researches on mooring systems which include the effect of line pretension, diameter and azimuth angles. The platform is modelled as a rigid body with three degrees-of-freedom and its motions are analyzed in time-domain using the implicit Newmark Beta technique. The mooring lines restoring force-excursion relationship is evaluated using a quasi-static approach. It is observed that the mooring system with lines arranged in less number of groups exhibits better performance in terms of the restoring forces as well as mean position of platform. However, the dynamic motions of platform remain unaffected for different line groups.

Multi-Branch Fault Line Location Method Based on Time Difference Matrix Fitting

The distribution network exhibits complex structural characteristics,which makes fault localization a challenging task.Especially when a branch of the multi-branch distribution network fails,the traditional multi-branch fault location algorithm makes it difficult to meet the demands of high-precision fault localization in the multi-branch distribution network system.In this paper,the multi-branch mainline is decomposed into single branch lines,transforming the complex multi-branch fault location problem into a double-ended fault location problem.Based on the different transmission characteristics of the fault-traveling wave in fault lines and non-fault lines,the endpoint reference time difference matrix S and the fault time difference matrix G were established.The time variation rule of the fault-traveling wave arriving at each endpoint before and after a fault was comprehensively utilized.To realize the fault segment location,the least square method was introduced.It was used to find the first-order fitting relation that satisfies the matching relationship between the corresponding row vector and the first-order function in the two matrices,to realize the fault segment location.Then,the time difference matrix is used to determine the traveling wave velocity,which,combined with the double-ended traveling wave location,enables accurate fault location.