To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is perf...To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM "footprint" concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.展开更多
The equation of time-domain wave propagation in dispersive media and the explicit beam propagation method are presented in this paper.This method is demonstrated by the short optical pulses in a directional coupler wi...The equation of time-domain wave propagation in dispersive media and the explicit beam propagation method are presented in this paper.This method is demonstrated by the short optical pulses in a directional coupler with second order dispersive effect and shows to be in full agreement with former references.This method is simple,easy and practical.展开更多
The main objective of this paper is to propose a two-phase solution algorithm for solving the Inventory Routing Problem with Time Windows (IRPTW), which has not been excessively researched in the literature. The sol...The main objective of this paper is to propose a two-phase solution algorithm for solving the Inventory Routing Problem with Time Windows (IRPTW), which has not been excessively researched in the literature. The solution approach is based on (a) a simple simulation for the planning phase (Phase I) and (b) the Variable Neighborhood Search Algorithm (VNS) for the routing phase (Phase II). Testing instances are established to investigate algorithmic performance, and the computational results are then reported. The computational study underscores the importance of integrating the inventory and vehicle routing decisions. Graphical presentation formats are provided to convey meaningful insights into the problem.展开更多
In this paper, we present an acceleration strategy for Smoothed Particle Hydrodynamics (SPH) on multi-GPU platform. For single-GPU, we first use a neighborhood search algorithm of compacting cell index combined with...In this paper, we present an acceleration strategy for Smoothed Particle Hydrodynamics (SPH) on multi-GPU platform. For single-GPU, we first use a neighborhood search algorithm of compacting cell index combined with spatial domain characteristics For multi-GPU, we focus on the changing patterns of SPH's computational time. Simple dynamic load balancing algorithm works well because the computational time of each time step changes slowly compared to previous time step. By further optimizing dynamic load balancing algorithm and the communication strategy among GPUs, a nearly linear speedup is achieved in different scenarios with a scale of millions particles. The quality and efficiency of our methods are demonstrated using multiple scenes with different particle numbers.展开更多
An analytical approach for probabilistic evaluation of transient stability of a power system incorporating a wind farm is presented in this study. Based on the fact that the boundary of practical dynamic security regi...An analytical approach for probabilistic evaluation of transient stability of a power system incorporating a wind farm is presented in this study. Based on the fact that the boundary of practical dynamic security region(PDSR) of a power system with double fed induction generators(DFIG) can be approximated by one or few hyper-planes in nodal power injection space, transient stability criterion for given configurations of pre-fault, fault-on and post-fault of a power system is to be expressed by certain expressions of linear combination of nodal injection vector and the transient stability probability(TSP) is further obtained with a much more simplified expression than the complex integral. Furthermore, considering uncertainties of nodal injection power including wind power and load, TSP is calculated analytically by Cornish-Fisher expansion, which can provide reliable evaluation results with high accuracy and much less computing time compared with Monte Carlo simulation. TSP and its visualization can further help operators and planners be aware of the degree of stability or instability and find critical components to monitor and reinforce. Test results on the New England 10-generators and 39-buses power system show the method's effectiveness and significance for probabilistic security assessment.展开更多
Modal identification involves estimating the modal parameters, such as modal frequencies, damping ratios, and mode shapes, of a structural system from measured data. Under the condition that noisy impulse response sig...Modal identification involves estimating the modal parameters, such as modal frequencies, damping ratios, and mode shapes, of a structural system from measured data. Under the condition that noisy impulse response signals associated with multiple input and output locations have been measured, the primary objective of this study is to apply the local or global noise removal technique for improving the modal identification based on the polyreference time domain (PTD) method. While the traditional PTD method improves modal parameter estimation by over-specifying the computational model order to absorb noise, this paper proposes an approach using the actual system order as the computational model order and rejecting much noise prior to performing modal parameter estimation algorithms. Two noise removal approaches are investigated: a "local" approach which removes noise from one signal at a time, and a "global" approach which removes the noise of multiple measured signals simultaneously. The numerical investigation in this article is based on experimental measurements from two test setups: a cantilever beam with 3 inputs and 10 outputs, and a hanged plate with 4 inputs and 32 outputs. This paper demonstrates that the proposed noise-rejection method outperforms the traditional noise-absorption PTD method in several crucial aspects.展开更多
基金supported by the Key Natural Science Foundation(No.41530320)Natural Science Foundation(No.41274121)+1 种基金Natural Science Foundation for young scientist(No.41404093)the Projects on the Development of the Key Equipment of Chinese Academy of Science(No.ZDYZ2012-1-03)
文摘To improve the inversion accuracy of time-domain airborne electromagnetic data, we propose a parallel 3D inversion algorithm for airborne EM data based on the direct Gauss-Newton optimization. Forward modeling is performed in the frequency domain based on the scattered secondary electrical field. Then, the inverse Fourier transform and convolution of the transmitting waveform are used to calculate the EM responses and the sensitivity matrix in the time domain for arbitrary transmitting waves. To optimize the computational time and memory requirements, we use the EM "footprint" concept to reduce the model size and obtain the sparse sensitivity matrix. To improve the 3D inversion, we use the OpenMP library and parallel computing. We test the proposed 3D parallel inversion code using two synthetic datasets and a field dataset. The time-domain airborne EM inversion results suggest that the proposed algorithm is effective, efficient, and practical.
文摘The equation of time-domain wave propagation in dispersive media and the explicit beam propagation method are presented in this paper.This method is demonstrated by the short optical pulses in a directional coupler with second order dispersive effect and shows to be in full agreement with former references.This method is simple,easy and practical.
文摘The main objective of this paper is to propose a two-phase solution algorithm for solving the Inventory Routing Problem with Time Windows (IRPTW), which has not been excessively researched in the literature. The solution approach is based on (a) a simple simulation for the planning phase (Phase I) and (b) the Variable Neighborhood Search Algorithm (VNS) for the routing phase (Phase II). Testing instances are established to investigate algorithmic performance, and the computational results are then reported. The computational study underscores the importance of integrating the inventory and vehicle routing decisions. Graphical presentation formats are provided to convey meaningful insights into the problem.
文摘In this paper, we present an acceleration strategy for Smoothed Particle Hydrodynamics (SPH) on multi-GPU platform. For single-GPU, we first use a neighborhood search algorithm of compacting cell index combined with spatial domain characteristics For multi-GPU, we focus on the changing patterns of SPH's computational time. Simple dynamic load balancing algorithm works well because the computational time of each time step changes slowly compared to previous time step. By further optimizing dynamic load balancing algorithm and the communication strategy among GPUs, a nearly linear speedup is achieved in different scenarios with a scale of millions particles. The quality and efficiency of our methods are demonstrated using multiple scenes with different particle numbers.
基金supported by the National Basic Research Program of China("973"Project)(Grant No.2013CB228204)the National Natural Science Foundation of China(Grant No.51407126)Tianjin Natural Science Foundation(Grant No.15JCQNJC07000)
文摘An analytical approach for probabilistic evaluation of transient stability of a power system incorporating a wind farm is presented in this study. Based on the fact that the boundary of practical dynamic security region(PDSR) of a power system with double fed induction generators(DFIG) can be approximated by one or few hyper-planes in nodal power injection space, transient stability criterion for given configurations of pre-fault, fault-on and post-fault of a power system is to be expressed by certain expressions of linear combination of nodal injection vector and the transient stability probability(TSP) is further obtained with a much more simplified expression than the complex integral. Furthermore, considering uncertainties of nodal injection power including wind power and load, TSP is calculated analytically by Cornish-Fisher expansion, which can provide reliable evaluation results with high accuracy and much less computing time compared with Monte Carlo simulation. TSP and its visualization can further help operators and planners be aware of the degree of stability or instability and find critical components to monitor and reinforce. Test results on the New England 10-generators and 39-buses power system show the method's effectiveness and significance for probabilistic security assessment.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51079134 and 51009124)the NSFC Major International Joint Research Project (Grant No. 51010009)+2 种基金the Program for Changjiang Scholars and Innovative Research Team in University (Grant No. PCSIRT 1086)the Natural Science Foundation of Shandong Province(Grant Nos. ZR2011EEQ022 and 2009ZRA05100)the Fundamental Research Funds for the Central Universities (Grant Nos. 27R1202008A and27R1002076A)
文摘Modal identification involves estimating the modal parameters, such as modal frequencies, damping ratios, and mode shapes, of a structural system from measured data. Under the condition that noisy impulse response signals associated with multiple input and output locations have been measured, the primary objective of this study is to apply the local or global noise removal technique for improving the modal identification based on the polyreference time domain (PTD) method. While the traditional PTD method improves modal parameter estimation by over-specifying the computational model order to absorb noise, this paper proposes an approach using the actual system order as the computational model order and rejecting much noise prior to performing modal parameter estimation algorithms. Two noise removal approaches are investigated: a "local" approach which removes noise from one signal at a time, and a "global" approach which removes the noise of multiple measured signals simultaneously. The numerical investigation in this article is based on experimental measurements from two test setups: a cantilever beam with 3 inputs and 10 outputs, and a hanged plate with 4 inputs and 32 outputs. This paper demonstrates that the proposed noise-rejection method outperforms the traditional noise-absorption PTD method in several crucial aspects.