CN113191105A - Electrical simulation method based on distributed parallel operation method - Google Patents

Electrical simulation method based on distributed parallel operation method Download PDF

Info

Publication number
CN113191105A
CN113191105A CN202110302150.8A CN202110302150A CN113191105A CN 113191105 A CN113191105 A CN 113191105A CN 202110302150 A CN202110302150 A CN 202110302150A CN 113191105 A CN113191105 A CN 113191105A
Authority
CN
China
Prior art keywords
matrix
iteration
solution
convergence
calculation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110302150.8A
Other languages
Chinese (zh)
Inventor
梁文毅
高秋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN202110302150.8A priority Critical patent/CN113191105A/en
Publication of CN113191105A publication Critical patent/CN113191105A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/32Circuit design at the digital level
    • G06F30/33Design verification, e.g. functional simulation or model checking
    • G06F30/3308Design verification, e.g. functional simulation or model checking using simulation
    • G06F30/331Design verification, e.g. functional simulation or model checking using simulation with hardware acceleration, e.g. by using field programmable gate array [FPGA] or emulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46Multiprogramming arrangements
    • G06F9/50Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5005Allocation of resources, e.g. of the central processing unit [CPU] to service a request
    • G06F9/5027Allocation of resources, e.g. of the central processing unit [CPU] to service a request the resource being a machine, e.g. CPUs, Servers, Terminals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46Multiprogramming arrangements
    • G06F9/50Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5061Partitioning or combining of resources
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/02CAD in a network environment, e.g. collaborative CAD or distributed simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses an electrical simulation method based on a distributed parallel operation method, which specifically comprises the following steps: s1, updating Jacobian matrix elements and right-end items; s2, preprocessing the Jacobian matrix to reduce the condition number of the coefficient matrix; s3, solving a matrix equation by adopting a Krylov subspace method; s4, in each process, according to the task distributed by the main process, corresponding calculation is carried out; s5, collecting the data after the calculation of each process into a main process, and processing to obtain the current iteration solution vector; s6, convergence judgment is carried out, when the matrix is solved for convergence, the step S7 is carried out, otherwise, the step S3 is carried out, and solution vectors are obtained according to the current calculation; and S7, judging whether the current nonlinear iterative solution convergence condition is satisfied, and when the convergence condition is not satisfied, re-solving until the nonlinear solution convergence. The invention optimizes the design, and can obviously improve the solving efficiency of the simulation of a complex large system by fully utilizing hardware computing resources in the local area network.

Description

Electrical simulation method based on distributed parallel operation method
Technical Field
The invention relates to the technical field of computer virtual simulation, in particular to an electrical simulation method based on a distributed parallel operation method.
Background
The electrical simulation technology is widely applied to the design of complex electrical systems such as aviation, aerospace, ships, weaponry, rail transit and the like, and currently, a general electrical system simulator generally adopts a direct method based on LU triangular decomposition to perform system analog simulation. For a large sparse matrix, there is a great technical obstacle to realizing parallel operation of LU triangular decomposition, so that a current general solver generally adopts a single-node technology to perform electrical system simulation. With the development of multi-power and full-power technologies, independent power supply systems are more and more complex, and meanwhile, with the continuous progress of modeling technologies, the complexity of simulation models of electrical systems is higher and higher, and due to the fact that single-node hardware configuration is adopted, the requirements of higher and higher complex system simulation are difficult to meet due to low simulation efficiency.
Disclosure of Invention
The invention aims to provide an electrical simulation method based on a distributed parallel operation method, which aims to solve the problems in the background technology and improve the convergence performance and the solving efficiency of the simulation of a complex electrical system.
In order to achieve the purpose, the invention provides the following technical scheme:
an electrical simulation method based on a distributed parallel operation method specifically comprises the following steps:
s1, reading a netlist file by a solver during initialization, wherein the netlist file comprises detailed information of element types, element pin node connection information, element model parameters, simulation step lengths, convergence parameters and an integration method, and constructing a Jacobian initial matrix based on a node voltage method; when Newton-Raphson nonlinear iteration is carried out, updating the elements of the Jacobian matrix and the right-end term according to the current solution vector state at any iteration moment;
s2, preprocessing the Jacobian matrix to reduce the condition number of the coefficient matrix;
s3, solving a matrix equation by adopting a Krylov subspace method;
s4, in each process, according to the task distributed by the main process, corresponding calculation is carried out;
s5, collecting the data after the calculation of each process into a main process, and carrying out comprehensive processing on the data to obtain the current iteration solution vector;
s6, carrying out convergence judgment according to iteration results, entering the step S7 when the matrix solution is converged, otherwise returning to the step S3, updating the Jacobian matrix and the right-end item again according to the solution vector obtained by current calculation, and continuing to iteratively solve the current matrix until iteration is converged; when the iteration times exceed the maximum iteration times and still do not converge, warning information is given;
s7, comparing the state of the solution vector and the right-end term with a convergence criterion, judging whether the current nonlinear iteration solution convergence condition is satisfied, substituting the solution vector obtained by iteration into a Jacobian matrix and the right-end term when the convergence condition is not satisfied, updating equation coefficients and excitation elements again, and solving equations until the nonlinear solution convergence; and repeating the nonlinear iteration process until all the solving tasks are completed.
Preferably, the initial state of the solution vector in step S1 is obtained from historical solution data before entering the iteration, and the solution vector obtained for each iteration is obtained after entering the non-linear iteration.
Preferably, in step S1, in the newton-raphson nonlinear iterative calculation process, a step-changing algorithm is used.
Preferably, after the preprocessing algorithm is adopted in step S2, the right-end item is processed synchronously according to the selected preprocessor.
Preferably, in step S3, obtaining an expression of an equation iterative solution according to a Krylov subspace method, and performing solution operation by using a sparse matrix; when the solved object matrix is large, the calculation load is decomposed according to the hardware resources in the local area network and is distributed to each process of the local area network.
Compared with the prior art, the invention has the beneficial effects that: the invention provides an electrical simulation technology based on a distributed parallel operation method in order to improve the solving efficiency of the simulation of a complex electrical system. The method is based on a Krylov subspace iteration method to realize distributed parallel operation. Because the main computational load of the Krylov subspace iteration method is matrix vector multiplication, the sparse matrix can be fully utilized to further improve the computational efficiency. The distributed parallel computing realizes multi-process parallel computing based on MPI parallel environment. The invention optimizes the design, and can obviously improve the solving efficiency of the simulation of a complex large system by fully utilizing hardware computing resources in the local area network.
Drawings
Fig. 1 is a schematic flow chart of an electrical simulation method based on a distributed parallel operation method.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, the present invention provides a technical solution: an electrical simulation method based on a distributed parallel operation method specifically comprises the following steps:
an electrical simulation method based on a distributed parallel operation method specifically comprises the following steps:
s1, reading a netlist file by a solver during initialization, wherein the netlist file comprises detailed information of element types, element pin node connection information, element model parameters, simulation step lengths, convergence parameters and an integration method, and constructing a Jacobian initial matrix based on a node voltage method; when Newton-Raphson nonlinear iteration is carried out, updating the elements of the Jacobian matrix and the right-end term according to the current solution vector state at any iteration moment; in step S1, the initial state of the solution vector is obtained from the historical solution data before entering the iteration, and the solution vector is obtained for each iteration after entering the nonlinear iteration; in the Newton-Raphson nonlinear iterative computation process, a variable step length algorithm is adopted to improve the computation efficiency;
s2, preprocessing the Jacobian matrix to reduce the condition number of the coefficient matrix; after a preprocessing algorithm is adopted, the right-end item is synchronously processed according to the selected preprocessing son;
s3, solving a matrix equation by adopting a Krylov subspace method;
obviously, the numerical solution problem here can be converted to solve a general system of linear equations:
Ax=b,
wherein the matrix A is formed by Rn×nAnd vector b ∈ RnGiven, x ∈ RnIs an unknown vector; setting a coefficient matrix A as a nonsingular large sparse matrix;
taking an initial vector x0∈RnSolving the system of equations is equivalent to solving
Ay=r0,x=x0+y,r0=b-Ax0
Constructing a Krylov subspace
Figure RE-RE-GDA0003125766630000051
Then one can be obtained
Figure RE-RE-GDA0003125766630000052
In a sense that it is the system of equations Ay r0Is best approximated by
y=A-1r0≈yk=h0r0+h1Ar0+hkAk-1r0
The approximate solution of the corresponding original equation set Ax ═ b is xk=x0+yk
In the above calculation process, the matrix inversion operation is replaced by a matrix polynomial, and in the calculation process, the matrix A and the vector A are used as the basisk-1r0Multiplication recursion to obtain Akr0(ii) a Because the matrix A has very obvious sparse characteristics and the parallel calculation of the sparse matrix and the vector multiplication is easy to realize, the sparse matrix operation is adopted; by utilizing the characteristic of the sparse matrix, the storage space is reduced, and the calculation requirement is reduced;
the Sparse matrix Storage format adopted here is Row Compressed Storage (CSR); here, the CSR format is taken as an example to describe a storage method of the sparse matrix; setting a sparse matrix A, creating three arrays when adopting a CSR format, wherein one array is a floating-point array V, the other two arrays are integer arrays C and R, and the three arrays are respectively used for:
A) the V array dimension is the number of the non-zero elements of the matrix A, and the non-zero elements of the matrix A are stored in a row traversal mode from top to bottom and from left to right;
B) the dimension of the C array is the same as that of the V array and is used for storing the column index of the element in the V array;
C) the dimension of the R array is the row number of the matrix A, and the index of the first non-zero element in each row of the matrix A in V is stored;
in step S3, obtaining an expression of an equation iterative solution according to a Krylov subspace method, and performing solution operation by adopting a sparse matrix; when the solved object matrix is large, decomposing the calculation load according to the hardware resources in the local area network, and distributing the calculation load to each process of the local area network;
s4, in each process, according to the task distributed by the main process, corresponding calculation is carried out;
s5, collecting the data after the calculation of each process into a main process, and carrying out comprehensive processing on the data to obtain the current iteration solution vector;
s6, carrying out convergence judgment according to iteration results, entering the step S7 when the matrix solution is converged, otherwise returning to the step S3, updating the Jacobian matrix and the right-end item again according to the solution vector obtained by current calculation, and continuing to iteratively solve the current matrix until iteration is converged; when the iteration times exceed the maximum iteration times and still do not converge, warning information is given;
s7, comparing the state of the solution vector and the right-end term with a convergence criterion, judging whether the current nonlinear iteration solution convergence condition is satisfied, substituting the solution vector obtained by iteration into a Jacobian matrix and the right-end term when the convergence condition is not satisfied, updating equation coefficients and excitation elements again, and solving equations until the nonlinear solution convergence; and repeating the nonlinear iteration process until all the solving tasks are completed.
The invention provides an electrical simulation technology based on a distributed parallel operation method in order to improve the solving efficiency of the simulation of a complex electrical system. The method is based on a Krylov subspace iteration method to realize distributed parallel operation. By fully utilizing hardware computing resources in the local area network, the method can obviously improve the solving efficiency of the simulation of the complex large system. Because the main computational load of the Krylov subspace iteration method is matrix vector multiplication, the sparse matrix can be fully utilized to further improve the computational efficiency, and therefore in the method, the comprehensive simulation of the complex electrical system is realized by adopting a parallel computation algorithm based on the sparse matrix. The distributed parallel computing realizes multi-process parallel computing based on MPI parallel environment.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (5)

1. An electrical simulation method based on a distributed parallel operation method is characterized in that: the method specifically comprises the following steps:
s1, reading a netlist file by a solver during initialization, wherein the netlist file comprises detailed information of element types, element pin node connection information, element model parameters, simulation step lengths, convergence parameters and an integration method, and constructing a Jacobian initial matrix based on a node voltage method; when Newton-Raphson nonlinear iteration is carried out, updating the elements of the Jacobian matrix and the right-end term according to the current solution vector state at any iteration moment;
s2, preprocessing the Jacobian matrix to reduce the condition number of the coefficient matrix;
s3, solving a matrix equation by adopting a Krylov subspace method;
s4, in each process, according to the task distributed by the main process, corresponding calculation is carried out;
s5, collecting the data after the calculation of each process into a main process, and carrying out comprehensive processing on the data to obtain the current iteration solution vector;
s6, carrying out convergence judgment according to iteration results, entering the step S7 when the matrix solution is converged, otherwise returning to the step S3, updating the Jacobian matrix and the right-end item again according to the solution vector obtained by current calculation, and continuing to iteratively solve the current matrix until iteration is converged; when the iteration times exceed the maximum iteration times and still do not converge, warning information is given;
s7, comparing the state of the solution vector and the right-end term with a convergence criterion, judging whether the current nonlinear iteration solution convergence condition is satisfied, substituting the solution vector obtained by iteration into a Jacobian matrix and the right-end term when the convergence condition is not satisfied, updating equation coefficients and excitation elements again, and solving equations until the nonlinear solution convergence; and repeating the nonlinear iteration process until all the solving tasks are completed.
2. The electrical simulation method based on the distributed parallel operation method according to claim 1, wherein: the initial state of the solution vector in step S1 is obtained from the historical solution data before entering the iteration, and after entering the nonlinear iteration, the solution vector is obtained for each iteration.
3. The electrical simulation method based on the distributed parallel operation method according to claim 1, wherein: in the newton-raphson nonlinear iterative calculation process in step S1, a step-changing algorithm is employed.
4. The electrical simulation method based on the distributed parallel operation method according to claim 1, wherein: after the preprocessing algorithm is adopted in step S2, the right-end item is synchronously processed according to the selected preprocessor.
5. The electrical simulation method based on the distributed parallel operation method according to claim 1, wherein: in step S3, obtaining an expression of an equation iterative solution according to a Krylov subspace method, and performing solution operation by adopting a sparse matrix; when the solved object matrix is large, the calculation load is decomposed according to the hardware resources in the local area network and is distributed to each process of the local area network.
CN202110302150.8A 2021-03-22 2021-03-22 Electrical simulation method based on distributed parallel operation method Pending CN113191105A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110302150.8A CN113191105A (en) 2021-03-22 2021-03-22 Electrical simulation method based on distributed parallel operation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110302150.8A CN113191105A (en) 2021-03-22 2021-03-22 Electrical simulation method based on distributed parallel operation method

Publications (1)

Publication Number Publication Date
CN113191105A true CN113191105A (en) 2021-07-30

Family

ID=76973568

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110302150.8A Pending CN113191105A (en) 2021-03-22 2021-03-22 Electrical simulation method based on distributed parallel operation method

Country Status (1)

Country Link
CN (1) CN113191105A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722966A (en) * 2021-09-07 2021-11-30 芯和半导体科技(上海)有限公司 Integrated circuit board simulation multistage distributed parallel computing method
CN114329981A (en) * 2021-12-29 2022-04-12 西安交通大学 DGMRES double-iteration high-dimensional electromagnetic simulation method and system based on ADMM preconditions
CN115718986A (en) * 2022-10-31 2023-02-28 南方电网数字电网研究院有限公司 Multi-core parallel time domain simulation method based on distributed memory architecture

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103617150A (en) * 2013-11-19 2014-03-05 国家电网公司 GPU (graphic processing unit) based parallel power flow calculation system and method for large-scale power system
US20140114636A1 (en) * 2012-10-18 2014-04-24 Nanotropic S.A. Large-scale power grid analysis on parallel architectures
CN104182209A (en) * 2014-08-27 2014-12-03 中国科学院软件研究所 PETSc-based GCRO-DR algorithm parallel processing method
CN105045768A (en) * 2015-09-01 2015-11-11 浪潮(北京)电子信息产业有限公司 Method and system for achieving GMRES algorithm
CN105260342A (en) * 2015-09-22 2016-01-20 浪潮(北京)电子信息产业有限公司 Solving method and system for symmetric positive definite linear equation set
CN110649624A (en) * 2019-11-09 2020-01-03 国网吉林省电力有限公司 Power flow parallel computing method for electric power system
US10860767B1 (en) * 2018-01-02 2020-12-08 Cadence Design Systems, Inc. Systems and methods for transient simulation of circuits with mutual inductors

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140114636A1 (en) * 2012-10-18 2014-04-24 Nanotropic S.A. Large-scale power grid analysis on parallel architectures
CN103617150A (en) * 2013-11-19 2014-03-05 国家电网公司 GPU (graphic processing unit) based parallel power flow calculation system and method for large-scale power system
CN104182209A (en) * 2014-08-27 2014-12-03 中国科学院软件研究所 PETSc-based GCRO-DR algorithm parallel processing method
CN105045768A (en) * 2015-09-01 2015-11-11 浪潮(北京)电子信息产业有限公司 Method and system for achieving GMRES algorithm
CN105260342A (en) * 2015-09-22 2016-01-20 浪潮(北京)电子信息产业有限公司 Solving method and system for symmetric positive definite linear equation set
US10860767B1 (en) * 2018-01-02 2020-12-08 Cadence Design Systems, Inc. Systems and methods for transient simulation of circuits with mutual inductors
CN110649624A (en) * 2019-11-09 2020-01-03 国网吉林省电力有限公司 Power flow parallel computing method for electric power system

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
FENG TU,ETC.: "A Message-Passing Distributed-Memory Newton-GMRES Parallel Power Flow Algorithm", 《IEEE POWER ENGINEERING SUMMER MEETING》, pages 1477 - 1482 *
R. SREERAMA KUMAR.ETC.: "A parallel distributed computing framework for Newton–Raphson load flow analysis of large interconnected power systems", 《ELECTRICAL POWER AND ENERGY SYSTEMS》, pages 1 - 6 *
REIJER IDEMA,ETC.: "Scalable Newton-Krylov Solver for Very Large Power Flow Problems", 《IEEE TRANSACTIONS ON POWER SYSTEMS》, pages 390 - 396 *
唐坤杰;董树锋;宋永华;: "基于不完全LU分解预处理迭代法的电力系统潮流算法", 中国电机工程学报, no. 1, pages 55 - 62 *
张志立;吴晓文;李自品;: "基于可变预条件处理与Broyden修正技术的不精确牛顿法潮流计算", 武汉大学学报(工学版), no. 04, pages 530 - 533 *
王明清;李明;张清;张广勇;吴韶华;: "基于MIC集群平台的GMRES算法并行加速", 计算机科学, no. 04, pages 197 - 201 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113722966A (en) * 2021-09-07 2021-11-30 芯和半导体科技(上海)有限公司 Integrated circuit board simulation multistage distributed parallel computing method
CN114329981A (en) * 2021-12-29 2022-04-12 西安交通大学 DGMRES double-iteration high-dimensional electromagnetic simulation method and system based on ADMM preconditions
CN114329981B (en) * 2021-12-29 2024-04-02 西安交通大学 DGMRES double-iteration high-dimensional electromagnetic simulation method and system based on ADMM precondition
CN115718986A (en) * 2022-10-31 2023-02-28 南方电网数字电网研究院有限公司 Multi-core parallel time domain simulation method based on distributed memory architecture
CN115718986B (en) * 2022-10-31 2023-12-12 南方电网数字电网研究院有限公司 Multi-core parallel time domain simulation method based on distributed memory architecture

Similar Documents

Publication Publication Date Title
CN113191105A (en) Electrical simulation method based on distributed parallel operation method
Chen et al. A Jacobian-free Newton-GMRES (m) method with adaptive preconditioner and its application for power flow calculations
Kedward et al. Efficient and exact mesh deformation using multiscale RBF interpolation
Rendall et al. Parallel efficient mesh motion using radial basis functions with application to multi‐bladed rotors
CN114091363B (en) Quantum algorithm-based computational fluid dynamics simulation method, device and equipment
Wang et al. Delaunay graph and radial basis function for fast quality mesh deformation
Wang et al. PR-ELM: Parallel regularized extreme learning machine based on cluster
Zhao et al. Power grid analysis with hierarchical support graphs
US8533653B2 (en) Support apparatus and method for simplifying design parameters during a simulation process
US20170140072A1 (en) Method and system for determining a configuration of a model having a collection of entities and satisfying a set of constraints
Wang et al. An efficient maximum bound principle preserving p-adaptive operator-splitting method for three-dimensional phase field shape transformation model
Zhang et al. A block-diagonal structured model reduction scheme for power grid networks
CN114880971B (en) Implicit method adopted by computational fluid dynamics software
Guo et al. Improved approximate minimum degree ordering method and its application for electrical power network analysis and computation
Zhang et al. A fast direct solver for boundary value problems on locally perturbed geometries
Rai Three-dimensional aerodynamic design using artificial neural networks
Georgieva et al. An algorithm for low-rank approximation of bivariate functions using splines
CN112949238A (en) Electrical simulation method based on iterative method
Nabors et al. Multipole and precorrected-FFT accelerated iterative methods for solving surface integral formulations of three-dimensional Laplace problems
Richardson et al. Regularized spectral multipole BEM for plane elasticity
Marco et al. Parallel genetic algorithms applied to optimum shape design in aeronautics
Wu et al. Adjoint-based optimization for blended-wing-body underwater gliders' shape design
Badcock et al. The application of sparse matrix techniques for the CFD based aeroelastic bifurcation analysis of a symmetric aerofoil
Yetkin et al. Recycling Newton–Krylov algorithm for efficient solution of large scale power systems
Martin-Burgos et al. A novel surface mesh deformation method for handling wing-fuselage intersections

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination