CN116050196B

CN116050196B - Multi-dimensional simulation method, device, equipment and storage medium

Info

Publication number: CN116050196B
Application number: CN202310345970.4A
Authority: CN
Inventors: 杨灿群; 黄颖杰; 王伟; 仲彦旭; 郑伟龙; 卢海林; 夏梓峻; 段莉莉
Original assignee: Haihe Laboratory Of Advanced Computing And Key Software Xinchuang; National Supercomputer Center In Tianjin; National University of Defense Technology
Current assignee: Haihe Laboratory Of Advanced Computing And Key Software Xinchuang; National Supercomputer Center In Tianjin; National University of Defense Technology
Priority date: 2023-04-03
Filing date: 2023-04-03
Publication date: 2023-06-30
Anticipated expiration: 2043-04-03
Also published as: CN116050196A

Abstract

The embodiment of the disclosure relates to a multi-dimensional simulation method, a device, equipment and a storage medium, wherein the method comprises the steps of carrying out simulation initialization aiming at a target engineering problem; setting intermittent monitoring points in the low-dimensional simulation model and determining the positions of the intermittent monitoring points; aiming at the low-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of a substance in the low-dimensional simulation model until the current position of the target discontinuity and the position of a discontinuity monitoring point meet preset conditions; high-dimensional reconstruction is carried out on the low-dimensional simulation model and the high-dimensional reconstruction is mapped into the high-dimensional simulation model, so that a new high-dimensional simulation model is obtained; aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration. According to the embodiment of the disclosure, the calculated amount can be reduced, the calculation efficiency is improved, the calculation error can be reduced, and the usability of multi-dimensional simulation is improved.

Description

Multi-dimensional simulation method, device, equipment and storage medium

Technical Field

The embodiment of the disclosure relates to the technical field of CAE simulation, in particular to a multi-dimensional simulation method, a device, equipment and a storage medium.

Background

In the process of equipment research and development and design in various fields, the problems of long period, large risk coefficient, high cost and the like exist in research and development through a test method, and the CAE simulation technology can effectively replace a test, greatly shorten the research and development period and reduce the research and development cost, so that the CAE simulation technology plays an increasingly important role in the process of equipment research and development.

At present, the simulation process for actual engineering problems is generally as follows: and establishing a high-dimensional simulation model corresponding to the actual engineering problem, and then simulating the dynamic interaction based on the high-dimensional simulation model. However, for some engineering problems (such as energy-accumulating jet penetration, underwater/aerial explosion, etc.), the calculation results of a large number of areas are not practical, because the attention content is less (such as only the dynamic mechanical response process of a partial area is concerned), and the partial model can be simplified (such as axisymmetry, plane symmetry, etc.), so that the problems of large calculation amount and low resource utilization rate exist, and a brake elbow is formed for engineering application.

Disclosure of Invention

In order to solve the above technical problems or at least partially solve the above technical problems, embodiments of the present disclosure provide a multi-dimensional simulation method, apparatus, device, and storage medium.

A first aspect of an embodiment of the present disclosure provides a multi-dimensional simulation method, including:

performing simulation initialization for the target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem;

setting intermittent monitoring points for a target substance in a low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substance interacts with a target intermittent, and the target intermittent comprises a contact intermittent and/or a shock wave intermittent;

aiming at the low-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of a substance in the low-dimensional simulation model until the current position of the target discontinuity and the position of a discontinuity monitoring point meet preset conditions;

high-dimensional reconstruction is carried out on the low-dimensional simulation model and the high-dimensional reconstruction is mapped into the high-dimensional simulation model, so that a new high-dimensional simulation model is obtained;

aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

A second aspect of an embodiment of the present disclosure provides a multi-dimensional simulation apparatus, the apparatus comprising:

The initialization module is used for carrying out simulation initialization aiming at the target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem;

the setting and determining module is used for setting intermittent monitoring points for target substances in the low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substances interact with the target intermittent, and the target intermittent comprises a contact intermittent and/or a shock wave intermittent;

the first determining module is used for determining the current position of the target discontinuity and the current physical state of a substance in the low-dimensional simulation model aiming at the low-dimensional simulation model until the current position of the target discontinuity and the position of a discontinuity monitoring point meet preset conditions;

the reconstruction module is used for carrying out high-dimensional reconstruction on the low-dimensional simulation model and mapping the high-dimensional reconstruction to the high-dimensional simulation model to obtain a new high-dimensional simulation model;

the second determining module is used for determining the current position of the target break and determining the current physical state of the substances in the new high-dimensional simulation model aiming at the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

A third aspect of the disclosed embodiments provides an electronic device, the server comprising: a processor and a memory, wherein the memory has stored therein a computer program which, when executed by the processor, performs the method of the first aspect described above.

A fourth aspect of the disclosed embodiments provides a computer readable storage medium having stored therein a computer program which, when executed by a processor, can implement the method of the first aspect described above.

Compared with the prior art, the technical scheme provided by the embodiment of the disclosure has the following advantages:

according to the embodiment of the disclosure, simulation initialization can be performed aiming at the target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem; setting intermittent monitoring points for a target substance in a low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substance interacts with a target intermittent, and the target intermittent comprises a contact intermittent and/or a shock wave intermittent; aiming at the low-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of a substance in the low-dimensional simulation model until the current position of the target discontinuity and the position of a discontinuity monitoring point meet preset conditions; high-dimensional reconstruction is carried out on the low-dimensional simulation model and the high-dimensional reconstruction is mapped into the high-dimensional simulation model, so that a new high-dimensional simulation model is obtained; aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration. Therefore, by adopting the technical scheme, part of problems (namely in an initial period of time) can be calculated in the low-dimensional simulation model, so that the high-dimensional problems are simplified into the low-dimensional problems to be calculated, and the high-dimensional simulation model is reconstructed in a high-dimensional mode and mapped into the high-dimensional simulation model to be continuously solved. And moreover, the accurate capture of the target discontinuities such as shock wave discontinuities, contact discontinuities and the like in a calculation region can be realized by setting the discontinuity monitoring points, so that the accurate control of the conversion opportunity among different dimension models is realized, the low-dimension simulation model is automatically mapped to the high-dimension simulation model at a proper opportunity, the calculation error is reduced, and the usability of multi-dimension simulation is improved.

Drawings

The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description, serve to explain the principles of the disclosure.

In order to more clearly illustrate the embodiments of the present disclosure or the solutions in the prior art, the drawings that are required for the description of the embodiments or the prior art will be briefly described below, and it will be obvious to those skilled in the art that other drawings can be obtained from these drawings without inventive effort.

FIG. 1 is a flow chart of a multi-dimensional simulation method provided by an embodiment of the present disclosure;

FIG. 2 is a schematic cross-sectional view of a high-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the present disclosure;

FIG. 3 is a schematic view of a partial section of a low-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the present disclosure;

FIG. 4 is a schematic illustration of a setup of intermittent monitoring points provided by an embodiment of the present disclosure;

FIG. 5 is a schematic diagram of high-dimensional reconstruction and mapping of a low-dimensional material model provided by embodiments of the present disclosure;

FIG. 6 is a schematic diagram of decoupling a multi-media interaction problem provided by an embodiment of the present disclosure;

FIG. 7 is a flow diagram of a multi-dimensional simulation process provided by an embodiment of the present disclosure;

fig. 8 is a schematic diagram of calculation results corresponding to a time, b time, c time and d time of a low-dimensional simulation model according to an embodiment of the present disclosure;

FIG. 9 is a schematic cross-sectional view of a new high-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the present invention;

FIG. 10 is a schematic diagram of the calculation results of a new high-dimensional simulation model provided by an embodiment of the present disclosure;

FIG. 11 is a simulated computation time contrast diagram provided by an embodiment of the present disclosure;

FIG. 12 is a schematic diagram of a multi-dimensional simulation device according to an embodiment of the present disclosure;

fig. 13 is a schematic structural view of an electronic device in an embodiment of the present disclosure.

Detailed Description

In order that the above objects, features and advantages of the present disclosure may be more clearly understood, a further description of aspects of the present disclosure will be provided below. It should be noted that, without conflict, the embodiments of the present disclosure and features in the embodiments may be combined with each other.

In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present disclosure, but the present disclosure may be practiced otherwise than as described herein; it will be apparent that the embodiments in the specification are only some, but not all, embodiments of the disclosure.

Fig. 1 is a flow chart of a multi-dimensional simulation method provided by an embodiment of the present disclosure, which may be performed by an electronic device. The electronic device may be exemplarily understood as a device having a page presentation function, such as a mobile phone, a tablet computer, a notebook computer, a desktop computer, a smart television, etc. As shown in fig. 1, the method provided in this embodiment includes the following steps:

s110, carrying out simulation initialization on the target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem.

In the embodiment of the disclosure, the simulation of the target engineering problem needs to be initiated by first performing simulation on the target engineering problem, so that preparation is made for subsequent simulation calculation of interaction of substances in the target engineering problem.

Specifically, the target engineering problems may include, but are not limited to, engineering problems in the field of energy-gathering jet damage and protection, engineering problems in the field of underwater/aerial explosion damage and protection, engineering problems in the field of hypersonic advanced weapon development and design, engineering problems in the field of armor protection design and evaluation, engineering problems in the field of hazardous chemicals production/storage/use, engineering problems in the field of ship design optimization, engineering problems in the field of bubble dynamics, engineering problems in the field of aerospace, engineering problems in the field of risk assessment and early warning, and the like.

Optionally, the simulation initialization may include establishing an euler coordinate system, establishing a calculation region and grid division (corresponding to S111), establishing a model (corresponding to S112), defining a target discontinuity (corresponding to S113), setting initial conditions (corresponding to S114), selecting a control equation (corresponding to S115), and setting a solution time (corresponding to S116). Accordingly, S110 may include S111-S116 as follows:

s111, establishing an Euler coordinate system according to the target engineering problem, defining the size of a calculation region, and carrying out grid division on the calculation region according to the preset Euler grid cell size.

The specific value of the preset euler mesh unit size may be set according to practical situations, and is not limited herein.

S112, establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem.

Optionally, S112 may include: establishing a corresponding high-dimensional material model aiming at a material in a target engineering problem so as to obtain a high-dimensional simulation model, wherein the high-dimensional material model corresponding to the material with a fluid phase is a Gao Weiou Law model, and the high-dimensional material model corresponding to the material with a solid phase is a Gao Weiou Law model or a high-dimensional Lagrange model; and performing dimension reduction modeling on the high-dimensional material model with the spherical symmetry and/or the axial symmetry to obtain a corresponding low-dimensional material model, thereby obtaining a low-dimensional simulation model.

Specifically, the high-dimensional simulation model includes a high-dimensional material model corresponding to a plurality of materials in the target engineering problem, and the high-dimensional material model may include a three-dimensional model. The low-dimensional simulation model comprises a low-dimensional material model corresponding to a plurality of materials in the target engineering problem, and the low-dimensional material model can comprise a one-dimensional model and/or a two-dimensional model.

Specifically, a gas-liquid-solid full-size three-dimensional model is built for the whole target engineering problem, a Gao Weili Euler model (i.e., a high-dimensional substance model) is built for a substance whose phase is fluid, and a high-dimensional Lagrange model or a high-dimensional Euler model (i.e., a high-dimensional substance model) is built for a substance whose phase is solid. And according to the spherical symmetry or axisymmetric characteristic of the high-dimensional material model, performing dimension reduction modeling, establishing a one-dimensional model (namely a low-dimensional material model) aiming at the spherical symmetry high-dimensional simulation model, and establishing a two-dimensional model (namely a low-dimensional material model) aiming at the axisymmetry high-dimensional material model.

Exemplary, fig. 2 is a schematic cross-sectional view of a high-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the disclosure, and fig. 3 is a schematic cross-sectional view of a low-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the disclosure. Referring to fig. 2 and 3, the engineering problem of penetration of the focused jet into the steel target plate includes the following: explosive 210, liner 220, target plate 230, and air 240, the high-dimensional simulation model includes high-dimensional material models corresponding to explosive 210, liner 220, target plate 230, and air 240, and the low-dimensional simulation model includes low-dimensional material models corresponding to explosive 210, liner 220, and air 240.

It should be noted that, the low-dimensional simulation model includes a low-dimensional material model corresponding to which materials in the target engineering problem, which can be determined according to the actual engineering problem, and is not limited herein.

It can be understood that the above-mentioned method for establishing the high-dimensional simulation model and the low-dimensional simulation model is simple and easy to implement, and establishes a one-dimensional low-dimensional material model according to the spherical symmetry characteristics, and establishes a two-dimensional low-dimensional material model according to the axisymmetric characteristics, so that the dimension of the low-dimensional material model is as low as possible, which is beneficial to reducing the calculation amount.

And S113, when the target discontinuity is a contact discontinuity (namely a multi-medium interface), defining an initial position of the multi-medium interface corresponding to the fluid by adopting a multi-medium interface tracking method.

Alternatively, an interface tracking method based on a level set method may be used to determine the initial position of the multi-media interface according to a symbolic distance function.

Alternatively, when the target discontinuity is a shockwave discontinuity, a shockwave discontinuity indicator may be used to determine an initial location of the shockwave discontinuity.

S114, setting initial physical states of substances in the high-dimensional simulation model and the low-dimensional simulation model.

Optionally, the physical state may include at least one of material properties, speed, pressure, material model, boundary conditions, detonation point location, detonation type, and the like.

S115, selecting a proper control equation (namely a preset low-dimensional control equation and a preset high-dimensional control equation) for the target engineering problem.

Alternatively, the control equation may include an N-S equation, a non-stick euler equation, an elastoplastic euler equation that considers the strength of the solid, and the like.

S116, setting calculation termination time (namely preset termination time).

S120, setting intermittent monitoring points for a target substance in the low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substance interacts with the target intermittent, and the target intermittent comprises a contact intermittent and/or a shock wave intermittent.

In the embodiment of the disclosure, in the process of simulating the target engineering problem, different parts of the model are distinguished according to whether the dimension reduction simplification can be performed, wherein the parts which cannot be subjected to the dimension reduction, namely the target substance, are not subjected to the dimension reduction. When the target substance and the target discontinuity are close to each other, the target substance and the target discontinuity are about to interact with each other, so that a discontinuity monitoring point can be set for the target substance and the position of the discontinuity monitoring point can be determined, whether the interaction between the target substance and the target discontinuity is about to occur or not can be judged according to the position of the target discontinuity and the position of the discontinuity monitoring point, and the control is accurately converted from the low-dimensional calculation to the high-dimensional calculation.

Optionally, setting intermittent monitoring points for the target substance in the low-dimensional simulation model may include: in the low-dimensional simulation model, intermittent monitoring points are set at positions spaced apart from a target substance by a preset number of Euler grid cells. Therefore, the monitoring discontinuity point can be arranged at a position close to the target substance, so that the distance between the target discontinuity and the discontinuity monitoring point is approximately equal to the distance between the target discontinuity and the target substance, the distance between the target discontinuity and the target substance can be judged according to the position of the target discontinuity and the position of the discontinuity monitoring point, the accuracy of simulation calculation is ensured, and whether interaction between the target substance and the target discontinuity is about to occur is judged, so that the control is switched from low-dimensional calculation to high-dimensional calculation. In addition, the monitoring discontinuity point is not a substance in the target engineering problem, is only used for monitoring the degree of the target discontinuity distance from the target substance and does not interact with other substances in the target engineering problem, so that the degree of the target discontinuity approaching the target substance is judged to be more accurate by taking the position of the monitoring discontinuity point as a reference, and the time for accurately converting control from low-dimensional calculation to high-dimensional calculation is facilitated.

Specifically, the specific values of the preset number may be set by those skilled in the art according to actual circumstances, and are not limited herein. For example, the preset number is 2, 3, 4, etc., but is not limited thereto.

Alternatively, intermittent monitoring points may be provided at positions spaced apart from the target substance by a preset number of euler mesh cells at the periphery of the target substance.

Further optionally, the number of intermittent monitoring points is multiple, and the multiple intermittent monitoring points are at least arranged on a side of the target substance facing the target interruption. Therefore, a plurality of intermittent monitoring points can be arranged on the path of the target intermittent moving towards the target substance, which is beneficial to more accurately monitoring the target intermittent.

Further optionally, a plurality of intermittent monitoring points encircle the target substance. Thus, the method is beneficial to monitor target interruption in multiple directions and more accurately.

Specifically, the specific number of intermittent monitoring points can be set by those skilled in the art according to actual situations, and is not limited herein.

Specifically, the distance between adjacent intermittent monitoring points may be equal to 0 (i.e., no gap between adjacent intermittent monitoring points), or may be greater than 0 (i.e., no gap between adjacent intermittent monitoring points, discontinuity), which is not limited herein.

Fig. 4 is a schematic diagram illustrating a setting of intermittent monitoring points according to an embodiment of the present disclosure. Referring to fig. 4, for a shock wave break or contact break, break monitoring points are set at a grid of the peripheral n (n=3 shown in fig. 4) layers where the target substance is located. According to the actual problem type, the intermittent monitoring points can be shock wave intermittent monitoring points used for monitoring and identifying shock wave intermittent or contact intermittent monitoring points used for monitoring and identifying contact intermittent.

Illustratively, with continued reference to fig. 2 and 3, the target material is a target plate 230 and the target discontinuity is a multi-media interface between the explosive 210 and liner 220, and therefore, a plurality of discontinuity monitoring points 250 are provided on a side of the multi-media interface facing the target plate 230.

Specifically, after the intermittent monitoring points are selected, the positions of the intermittent monitoring points in the low-dimensional simulation model can be determined, and details are omitted here.

S130, aiming at the low-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substance in the low-dimensional simulation model until the current position of the target discontinuity and the position of the discontinuity monitoring point meet preset conditions.

And S140, performing high-dimensional reconstruction on the low-dimensional simulation model and mapping the high-dimensional reconstruction to the high-dimensional simulation model to obtain a new high-dimensional simulation model.

In the embodiment of the present disclosure, in order to reduce the amount of calculation in the process of simulating the target engineering problem, when the target substance and the target discontinuity are far apart (i.e., the target substance and the target discontinuity are not about to interact with each other), the calculation may be performed in a low dimension, and specifically, the following steps may be periodically performed until the current position of the target discontinuity and the position of the discontinuity monitoring point satisfy the preset condition: in the current computing period, for the low-dimensional simulation model, determining a current position of the target discontinuity and determining a current physical state of the substance in the low-dimensional simulation model. When the current position of the target discontinuity and the position of the discontinuity monitoring point meet the preset condition, the interaction between the target discontinuity and the target substance is indicated, and at the moment, high-dimensional reconstruction and mapping can be performed on the low-dimensional simulation model so as to perform calculation in high dimensions.

Optionally, the preset condition is that a distance between the target discontinuity and the discontinuity monitoring point is smaller than a preset distance, wherein the distance between the target discontinuity and the discontinuity monitoring point is determined by a current position of the target discontinuity and a position of the discontinuity monitoring point.

Specifically, the specific value of the preset distance may be set by those skilled in the art according to the actual situation, and is not limited herein.

Optionally, the preset distance is a preset euler mesh cell size.

In some embodiments, the target discontinuity is a contact discontinuity, and the discontinuity monitoring point is a contact discontinuity monitoring point, at this time, a distance between the contact discontinuity monitoring point and the contact discontinuity can be determined according to a current position of the contact discontinuity and a position of the contact discontinuity monitoring point, so as to determine whether high-dimensional reconstruction and mapping are needed. Specifically, S1311, determining, for the low-dimensional simulation model, a current position of the contact discontinuity and determining a current physical state of a substance in the low-dimensional simulation model; s1312, judging whether the distance between the contact discontinuity and the contact discontinuity monitoring point is smaller than the preset Euler mesh unit size, if yes, executing S140, and if not, returning to execute the step S1311 until the distance between the contact discontinuity and the contact discontinuity monitoring point is determined to be smaller than the preset Euler mesh unit size according to the current position of the contact discontinuity and the position of the contact discontinuity monitoring point.

In particular, the time advance step Δt can be determined according to the hyperbolic conservation law stabilization condition, the understanding about the time advance step being as follows: if the last calculation period calculates the target interruption at t _i-1 The position of the target break calculated in the current calculation period (i.e. the current position) is the position of the target break from t _i-1 Starting at moment, the movement time advances by a step delta t and then at t _i Time (t) _i = t _i-1 A + [ delta ] t) the location reached; if the last calculation period calculates to obtain the substance in the low-dimensional simulation model at t _i-1 The physical state of the material in the low-dimensional simulation model calculated in the current calculation period (i.e. the current physical state) is the physical state of the material in the low-dimensional simulation model from t _i-1 Starting at the moment, the time advances by a step delta t and then at t _i Physical state of time of day. Since the contact discontinuity moves a distance per time step less than a preset Euler mesh cell size, a high-dimensional reconstruction and mapping bar can be determinedThe parts are as follows:

wherein, the liquid crystal display device comprises a liquid crystal display device,

for the distance between the contact discontinuity monitoring point and the contact discontinuity, < > j->

Is a preset euler mesh cell size. If the distance between the contact intermittent monitoring point and the contact intermittent is smaller than the size of a preset Euler grid unit according to the current position of the contact intermittent and the position of the contact intermittent monitoring point calculated in the current calculation period, high-dimensional reconstruction and mapping are carried out, and otherwise, the next calculation period is entered.

In other embodiments, the target discontinuity is a shockwave discontinuity, and the discontinuity monitoring point is a shockwave discontinuity monitoring point, and at this time, the distance between the shockwave discontinuity monitoring point and the shockwave discontinuity can be determined according to the current position of the shockwave discontinuity and the position contacting the discontinuity monitoring point, so as to determine whether high-dimensional reconstruction and mapping are needed. Specifically, S1321, determining, for the low-dimensional simulation model, a current position of the shock wave discontinuity and determining a current physical state of a substance in the low-dimensional simulation model; s1322, judging whether the distance between the shock wave interruption and the shock wave interruption monitoring point is smaller than the preset Euler mesh unit size, if yes, executing S140, and if not, returning to the step S1321 until the distance between the shock wave interruption and the shock wave interruption monitoring point is smaller than the preset Euler mesh unit size according to the current position of the shock wave interruption and the position of the shock wave interruption monitoring point.

Specifically, the conditions for high-dimensional reconstruction and mapping are as follows:

for the distance between the shock wave break monitoring point and the shock wave break, < >>

Is a preset euler mesh cell size. If the distance between the shock wave interruption monitoring point and the shock wave interruption is smaller than the preset Euler grid cell size according to the current position of the shock wave interruption and the position of the shock wave interruption monitoring point obtained by calculation in the current calculation period, high-dimensional reconstruction and mapping are carried out, and otherwise, the next calculation period is entered.

Optionally, S140 may include: performing high-dimensional reconstruction on the low-dimensional simulation model to obtain a Gao Weichong constructed model; and mapping the Gao Weichong structural model to an Euler grid in the high-dimensional simulation model to obtain a new high-dimensional simulation model.

Specifically, the low-dimensional data of the low-dimensional simulation model is subjected to high-dimensional reconstruction in different directions to obtain high-dimensional data corresponding to the high-dimensional reconstruction model, for example, for the low-dimensional simulation model, if the high-dimensional material model corresponding to the low-dimensional material model is in spherical symmetry, the low-dimensional material model is reconstructed into a three-dimensional model with spherical symmetry, and if the high-dimensional material model corresponding to the low-dimensional material model is in axial symmetry, the low-dimensional material model is reconstructed into a three-dimensional model with axial symmetry, so that the high-dimensional reconstruction model can be obtained. Because the high-dimensional data corresponding to the Gao Weichong constructed model is data which does not contain a grid structure, the high-dimensional data obtained through reconstruction can be mapped to a high-dimensional Euler grid in the high-dimensional simulation model, and a new high-dimensional simulation model containing a low-dimensional calculation result is obtained.

Illustratively, FIG. 5 is a schematic diagram of high-dimensional reconstruction and mapping of a low-dimensional material model provided by embodiments of the present disclosure. Referring to fig. 5, a high-dimensional reconstruction is performed on a one-dimensional low-dimensional substance model to obtain a corresponding ball-symmetrical three-dimensional model (only a hemisphere is shown in fig. 5 for convenience of drawing), and then the three-dimensional model is mapped into a high-dimensional euler mesh.

It can be understood that the embodiment of the disclosure can automatically reconstruct the low-dimensional data in different directions to obtain the high-dimensional data, map the high-dimensional data to Euler grids in the high-dimensional simulation model, ensure accurate matching between the grids in different dimensions and the state data, and reduce calculation errors.

S150, aiming at the new high-dimensional simulation model, determining the current position of the target break and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

In the embodiment of the disclosure, as described above, in the process of simulating the target engineering problem, from the point that the interaction between the target substance and the target discontinuity is about to start, the calculation may be performed in a high dimension, and specifically, the following steps may be periodically performed until the simulation duration reaches the preset termination duration: in the current computing period, aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model. And when the simulation duration reaches the preset termination duration, the simulation is ended.

Specifically, when the simulation duration reaches the preset termination duration, the simulation motion of the target is started from the position where the target is positioned at the time of initialization, and after the simulation motion reaches the preset termination duration, the whole multi-dimensional simulation process is ended.

Specifically, the understanding of the "current position" and the "current physical state" for the new high-dimensional simulation model is similar to that for the new high-dimensional simulation model, and is not repeated here.

It can be appreciated that when the interval between the target substance and the target is far, the embodiment of the disclosure can convert the three-dimensional model into a one-dimensional, two-dimensional and other low-dimensional model for calculation, compared with the traditional full three-dimensional calculation, the method greatly saves calculation time, improves calculation efficiency, saves a large amount of calculation resources, avoids the situation of wasting a large amount of calculation resources in the three-dimensional calculation, and has obvious advantages in the engineering field. Moreover, the embodiment of the disclosure can realize accurate capture of shock wave interruption, contact interruption and the like in a calculation area based on the arrangement of the interruption monitoring points, and further realize accurate control of conversion opportunities among different dimension models. In addition, the embodiment of the disclosure can automatically reconstruct the low-dimensional data in different directions to obtain the high-dimensional data, map the high-dimensional data to the high-dimensional calculation grid, ensure accurate matching between the grids in different dimensions and the state data, reduce calculation errors and promote usability of multi-dimensional simulation.

According to the embodiment of the disclosure, part of problems (namely, in an initial period of time) can be calculated in the low-dimensional simulation model, so that the high-dimensional problems are simplified into the low-dimensional problems to be calculated, and the low-dimensional simulation model is reconstructed in a high-dimensional mode and mapped into the high-dimensional simulation model to be solved continuously. And moreover, the accurate capture of the target discontinuities such as shock wave discontinuities, contact discontinuities and the like in a calculation region can be realized by setting the discontinuity monitoring points, so that the accurate control of the conversion opportunity among different dimension models is realized, the low-dimension simulation model is automatically mapped to the high-dimension simulation model at a proper opportunity, the calculation error is reduced, and the usability of multi-dimension simulation is improved.

In another embodiment of the present disclosure, the target discontinuity is a contact discontinuity; wherein determining, for the low-dimensional simulation model, a current location of the target discontinuity and determining a current physical state of the substance in the low-dimensional simulation model comprises: s131, after interaction of different substances in the low-dimensional simulation model is decoupled, solving a preset low-dimensional control equation to obtain a corresponding current physical state aiming at each substance and target discontinuity in the low-dimensional simulation model; s132, determining the current position of the target break according to the current physical state of the low-dimensional simulation model.

Specifically, S131 may include: s1301, adopting a multi-medium processing method to process different medium interactions in the low-dimensional simulation model.

Optionally, a modified virtual fluid method (MGFM) may be used to process the multi-media interaction physical state solution problem, decouple the multi-media interaction problem into a plurality of single-media problems for processing, avoid non-physical oscillation on both sides of the contact discontinuity, and reduce the calculation error.

Illustratively, FIG. 6 is a schematic diagram of decoupling a multi-media interaction problem provided by embodiments of the present disclosure. Referring to fig. 6, the interaction problem of medium 1 and medium 2 is decoupled.

S1302, solving a preset low-dimensional control equation, and realizing space calculation and time propulsion of the low-dimensional simulation model to obtain the current physical state of the substances in the low-dimensional simulation model.

Alternatively, a finite difference method may be used for discrete computation, where the spatial computation format uses a 5-order WENO format and the time-marching format uses a 3-order R-K format.

Specifically, S132 may include: and calculating a level set function, pushing interface motion in the low-dimensional simulation model, and determining the current position of the contact discontinuity. Wherein the level set function is as follows:

u is the velocity of the substance in the x-direction, v is the velocity of the substance in the y-direction, and t is the time, as a function of the symbolic distance. Alternatively, the space discrete and time propulsion can be performed by adopting the five-order HJ-WENO and the three-order TVD Runge-Kutta format respectively, and the current position of the contact break can be obtained according to the value of the symbol distance function.

Accordingly, for the new high-dimensional simulation model, determining the current location of the target discontinuity and determining the current physical state of the substance in the new high-dimensional simulation model includes: s151, after interaction of different substances in the new high-dimensional simulation model is decoupled, aiming at each substance and target interruption in the new high-dimensional simulation model, a preset high-dimensional control equation is adopted to solve to obtain a corresponding current physical state, wherein the current physical state is a physical state when a time propulsion step length is passed after the physical state is obtained relative to the last solution; s152, determining the current position of the target break according to the current physical state of the new high-dimensional simulation model, wherein the current position is the position when the time pushing step length is passed after the position is obtained by solving for the last time.

Specifically, S151 may include: s1511, adopting a multi-medium processing method to process the interaction of different media in the high-dimensional simulation model.

S1512, solving a preset high-dimensional control equation, and realizing space calculation and time propulsion of the high-dimensional model to obtain the current physical state of the substance in the high-dimensional simulation model.

Optionally, for the high-dimensional simulation model, the Euler model can adopt a finite difference method for discrete calculation, a space calculation format adopts a 5-order WENO format, and a time propulsion format adopts a 3-order R-K format; the Lagrange model can adopt a nonlinear explicit finite element method to perform discrete calculation, and the time propulsion format adopts a 2-order central difference format.

S1513, if the high-dimensional simulation model comprises a Lagrangian model, a coupled Lagrangian-Euler (CEL) processing method can be adopted to process the interaction between Euler and Lagrangian model in the calculation domain so as to correct the current physical state obtained by calculation in S1512.

Optionally, a virtual Euler-Lagrange (GEL) method is adopted to process interaction between the Euler model and the Lagrange model, so that effective strong coupling between the Euler model and the Lagrange model is realized, and the calculation accuracy is improved.

Specifically, S152 may include: and calculating a level set function, pushing the interface in the new high-dimensional simulation model to move, and determining the current position of the contact discontinuity. Wherein the level set function is as follows:

as a function of the symbolic distance, u is the velocity of the substance in the x-direction, v is the velocity of the substance in the y-direction, w is the velocity of the substance in the z-direction, and t is the time. Alternatively, five-order HJ-WENO and three-order TVD Runge-Kutta formats can be adopted respectivelyAnd obtaining the current position of the contact discontinuity according to the value of the symbol distance function by space dispersion and time promotion.

It can be appreciated that by adopting an accurate contact discontinuity tracking method, accurate capture of contact discontinuities in a calculation region can be realized, so that the effectiveness and accuracy of discontinuity monitoring can be ensured, and the method is the basis and guarantee of multidimensional adaptive conversion.

In yet another embodiment of the present disclosure, the target discontinuity is a shockwave discontinuity, wherein determining the current location of the target discontinuity comprises: a current location of the target discontinuity is determined using the shock wave discontinuity indicator.

Alternatively, a solution-based average global variance (Average total variation of the solution, ATV) discontinuity indicator may be used to effectively identify the shockwave discontinuity, resulting in the location of the shockwave discontinuity in the calculation region.

It can be appreciated that by adopting the accurate shock wave break indicator method, accurate capture of shock wave break in the calculation area can be realized, so that the effectiveness and accuracy of break monitoring can be ensured, and the method is the basis and guarantee of multi-dimensional self-adaptive conversion.

For better illustrating the objects and advantages of the embodiments of the present disclosure, a detailed description will be made of the multi-dimensional simulation method provided by the embodiments of the present disclosure based on a specific example.

Fig. 7 is a flow chart of a multi-dimensional simulation process provided by an embodiment of the present disclosure. Referring to fig. 7, the multi-dimensional simulation process includes the steps of:

s710, CEA simulation initialization, intermittent monitoring points setting and intermittent monitoring point position determination.

Specifically, CEA simulation initialization may be initialized as described above with reference to S111-S116.

Specifically, setting the intermittent monitoring point and determining the position of the intermittent monitoring point may refer to the content described in the foregoing S120, and will not be described herein.

Illustratively, aiming at the engineering problem of penetration of the energy-accumulating jet flow into the steel target plate, a corresponding high-dimensional simulation model and a low-dimensional simulation are establishedA true model is shown in fig. 2 and 3. Wherein the low-dimensional simulation model is an axisymmetric model, the target substance is a target plate 230, the target discontinuity is a contact discontinuity between the explosive 210 and the liner 220, and thus, the discontinuity monitoring point 250 is a contact discontinuity monitoring point, and the calculation termination time (i.e., t is set ^end ) Set to 250 mus.

S720, calculating and solving the low-dimensional simulation model.

Specifically, for a low-dimensional simulation model, a current location of a target discontinuity is determined and a current physical state of a substance in the low-dimensional simulation model is determined. The low-dimensional simulation model calculation may be solved with reference to the descriptions of S131 and S132, and will not be described here.

S730, monitoring different discontinuities by the discontinuity monitoring points.

S740, judging whether the intermittent monitoring point is intermittent. If not, the process returns to S720, and if yes, S750 is executed.

Specifically, it may be determined whether the intermittent monitoring point has an interruption according to the specific method described in S1312, which is not described herein.

S750, constructing a low-dimensional simulation model Gao Weichong.

S760, high-dimensional data are mapped and interpolated to a high-dimensional grid.

Specifically, the high-dimensional reconstruction and mapping may refer to the content described in the foregoing S140, and will not be described herein.

Exemplary, fig. 9 is a schematic cross-sectional view of a new high-dimensional simulation model corresponding to a focused jet penetration steel target plate according to an embodiment of the present invention. Referring to fig. 9, when there is a break at the time-break monitoring point with the simulation duration of 208.1 μs, the high-dimensional reconstruction and mapping conditions are satisfied, at this time, S750 is executed, the high-dimensional reconstruction is performed on the low-dimensional data corresponding to the low-dimensional simulation model, and the high-dimensional data of the high-dimensional reconstruction model obtained by the reconstruction is mapped to the high-dimensional euler mesh, which contains the new high-dimensional simulation model of the low-dimensional calculation information.

S770, calculating and solving a new high-dimensional simulation model.

Specifically, for a new high-dimensional simulation model, a current location of the target discontinuity is determined and a current physical state of the material in the low-dimensional simulation model is determined. The low-dimensional simulation model calculation may be solved with reference to the descriptions of S151 and S152, and will not be described here.

Exemplary, the simulation result of far-field interaction of the energy-gathering jet and the steel target plate is obtained by carrying out calculation solution on a high-dimensional model containing low-dimensional calculation information.

After S770, it is determined whether the simulation time period t is longer than the preset termination time period t ^end . If yes, the physical state (i.e., simulation result) is output, if not, the process returns to S770.

For example, fig. 8 is a schematic diagram of calculation results corresponding to time a, time b, time c and time d of a low-dimensional simulation model provided by the embodiment of the present disclosure, according to fig. 8, it can be seen that an explosive explosion extrudes a liner to form an energy-accumulating jet, a high-speed jet impacts air to form a front shock wave break, then the jet is further stretched, the jet head becomes blunt due to air resistance, a numerical simulation result accords with an actual engineering problem, and accurate simulation of interactions of various mediums such as explosives, jet, air and the like can be realized, and a contact break is clearly and accurately distinguished. Fig. 10 is a schematic diagram of a calculation result of a new high-dimensional simulation model provided by the embodiment of the present disclosure, according to fig. 10, it can be known that after jet flow is further stretched, the target plate is penetrated, and the target plate forms penetration holes, which accords with the actual situation, and proves the effectiveness of a multi-dimensional adaptive coupling solving technique of the present invention. Fig. 11 is a comparison diagram of simulation calculation time provided by the embodiment of the present disclosure, referring to fig. 11, when a conventional full three-dimensional calculation method is used 6422S, and when a multi-dimensional simulation method implemented by the present disclosure is used 1573S, it can be seen that the embodiment of the present disclosure can greatly shorten calculation time, promote calculation speed, and has obvious advantages compared with the conventional full three-dimensional calculation method, thereby proving good calculation performance of the method of the present invention.

According to the embodiment of the disclosure, based on discontinuous monitoring points and self-adaptive coupling low-dimensional and high-dimensional calculation, the accuracy of coupling calculation of different dimensions can be greatly improved, the calculation efficiency is improved, and the errors of CAE simulation and test are reduced. By setting the intermittent monitoring points, the shock wave interruption and contact interface positions in the calculation domain can be accurately captured, and the method is a basis for self-adaptive conversion among models with different dimensions; the self-adaptive coupling can automatically reconstruct a low-dimensional model in high dimensions and map the model to a Euler grid in high dimensions, so that accurate reconstruction and mapping between grids in different dimensions and state data are realized, the prediction precision and calculation efficiency of simulation are improved, and the engineering application level is improved. Thus, the embodiment of the disclosure can effectively realize high-precision and high-efficiency simulation prediction on the complex engineering problem of far-field interaction.

Fig. 12 is a schematic structural diagram of a multi-dimensional simulation apparatus according to an embodiment of the present disclosure, which may be understood as the above-mentioned electronic device or a part of functional modules in the above-mentioned electronic device. As shown in fig. 12, the multi-dimensional simulation apparatus 1200 includes: the initialization module 1210 is configured to perform simulation initialization for the target engineering problem, where the simulation initialization includes establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem;

A setting and determining module 1220 for setting intermittent monitoring points for a target substance in the low-dimensional simulation model and determining the locations of the intermittent monitoring points, wherein the target substance interacts with a target intermittent, the target intermittent comprising a contact intermittent and/or a shock wave intermittent;

a first determining module 1230, configured to determine, for the low-dimensional simulation model, a current position of the target discontinuity and determine a current physical state of a substance in the low-dimensional simulation model, until the current position of the target discontinuity and a position of a discontinuity monitoring point meet a preset condition;

a reconstruction module 1240, configured to perform high-dimensional reconstruction on the low-dimensional simulation model and map the low-dimensional simulation model into the high-dimensional simulation model, so as to obtain a new high-dimensional simulation model;

a second determining module 1250 is configured to determine, for the new high-dimensional simulation model, a current position of the target discontinuity and a current physical state of the substance in the new high-dimensional simulation model until the simulation duration reaches a preset termination duration.

In another embodiment of the present disclosure, the preset condition is that a distance between the target discontinuity and the discontinuity monitoring point is less than a preset euler mesh cell size, wherein the distance between the target discontinuity and the discontinuity monitoring point is determined by a current position of the target discontinuity and a position of the discontinuity monitoring point.

In yet another embodiment of the present disclosure, the setup and determination module 1220 includes a setup submodule for setting intermittent monitoring points for a target substance in a low-dimensional simulation model, which may include:

the setting unit is used for setting intermittent monitoring points at the positions of a preset number of Euler grid units spaced from the target substance in the low-dimensional simulation model.

In yet another embodiment of the present disclosure, the number of intermittent monitoring points is a plurality, and the plurality of intermittent monitoring points surrounds the circumference of the target substance.

In yet another embodiment of the present disclosure, the initialization module 1210 includes a building sub-module for building a high-dimensional simulation model and a low-dimensional simulation model corresponding to a target engineering problem, where the building sub-module may include:

the building unit is used for building a corresponding high-dimensional substance model aiming at a substance in a target engineering problem so as to obtain a high-dimensional simulation model, wherein the high-dimensional substance model corresponding to a substance with a phase state of fluid is a Gao Weiou Law model, and the high-dimensional substance model corresponding to a substance with a phase state of solid is a Gao Weiou Law model or a high-dimensional Lagrange model;

and the dimension reduction unit is used for carrying out dimension reduction modeling on the high-dimension material model with the spherical symmetry and/or the axial symmetry to obtain a corresponding low-dimension material model, thereby obtaining a low-dimension simulation model.

In yet another embodiment of the present disclosure, the target discontinuity is a contact discontinuity; wherein the first determining module 1230 may include:

the first determining submodule is used for obtaining a corresponding current physical state by solving a preset low-dimensional control equation aiming at each substance and target interruption in the low-dimensional simulation model after decoupling interaction of different substances in the low-dimensional simulation model;

the second determining submodule is used for determining the current position of the target break according to the current physical state of the low-dimensional simulation model;

wherein the second determining module 1250 may include:

the third determining submodule is used for decoupling interaction of different substances in the new high-dimensional simulation model, and solving a corresponding current physical state by adopting a preset high-dimensional control equation according to each substance and target interruption in the new high-dimensional simulation model, wherein the current physical state is a physical state when a time propulsion step length is passed after the physical state is obtained by solving the current physical state relatively to the last time;

and the fourth determination submodule is used for determining the current position of the target break according to the current physical state of the new high-dimensional simulation model, wherein the current position is the position when the step length is advanced by time after the position is obtained by solving the current position for the last time.

In yet another embodiment of the present disclosure, the target discontinuity is a shockwave discontinuity, wherein the first determining module 1230 may include:

a fifth determination sub-module for determining a current location of the target discontinuity using the shock wave discontinuity indicator.

In yet another embodiment of the present disclosure, the reconstruction module 1240 may include:

the reconstruction sub-module is used for carrying out high-dimensional reconstruction on the low-dimensional simulation model to obtain a Gao Weichong structural model;

and the mapping sub-module is used for mapping the high-dimensional reconstruction model to the Euler grid in the high-dimensional simulation model to obtain a new high-dimensional simulation model.

The device provided in this embodiment can execute the method of any one of the above embodiments, and the execution mode and the beneficial effects thereof are similar, and are not described herein again.

The embodiment of the disclosure also provides an electronic device, which comprises: a memory in which a computer program is stored; a processor for executing the computer program, which when executed by the processor can implement the method of any of the above embodiments.

By way of example, fig. 13 is a schematic structural diagram of an electronic device in an embodiment of the present disclosure. Referring now in particular to fig. 13, a schematic diagram of an electronic device 1300 suitable for use in implementing embodiments of the present disclosure is shown. The electronic device 1300 in the embodiments of the present disclosure may include, but is not limited to, mobile terminals such as mobile phones, notebook computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), in-vehicle terminals (e.g., in-vehicle navigation terminals), and the like, and stationary terminals such as digital TVs, desktop computers, and the like. The electronic device shown in fig. 13 is merely an example and should not impose any limitations on the functionality and scope of use of embodiments of the present disclosure.

As shown in fig. 13, the electronic device 1300 may include a processing means (e.g., a central processor, a graphics processor, etc.) 1301, which may perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 1302 or a program loaded from a storage means 1308 into a Random Access Memory (RAM) 1303. In the RAM 1303, various programs and data necessary for the operation of the electronic apparatus 1300 are also stored. The processing device 1301, the ROM 1302, and the RAM 1303 are connected to each other through a bus 1304. An input/output (I/O) interface 1305 is also connected to bus 1304.

In general, the following devices may be connected to the I/O interface 1305: input devices 1306 including, for example, a touch screen, touchpad, keyboard, mouse, camera, microphone, accelerometer, gyroscope, and the like; an output device 1307 including, for example, a Liquid Crystal Display (LCD), a speaker, a vibrator, and the like; storage 1308 including, for example, magnetic tape, hard disk, etc.; and communication means 1309. The communication means 1309 may allow the electronic device 1300 to communicate with other devices wirelessly or by wire to exchange data. While fig. 13 shows an electronic device 1300 having various means, it is to be understood that not all illustrated means are required to be implemented or provided. More or fewer devices may be implemented or provided instead.

In particular, according to embodiments of the present disclosure, the processes described above with reference to flowcharts may be implemented as computer software programs. For example, embodiments of the present disclosure include a computer program product comprising a computer program embodied on a non-transitory computer readable medium, the computer program comprising program code for performing the method shown in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network via the communications device 1309, or installed from the storage device 1308, or installed from the ROM 1302. The above-described functions defined in the methods of the embodiments of the present disclosure are performed when the computer program is executed by the processing device 1301.

It should be noted that the computer readable medium described in the present disclosure may be a computer readable signal medium or a computer readable storage medium, or any combination of the two. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples of the computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this disclosure, a computer-readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present disclosure, however, the computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, with the computer-readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: electrical wires, fiber optic cables, RF (radio frequency), and the like, or any suitable combination of the foregoing.

In some implementations, the clients, servers may communicate using any currently known or future developed network protocol, such as HTTP (HyperText Transfer Protocol ), and may be interconnected with any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network ("LAN"), a wide area network ("WAN"), the internet (e.g., the internet), and peer-to-peer networks (e.g., ad hoc peer-to-peer networks), as well as any currently known or future developed networks.

The computer readable medium may be contained in the electronic device; or may exist alone without being incorporated into the electronic device.

The computer readable medium carries one or more programs which, when executed by the electronic device, cause the electronic device to: performing simulation initialization for the target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem; setting intermittent monitoring points for a target substance in a low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substance interacts with a target intermittent, and the target intermittent comprises a contact intermittent and/or a shock wave intermittent; aiming at the low-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of a substance in the low-dimensional simulation model until the current position of the target discontinuity and the position of a discontinuity monitoring point meet preset conditions; high-dimensional reconstruction is carried out on the low-dimensional simulation model and the high-dimensional reconstruction is mapped into the high-dimensional simulation model, so that a new high-dimensional simulation model is obtained; aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

Computer program code for carrying out operations of the present disclosure may be written in one or more programming languages, including, but not limited to, an object oriented programming language such as Java, smalltalk, C ++ and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (for example, through the Internet using an Internet service provider).

The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.

The units involved in the embodiments of the present disclosure may be implemented by means of software, or may be implemented by means of hardware. Wherein the names of the units do not constitute a limitation of the units themselves in some cases.

The functions described above herein may be performed, at least in part, by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: a Field Programmable Gate Array (FPGA), an Application Specific Integrated Circuit (ASIC), an Application Specific Standard Product (ASSP), a system on a chip (SOC), a Complex Programmable Logic Device (CPLD), and the like.

In the context of this disclosure, a machine-readable medium may be a tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. The machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. The machine-readable medium may include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples of a machine-readable storage medium would include an electrical connection based on one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

The embodiments of the present disclosure further provide a computer readable storage medium, where a computer program is stored, where the computer program, when executed by a processor, may implement a method according to any one of the foregoing embodiments, and the implementation manner and beneficial effects of the method are similar, and are not described herein again.

It should be noted that in this document, relational terms such as "first" and "second" and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.

The foregoing is merely a specific embodiment of the disclosure to enable one skilled in the art to understand or practice the disclosure. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the disclosure. Thus, the present disclosure is not intended to be limited to the embodiments shown and described herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A multi-dimensional simulation method, comprising:

performing simulation initialization for a target engineering problem, wherein the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem; the establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem comprises the following steps: establishing a corresponding high-dimensional material model aiming at the material in the target engineering problem, thereby obtaining the high-dimensional simulation model; performing dimension reduction modeling on the high-dimensional material model with the spherical symmetry and/or the axial symmetry to obtain a corresponding low-dimensional material model, thereby obtaining the low-dimensional simulation model;

Setting intermittent monitoring points for target substances in the low-dimensional simulation model and determining the positions of the intermittent monitoring points, wherein the target substances interact with target discontinuities, and the target discontinuities comprise contact discontinuities and/or shock wave discontinuities;

determining the current position of the target discontinuity and the current physical state of a substance in the low-dimensional simulation model aiming at the low-dimensional simulation model until the current position of the target discontinuity and the position of the discontinuity monitoring point meet preset conditions; the preset condition is that the distance between the target discontinuity and the discontinuity monitoring point is smaller than a preset distance;

performing high-dimensional reconstruction on the low-dimensional simulation model and mapping the high-dimensional reconstruction to the high-dimensional simulation model to obtain a new high-dimensional simulation model;

and aiming at the new high-dimensional simulation model, determining the current position of the target discontinuity and determining the current physical state of the substances in the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

2. The method of claim 1, wherein the preset condition is that a distance between the target discontinuity and the discontinuity monitoring point is less than a preset euler mesh cell size, wherein the distance between the target discontinuity and the discontinuity monitoring point is determined by a current location of the target discontinuity and a location of the discontinuity monitoring point.

3. The method of claim 1, wherein the setting intermittent monitoring points in the low-dimensional simulation model comprises:

in the low-dimensional simulation model, the intermittent monitoring points are set at positions spaced apart from the target substance by a preset number of Euler grid cells.

4. A method according to claim 3, wherein the number of intermittent monitoring points is a plurality, and a plurality of intermittent monitoring points are circumferentially around the target substance.

5. The method of claim 1, wherein the step of determining the position of the substrate comprises,

the high-dimensional material model corresponding to the material with the fluid phase is a Gao Weiou Law model, and the high-dimensional material model corresponding to the material with the solid phase is a Gao Weiou Law model or a high-dimensional Lagrange model.

6. The method of claim 1, wherein the target discontinuity is a contact discontinuity;

wherein the determining, for the low-dimensional simulation model, the current location of the target discontinuity and determining the current physical state of the substance in the low-dimensional simulation model includes:

after interaction of different substances in the low-dimensional simulation model is decoupled, aiming at each substance in the low-dimensional simulation model and the target discontinuity, a preset low-dimensional control equation is adopted to solve to obtain a corresponding current physical state;

Determining the current position of the target discontinuity according to the current physical state of the low-dimensional simulation model;

wherein determining, for the new high-dimensional simulation model, a current location of the target discontinuity and determining a current physical state of a substance in the new high-dimensional simulation model comprises:

after decoupling interaction of different substances in the new high-dimensional simulation model, solving a corresponding current physical state by adopting a preset high-dimensional control equation according to each substance in the new high-dimensional simulation model and the target discontinuity, wherein the current physical state is a physical state when a time propulsion step length is passed after the physical state is obtained by solving the current physical state relatively to the last time;

and determining the current position of the target break according to the current physical state of the new high-dimensional simulation model, wherein the current position is the position when the time propulsion step length passes after the position is obtained by solving the current position for the last time.

7. The method of claim 1, wherein the target discontinuity is a shock wave discontinuity;

wherein said determining the current location of the target discontinuity comprises:

a current location of the target discontinuity is determined using a shockwave discontinuity indicator.

8. The method of claim 1, wherein performing the high-dimensional reconstruction of the low-dimensional simulation model to obtain a high-dimensional simulation model comprises:

performing high-dimensional reconstruction on the low-dimensional simulation model to obtain a Gao Weichong structural model;

mapping the Gao Weichong structural model to Euler grids in the high-dimensional simulation model to obtain the new high-dimensional simulation model.

9. A multi-dimensional simulation apparatus, comprising:

the system comprises an initialization module, a simulation module and a control module, wherein the initialization module is used for performing simulation initialization aiming at a target engineering problem, and the simulation initialization comprises the steps of establishing a high-dimensional simulation model and a low-dimensional simulation model corresponding to the target engineering problem; the initialization module comprises: the building unit is used for building a corresponding high-dimensional substance model aiming at the substances in the target engineering problem so as to obtain the high-dimensional simulation model; the dimension reduction unit is used for carrying out dimension reduction modeling on the high-dimension material model with the spherical symmetry and/or the axial symmetry to obtain a corresponding low-dimension material model, so as to obtain the low-dimension simulation model;

the setting and determining module is used for setting intermittent monitoring points for target substances in the low-dimensional simulation model and determining positions of the intermittent monitoring points, wherein the target substances interact with target discontinuities, and the target discontinuities comprise contact discontinuities and/or shock wave discontinuities;

The first determining module is used for determining the current position of the target discontinuity and determining the current physical state of a substance in the low-dimensional simulation model aiming at the low-dimensional simulation model until the current position of the target discontinuity and the position of the discontinuity monitoring point meet preset conditions; the preset condition is that the distance between the target discontinuity and the discontinuity monitoring point is smaller than a preset distance;

and the second determining module is used for determining the current position of the target break and the current physical state of the substances in the new high-dimensional simulation model aiming at the new high-dimensional simulation model until the simulation duration reaches the preset termination duration.

10. An electronic device, comprising:

a processor and a memory, wherein the memory has stored therein a computer program which, when executed by the processor, performs the method of any of claims 1-8.

11. A computer readable storage medium, characterized in that the storage medium has stored therein a computer program which, when executed by a processor, implements the method according to any of claims 1-8.