CN113158520B - Fuel ice layer interface tracking simulation method for freezing target system - Google Patents
Fuel ice layer interface tracking simulation method for freezing target system Download PDFInfo
- Publication number
- CN113158520B CN113158520B CN202110384783.8A CN202110384783A CN113158520B CN 113158520 B CN113158520 B CN 113158520B CN 202110384783 A CN202110384783 A CN 202110384783A CN 113158520 B CN113158520 B CN 113158520B
- Authority
- CN
- China
- Prior art keywords
- fuel
- ice layer
- moment
- temperature
- fuel ice
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/10—Nuclear fusion reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A fuel ice layer interface tracking simulation method for a freezing target system is characterized in that an Euler coordinate system and a Lagrange coordinate system are combined, data transfer between the two coordinate systems is completed through Gaussian kernel function estimation, a heat conduction differential equation is reconstructed by combining a Stefan equation, tracking calculation of a fuel ice layer interface is completed, profile distribution of a solid fuel ice layer in fuel target pellets in the freezing target system under different boundary loads is obtained, and requirements of low mode roughness of ice layer uniformity ignition being less than 0.5 mu m, ice layer thickness being +/-3 mu m and gas density are met; the method can more intuitively and comprehensively represent the influence of the temperature field on the morphology of the fuel ice layer in the target pellet and provide theoretical technical guidance for experiments.
Description
Technical Field
The invention belongs to the technical field of design simulation of a freezing target system, and particularly relates to a fuel ice layer interface tracking simulation method for the freezing target system.
Background
The inertial confinement nuclear fusion is to irradiate the surface of the frozen target pellet with laser uniformly so as to achieve the ignition condition of high temperature and high density and realize the fusion reaction. In order to meet the ignition requirement and avoid Rayleigh-Taylor instability, the thickness uniformity of a deuterium-deuterium fuel ice layer in a target pellet needs to be more than 99 percent, and the root-mean-square roughness of the inner surface of the fuel ice layer needs to be less than 1 mu m. The low mode roughness of the fuel ice layer is mainly determined by the temperature field around the target pellet, and therefore the importance of the control of the temperature field of the frozen target is particularly prominent.
At present, the uniformity of the fuel ice layer is mainly characterized by an experimental mode, and specific characterization methods comprise LED backlight characterization and X-ray phase contrast imaging characterization. Since the fuel ice layer profile within the target pellet can only be observed in a specific direction (north-south dipolar direction, equatorial direction), a large amount of information on the fuel ice layer profile is masked. In order to more intuitively and comprehensively characterize the influence of a temperature field on the morphology of a fuel ice layer in a target pellet and provide theoretical technical guidance for experiments, a fuel ice layer interface tracking simulation method for a freezing target system needs to be established urgently.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention aims to provide a fuel ice layer interface tracking simulation method for a freezing target system, which can more intuitively and comprehensively represent the influence of a temperature field on the appearance of a fuel ice layer in a target pellet and provide theoretical technical guidance for experiments.
In order to achieve the purpose, the invention adopts the technical scheme that:
a tracking simulation method for a fuel ice layer interface in a freezing target system is characterized in that an Euler coordinate system and a Lagrange coordinate system are combined, data transfer between the two coordinate systems is completed through Gaussian kernel function estimation, a heat conduction differential equation is reconstructed by combining a Stefan equation, tracking calculation of the fuel ice layer interface is completed, profile distribution of a solid fuel ice layer in fuel target pellets in the freezing target system under different boundary loads is obtained, and requirements of low mode roughness of ice layer uniformity ignition being less than 0.5 mu m, ice layer thickness being +/-3 mu m and gas density are met.
A method for simulating fuel ice layer interface tracking in a freeze target system, comprising the steps of:
(1) Setting the initial morphology of the fuel ice layer, the temperature distribution of the outer surface of the target pellet and the density parameter of the fuel gas;
(2) Calculating t m The actual volume of the fuel ice layer in the grid of the Euler coordinate system at the moment is beta:
assuming that the arrangement of the fuel ice layer contour particles within a single Euler grid is approximately a straight line, β is expressed as:
wherein, deltax, delay and Delaz are the grid sizes in the directions of x, y and z respectively, and the central point of the target pill is arranged asThe position coordinate of the central point of the Euler grid isThe mean position vector of all fuel ice layer profile particles in a single Euler grid;
(3) Obtaining the corresponding saturation state temperature according to the density of the fuel gas, and loading the saturation state temperature to the boundary of the inner surface of the fuel ice layer;
(4) Calculating t m Temperature field distribution under the time Euler coordinate system:
considering the interface temperature of the fuel ice layer to be uniform at each moment, wherein the value of the interface temperature is equal to the phase change saturation temperature corresponding to the density of the fuel gas at the moment; meanwhile, equivalent treatment is carried out on the fuel gas as a solid, under an Euler coordinate, only a temperature field of a solid domain needs to be solved, and the adopted control equation is an unsteady state heat conduction equation:
where T is the temperature, T is the time, α eff For the equivalent thermal diffusivity, the expression is as follows:
where k is the thermal conductivity, ρ is the density, c p Constant pressure heat capacity;
(5) According to the temperature field distribution and the Stefan equation, calculating t under the Euler coordinate system m The normal displacement speed of the fuel ice layer interface at the moment:
the solution of the normal moving speed of the interface of the fuel ice layer is based on a Stefan equation, and the form of the Stefan equation is as follows:
where v is the normal moving velocity, Δ H is the latent heat, subscript ice denotes the fuel solid, gas denotes the fuel gas;
(6) Calculating t under a Lagrange coordinate system through Gaussian kernel function estimation m Normal displacement speed of profile particles of different fuel ice layers at different moments:
according to the distribution condition of each fuel ice layer profile particle in an Euler grid, determining the velocity of each particle by a kernel function estimation method, wherein the kernel function is a Gaussian function and has the form:
wherein b is taken as the average size of the grid, i.e. b = (Δ x + Δ y + Δ z)/3, Δ x, Δ y and Δ z are the grid sizes in the x, y and z directions, respectively;
if t m The position coordinate of the contour particle of the s th fuel ice layer at the moment isThe position coordinate of the central point of the Euler grid isThe normal moving speed of the profile of the fuel ice layer obtained by calculation of the Euler grid is v i,j,k That t m At the time of the normal direction of the profile pointThe moving speed is as follows:
(7) Calculating t m+1 Position vector of the fuel ice layer profile particle at time:
fuel ice layer profile particle at t m After the determination of the normal displacement speed at the moment, it is at t m+1 The position vector of the time is:
(8) According to t m+1 Calculating t from the position vector of the profile particle of the fuel ice layer at the moment m+1 The mass of the fuel ice layer at the moment of time, thereby obtaining t m+1 The time-of-day fuel gas density; meanwhile, calculating the actual volume ratio beta of the fuel ice layer in the grid of the Euler coordinate system through Gaussian kernel function estimation;
the density of the fuel gas is indirectly obtained by the volume change of the fuel ice layer, the total number of the profile particles of the fuel ice layer is recorded as N, the radius of the inner surface of the target pellet is recorded as R, and the central point of the target pellet is positioned asThen t m The mass of the fuel ice layer at the moment is:
the mass of the fuel ice layer and the fuel gas at the initial time is recorded as m ice (0) And m gas (0) Then t is m The fuel gas density at that time is:
(9) Returning to the step (2), calculating the next moment until the morphology of the fuel ice layer is basically unchanged or the next moment is calculated to the specified moment, and ending the calculation; the profile distribution of the solid fuel ice layer in the fuel target pellets in the freezing target system under different boundary loads is obtained, and the requirements of low mode roughness of ice layer uniformity ignition being less than 0.5 mu m, ice layer thickness being +/-3 mu m and gas density are met.
In the step (1), the temperature distribution of the outer surface of the target pill is in a distribution trend of high temperature in the north and south and low temperature in the equator.
And (2) the initial appearance of the fuel ice layer in the step (1) is characterized by the profile particles of the fuel ice layer in a Lagrange coordinate system.
In the step (3), because the target pellet has a small size, the time constant of thermal diffusion is compared with the calculated time step length to be a small amount, and the temperature of the fuel gas in the target pellet is considered to be equal everywhere.
In the step (3), the fuel gas in the target pellet is considered to be in a phase-change saturation state all the time.
In the step (8), the temperature of the fuel gas in the saturation state is found by the NIST REFPROP database under the premise that the density of the fuel gas in the saturation state is known.
The invention has the beneficial effects that:
the invention discloses a fuel ice layer interface tracking simulation method for a freezing target system, which can calculate the outline migration condition of an ice layer in a target pellet according to different thermal boundary loads and solve the problem that the appearance of the ice layer cannot be completely observed and the uniformity of the ice layer cannot be evaluated in an experiment. By the method, the uniformity of the ice layer in the target pellet under different thermal boundary loads can be measured, a set of control flow of the freezing target temperature field can be reversely provided, and clear theoretical guidance is provided for experiments.
Drawings
Fig. 1 is a schematic view of a target pellet structure according to an embodiment of the present invention.
FIG. 2 is a flow chart of the method of the present invention.
Fig. 3 is a calculation result of the embodiment of the present invention.
Detailed Description
The technical solutions of the present invention are described below clearly and completely with reference to the following embodiments, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all 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.
The target pellet of this embodiment has a structure as shown in fig. 1, the outermost CH target shell 1, an outer diameter of 1.1mm and an inner diameter of 0.91mm is attached to the inner side of the CH target shell 1, a solid fuel ice layer 2 with a thickness of 68 μm is attached to the inner side of the CH target shell 1, and a fuel gas cavity 3 is provided in the center of the target pellet.
Referring to fig. 2, a method for fuel ice layer interface tracking simulation in a freezing target system includes the steps of:
(1) Setting the initial morphology of the fuel ice layer, the temperature distribution of the outer surface of the target pellet and the density parameter of the fuel gas;
setting the initial morphology of the fuel ice layer to be uniform distribution, setting the temperature distribution of the outer surface of the target pellet to be 18.05K at the north and south poles and 18.00K at the equator, and performing linear transition, wherein the temperature distribution of the outer surface of the target pellet is the distribution trend that the north and south poles are high and the equator is low; the initial appearance of the fuel ice layer is characterized by the profile particles of the fuel ice layer of a Lagrange coordinate system;
(2) Calculating t m The actual volume ratio beta of the fuel ice layer in the grid of the Euler coordinate system at the moment;
at known t m On the premise of position vector of contour particles of fuel ice layer at moment, t is obtained m The actual volume ratio beta of the fuel ice layer in the Euler grid of the fuel ice layer interface exists at any moment; assuming that the arrangement of the fuel ice layer contour particles within a single Euler grid is approximately a straight line, β is expressed as:
wherein, deltax, delay and Delaz are the grid sizes in the directions of x, y and z respectively, and the central point of the target pill is arranged asEuler gridThe position coordinate of the central point isThe mean position vector of all fuel ice layer profile particles in a single Euler grid; for a non-interface euler grid, β =1 indicates that all fuel ice layers are in the euler grid, and β =0 indicates that all fuel gas is in the euler grid;
(3) Obtaining the corresponding saturation state temperature according to the density of the fuel gas, and loading the saturation state temperature to the boundary of the inner surface of the fuel ice layer; considering that the fuel gas in the target pellet is always in a phase-change saturation state;
(4) Calculating t m Distributing the temperature field under the Euler coordinate system at the moment;
because the fuel gas has no pressure and temperature gradient, the interface temperature of the fuel ice layer can be considered to be uniform at each moment, and the value of the interface temperature is equal to the phase change saturation temperature corresponding to the density of the fuel gas at the moment; meanwhile, the fuel gas is treated equivalently as a solid, so that the calculation requirement is reduced, therefore, under the Euler coordinate, only the temperature field of the solid domain needs to be solved, and the adopted control equation is an unsteady state heat conduction equation:
where T is the temperature, T is the time, α eff For the equivalent thermal diffusivity, the expression is as follows:
where k is the thermal conductivity, ρ is the density, c p Constant pressure heat capacity;
(5) According to the temperature field distribution and the Stefan equation, calculating t under the Euler coordinate system m The normal displacement speed of the interface of the fuel ice layer at any moment;
the solution of the normal moving speed of the interface of the fuel ice layer is based on a Stefan equation, and the form of the Stefan equation is as follows:
where v is the normal moving velocity, Δ H is the latent heat, subscript ice denotes the fuel solid, gas denotes the fuel gas;
(6) Calculating t under a Lagrange coordinate system through Gaussian kernel function estimation m The normal displacement speeds of profile particles of different fuel ice layers at different moments;
calculating t from Stefan equation m After the normal moving speed of the fuel ice layer interface on the Euler grid node is reached, the speed of each particle is determined by a kernel function estimation method according to the distribution condition of each fuel ice layer profile particle in the Euler grid, wherein the kernel function is a Gaussian function with higher precision and has the form:
wherein b is taken as the average size of the grid, i.e. b = (Δ x + Δ y + Δ z)/3, Δ x, Δ y and Δ z are the grid sizes in the x, y and z directions, respectively;
if t m The position coordinate of the contour particle of the s th fuel ice layer at the moment isThe position coordinate of the central point of the Euler grid isThe normal moving speed of the profile of the fuel ice layer obtained by calculation of the Euler grid is v i,j,k That t m The normal moving speed of the profile particle at the moment is as follows:
(7) Calculating t m+1 The position vector of the profile particle of the fuel ice layer at the moment;
fuel ice layer profilePoint is at t m After the normal movement speed at time is determined, it is at t m+1 The position vector of the time is:
(8) According to t m+1 Calculating t from the position vector of the profile particle of the fuel ice layer at the moment m+1 The mass of the fuel ice layer at the moment of time, thereby obtaining t m+1 The time-of-day fuel gas density; meanwhile, calculating the actual volume ratio beta of the fuel ice layer in the grid of the Euler coordinate system through Gaussian kernel function estimation;
in the phase change process, the change of the shape of the fuel ice layer inevitably causes the change of the density of the fuel gas at the center of the target pellet, and the density of the fuel gas simultaneously determines the temperature of the inner surface of the fuel ice layer; the density of fuel gas is indirectly obtained by the volume change of the fuel ice layer, and on the premise that the thickness of the local fuel ice layer is far smaller than the inner surface radius of the target pellet, the mass of the fuel ice layer is obtained by approximating the position coordinates of the profile particles of the fuel ice layer, wherein the total number of the profile particles of the fuel ice layer is N, the inner surface radius of the target pellet is R, and the central point of the target pellet is positioned asThen t m The mass of the fuel ice layer at the moment is:
the mass of the fuel ice layer and the fuel gas at the initial time is recorded as m ice (0) And m gas (0) Then t is m The fuel gas density at that time is:
on the premise of knowing the density of the fuel gas in the saturation state, the temperature of the fuel gas in the saturation state is found through a NIST REFPROP database;
(9) Returning to the step (2), calculating the next moment until 1000s, and finishing the calculation; the profile distribution of the solid fuel ice layer in the fuel target pellets in the freezing target system under different boundary loads is obtained, and the requirements of the ice layer uniformity ignition low-mode roughness less than 0.5 mu m, the ice layer thickness +/-3 mu m and the gas density are met;
the thickness distribution curves of the fuel ice layers at different moments are shown in FIG. 3, wherein 90 degrees represents the south pole, 90 degrees represents the north pole, and 0 degrees represents the equator; because the temperature of the outer surface of the target pellets is in the distribution trend of high temperature in the south and north and low temperature in the equator, the fuel ice layer in the south and north poles is gradually thinned, and the fuel ice layer in the equator is gradually thickened.
Claims (6)
1. A method for simulating fuel ice layer interface tracking in a freezing target system, comprising the steps of:
(1) Setting the initial morphology of the fuel ice layer, the temperature distribution of the outer surface of the target pellet and the density parameter of the fuel gas;
(2) Calculating t m The actual volume of the fuel ice layer in the grid of the Euler coordinate system at the moment accounts for beta:
assuming that the arrangement of the fuel ice layer contour particles within a single Euler grid is approximately a straight line, β is expressed as:
wherein, deltax, delay and Delaz are the grid sizes in the directions of x, y and z respectively, and the central point of the target pill is arranged asThe position coordinate of the central point of the Euler grid isThe mean position vector of all fuel ice layer profile particles in a single Euler grid;
(3) Obtaining the corresponding saturation state temperature according to the density of the fuel gas, and loading the saturation state temperature to the boundary of the inner surface of the fuel ice layer;
(4) Calculating t m Temperature field distribution under the Euler coordinate system at the moment:
considering the interface temperature of the fuel ice layer to be uniform at each moment, wherein the value of the interface temperature is equal to the phase change saturation temperature corresponding to the density of the fuel gas at the moment; meanwhile, equivalent treatment is carried out on the fuel gas as a solid, under an Euler coordinate, only a temperature field of a solid domain needs to be solved, and the adopted control equation is an unsteady state heat conduction equation:
where T is the temperature, T is the time, α eff For the equivalent thermal diffusivity, the expression is as follows:
where k is the thermal conductivity, ρ is the density, c p Constant pressure heat capacity;
(5) According to the temperature field distribution and the Stefan equation, calculating t under the Euler coordinate system m The normal displacement speed of the interface of the fuel ice layer at the moment:
the solution of the normal moving speed of the interface of the fuel ice layer is based on a Stefan equation, and the form of the Stefan equation is as follows:
where v is the normal moving velocity, Δ H is the latent heat, subscript ice denotes the fuel solid, gas denotes the fuel gas;
(6) Calculating t under a Lagrange coordinate system through Gaussian kernel function estimation m Normal displacement velocity of profile particles of different fuel ice layers at different times:
according to the distribution condition of each fuel ice layer profile particle in an Euler grid, determining the velocity of each particle by a kernel function estimation method, wherein the kernel function is a Gaussian function and has the form:
wherein b is the average size of the grid, i.e. b = (Δ x + Δ y + Δ z)/3, and Δ x, Δ y and Δ z are the grid sizes in the x, y and z directions, respectively;
if t m The position coordinate of the contour particle of the s th fuel ice layer at the moment isThe position coordinate of the central point of the Euler grid isThe normal moving speed of the profile of the fuel ice layer obtained by calculation of the Euler grid is v i,j,k That t m The normal moving speed of the profile particle at the moment is as follows:
(7) Calculating t m+1 Position vector of the fuel ice layer profile particle at time:
fuel ice layer profile particle at t m After the normal movement speed at time is determined, it is at t m+1 The position vector of the time is:
(8) According to t m+1 Calculating t from the position vector of the profile particle of the fuel ice layer at the moment m+1 The mass of the fuel ice layer at the moment of time, thereby obtaining t m+1 The time-of-day fuel gas density; meanwhile, calculating the actual volume ratio beta of the fuel ice layer in the grid of the Euler coordinate system through Gaussian kernel function estimation;
the density of fuel gas is indirectly obtained by the volume change of the fuel ice layer, the total number of contour particles of the fuel ice layer is recorded as N, the radius of the inner surface of the target pellet is recorded as R, and the central point of the target pellet is positioned asThen t m The mass of the fuel ice layer at the moment is:
the mass of the fuel ice layer and the fuel gas at the initial time is recorded as m ice (0) And m gas (0) Then t is m The fuel gas density at that time is:
(9) Returning to the step (2), calculating the next moment until the morphology of the fuel ice layer is basically unchanged or the next moment is calculated to the specified moment, and ending the calculation; the profile distribution of the solid fuel ice layer in the fuel target pellets in the freezing target system under different boundary loads is obtained, and the requirements of low mode roughness of ice layer uniformity ignition less than 0.5 mu m, ice layer thickness +/-3 mu m and gas density are met.
2. The method for simulating the tracking of the fuel ice layer interface in the freezing target system according to claim 1, wherein: in the step (1), the temperature distribution of the outer surface of the target pill is in a distribution trend of high temperature in the north and south and low temperature in the equator.
3. The method for simulating the tracking of the fuel ice layer interface in the freezing target system according to claim 1, wherein: and (2) in the step (1), the initial appearance of the fuel ice layer is characterized by the profile particles of the fuel ice layer in a Lagrange coordinate system.
4. The method for simulating the tracking of the fuel ice layer interface in the freezing target system according to claim 1, wherein: in the step (3), because the target pellet has a small size, the time constant of thermal diffusion is compared with the calculated time step length to be a small amount, and the temperature of the fuel gas in the target pellet is considered to be equal everywhere.
5. The method for simulating the tracking of the fuel ice layer interface in the freezing target system according to claim 1, wherein: in the step (3), the fuel gas in the target pellet is considered to be in a phase-change saturation state all the time.
6. The method for simulating the tracking of the fuel ice layer interface in the freezing target system according to claim 1, wherein: in the step (8), the temperature of the fuel gas in the saturation state is found by the NIST REFPROP database under the premise that the density of the fuel gas in the saturation state is known.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110384783.8A CN113158520B (en) | 2021-04-09 | 2021-04-09 | Fuel ice layer interface tracking simulation method for freezing target system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110384783.8A CN113158520B (en) | 2021-04-09 | 2021-04-09 | Fuel ice layer interface tracking simulation method for freezing target system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113158520A CN113158520A (en) | 2021-07-23 |
CN113158520B true CN113158520B (en) | 2022-10-28 |
Family
ID=76889737
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110384783.8A Active CN113158520B (en) | 2021-04-09 | 2021-04-09 | Fuel ice layer interface tracking simulation method for freezing target system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113158520B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113506440B (en) * | 2021-09-08 | 2021-11-30 | 四川国蓝中天环境科技集团有限公司 | Traffic state estimation method for multi-source data fusion under Lagrange coordinate system |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093835A (en) * | 2013-01-25 | 2013-05-08 | 西安交通大学 | Fast ignition target fuel layering device based on fluid magnetic suspension and fast ignition target fuel layering method based on fluid magnetic suspension |
WO2013078629A1 (en) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | Numerical simulation method for aircraft flight icing |
CN104268943A (en) * | 2014-09-28 | 2015-01-07 | 北京航空航天大学 | Fluid simulation method based on Eulerian-Lagrangian coupling method |
CN106297902A (en) * | 2016-08-31 | 2017-01-04 | 西安交通大学 | A kind of filling aeroge ICF cryogenic target system |
CN108877960A (en) * | 2017-05-15 | 2018-11-23 | 陈鹏玮 | ICF freezes target assembly and ICF shielding case opening speed best practice |
WO2020056177A1 (en) * | 2018-09-12 | 2020-03-19 | The Regents Of The University Of Colorado, A Body Corporate | Cryogenically cooled vacuum chamber radiation shields for ultra-low temperature experiments and extreme high vacuum (xhv) conditions |
CN111241742A (en) * | 2019-12-27 | 2020-06-05 | 西安交通大学 | Multiphase flow calculation method |
RU2730003C1 (en) * | 2019-10-22 | 2020-08-14 | Федеральное государственное бюджетное учреждение "Арктический и Антарктический научно-исследовательский институт" (ФГБУ "ААНИИ") | Method of determining physical and mechanical and morphometric characteristics of ice hummock formations |
-
2021
- 2021-04-09 CN CN202110384783.8A patent/CN113158520B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013078629A1 (en) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | Numerical simulation method for aircraft flight icing |
CN103093835A (en) * | 2013-01-25 | 2013-05-08 | 西安交通大学 | Fast ignition target fuel layering device based on fluid magnetic suspension and fast ignition target fuel layering method based on fluid magnetic suspension |
CN104268943A (en) * | 2014-09-28 | 2015-01-07 | 北京航空航天大学 | Fluid simulation method based on Eulerian-Lagrangian coupling method |
CN106297902A (en) * | 2016-08-31 | 2017-01-04 | 西安交通大学 | A kind of filling aeroge ICF cryogenic target system |
CN108877960A (en) * | 2017-05-15 | 2018-11-23 | 陈鹏玮 | ICF freezes target assembly and ICF shielding case opening speed best practice |
WO2020056177A1 (en) * | 2018-09-12 | 2020-03-19 | The Regents Of The University Of Colorado, A Body Corporate | Cryogenically cooled vacuum chamber radiation shields for ultra-low temperature experiments and extreme high vacuum (xhv) conditions |
RU2730003C1 (en) * | 2019-10-22 | 2020-08-14 | Федеральное государственное бюджетное учреждение "Арктический и Антарктический научно-исследовательский институт" (ФГБУ "ААНИИ") | Method of determining physical and mechanical and morphometric characteristics of ice hummock formations |
CN111241742A (en) * | 2019-12-27 | 2020-06-05 | 西安交通大学 | Multiphase flow calculation method |
Non-Patent Citations (6)
Title |
---|
Numerical analysis of dynamic heating modulation during rapid cooling of fuel layer in an indirect-drive cryogenic target;Zhao, J等;《PROGRESS IN NUCLEAR ENERGY》;20190731;第22-30页 * |
Recent progress in cryogenic techniques applied to pulsed ion beam generation and target ablation experiments;K. Kasuya等;《Digest of Technical Papers. 12th IEEE International Pulsed Power Conference. (Cat. No.99CH36358)》;19991231;第2卷;第1014-1017页 * |
冷冻靶降温过程温度场分布及冰层生存时间研究;陈洵等;《原子能科学技术》;20201231(第07期);第183-190页 * |
基于温度冲击的冷冻靶氘氘冰层结晶生长行为研究;陶朝友等;《原子能科学技术》;20201231;第54卷(第11期);第2209-2216页 * |
基于温度梯度的ICF冷冻靶均化技术研究;王凯等;《强激光与粒子束》;20150331;第27卷(第03期);第032025:1-7页 * |
靶丸材料及氘氚冰层分布对冷冻靶温度场的影响;李翠等;《制冷学报》;20200430;第41卷(第02期);第34-40页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113158520A (en) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Eccentricity optimization of a horizontal shell-and-tube latent-heat thermal energy storage unit based on melting and melting-solidifying performance | |
CN108423154B (en) | Hypersonic aircraft leading edge thermal protection method based on gradient porous material | |
CN113158520B (en) | Fuel ice layer interface tracking simulation method for freezing target system | |
CN110955991B (en) | Fluid-solid coupling calculation method for interface bidirectional data exchange | |
Qian et al. | An experimental investigation on transpiration cooling performances using solid hydrogel as coolant | |
CN109376403A (en) | A kind of two-dimentional icing method for numerical simulation based on Descartes's self-adapting reconstruction technology | |
CN107833275A (en) | A kind of dynamic ice microstructure three-dimensional modeling method | |
CN106709137A (en) | Polynomial maximum likelihood estimation-based soft measurement method for temperature field of reactor core of nuclear power plant | |
He et al. | LBM prediction of effective thermal conductivity of lithium-ion battery graphite anode | |
Wang et al. | Penetrative turbulent Rayleigh–Bénard convection in two and three dimensions | |
US20240135057A1 (en) | Method for numerical reconstruction and heat transfer characteristics evaluation of microstructure of thermal barrier coatings containing microcracks | |
CN111753250B (en) | One-dimensional unsteady state heat conduction inverse problem method | |
CN112613243A (en) | Method and device for fluid mechanics simulation and computer readable storage medium | |
Tian et al. | Sea urchin skeleton-inspired triply periodic foams for fast latent heat storage | |
CN117332511A (en) | Self-adaptive coupling numerical simulation method for high-temperature unbalanced flow of hypersonic vehicle | |
CN109192251A (en) | A kind of calculation method of the three-dimensional fine power distribution of solid fuel molten salt reactor | |
Guo et al. | Studying the advantages of equal curvature curved fin to enhance phase change heat storage | |
Wang et al. | Extended development of a bubble percolation method to predict boiling crisis of flow boiling | |
CN111036909B (en) | Grain structure numerical value prediction method in metal rapid forming process | |
CN116542116B (en) | Hypersonic aircraft full flow field numerical simulation method considering micropore gaseous sweating | |
CN108920872A (en) | For the BCP particle localization realization method and system of DSMC method | |
CN116213762A (en) | Microstructure morphology prediction method for laser metal additive manufacturing fused deposition layer | |
CN115542746A (en) | Energy control reentry guidance method and device for hypersonic aircraft | |
CN114330080A (en) | Prediction method of aircraft surface-symmetric cross-basin flow field | |
CN110758775B (en) | Multi-pulse area hovering method based on asteroid surface observation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |