CN112651083A - Detonation laser thruster analysis method based on FLUENT simulation - Google Patents

Detonation laser thruster analysis method based on FLUENT simulation Download PDF

Info

Publication number
CN112651083A
CN112651083A CN202011621232.0A CN202011621232A CN112651083A CN 112651083 A CN112651083 A CN 112651083A CN 202011621232 A CN202011621232 A CN 202011621232A CN 112651083 A CN112651083 A CN 112651083A
Authority
CN
China
Prior art keywords
fluent
grid
boundary
setting
calculation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011621232.0A
Other languages
Chinese (zh)
Other versions
CN112651083B (en
Inventor
檀朝彬
展先彪
王正一
李凌霄
孟令坤
辛尚青
郑国浩
朱磊
付智超
梁爽
朱亚琦
曹烨彤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Institute of Aerospace Engineering
Original Assignee
North China Institute of Aerospace Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Institute of Aerospace Engineering filed Critical North China Institute of Aerospace Engineering
Priority to CN202011621232.0A priority Critical patent/CN112651083B/en
Publication of CN112651083A publication Critical patent/CN112651083A/en
Application granted granted Critical
Publication of CN112651083B publication Critical patent/CN112651083B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a detonation laser thruster analysis method based on FLUENT simulation, which comprises the following steps of: s1, determining a control equation according to the working process of the knock mode; s2, constructing a user-defined function UDF; s3, establishing a parabolic model, modeling by using Gambit software, selecting a grid type to divide a grid, setting boundary conditions, and outputting a case file; s4, opening FLUENT software, and importing the case file into the FLUENT software; s5, detecting grids and setting grid units; s6, importing a custom function UDF by FLUENT software; s7, selecting a correct radiation model, an energy model and a turbulence model, and setting material properties and flow field pressure; s8, setting a moving grid; s9, selecting a solving method, and determining solving precision; s10, initializing a boundary, setting iteration parameters and calculating; s11, look at residual map, velocity, pressure and temperature map. The invention discusses the change reason of the thrust in the laser propelling process by using FLUENT simulation software, and can ensure the stable speed in the rocket flying process.

Description

Detonation laser thruster analysis method based on FLUENT simulation
Technical Field
The invention relates to the technical field of laser propulsion, in particular to a detonation laser thruster analysis method based on FLUENT simulation.
Background
The propulsion technology is the core technology of a space carrier, most of the existing space propulsion devices are liquid or solid rocket engines, and the existing space propulsion devices have the common point that heat energy generated by chemical reaction is converted into jet kinetic energy through a jet pipe, and thrust is generated by utilizing the recoil action. Research and practice have shown that for conventional chemical rocket engines, the specific impulse of the engine does not exceed 5000m/s due to the limitations of combustion chamber temperature and combustion product molecular weight. In order to improve the specific impulse, various non-chemical propulsion modes are researched, such as electric propulsion, nuclear energy propulsion, laser propulsion and the like, and an electric propulsion system can obtain high specific impulse but has low thrust; the nuclear propulsion can obtain great thrust and high specific impulse, but the propellant is toxic and has pollution, and large energy equipment needs to be carried; the laser propulsion system forms high-temperature and high-pressure plasma by focusing laser breakdown working media, the plasma is ejected out from the spray pipe to generate recoil force on the engine, the working media of the laser propulsion system are separated from energy, high-energy laser transmitted from a long distance provides energy required by propulsion, and high specific impulse and thrust can be obtained. Aiming at the wide application prospect of the laser propulsion system, an analysis platform aiming at the laser propulsion process needs to be provided urgently, and convenience is provided for the process research of laser propulsion.
Disclosure of Invention
The invention aims to provide a detonation laser thruster analysis method based on FLUENT simulation, which is characterized in that dynamic mesh modeling is applied, a user-defined function UDF is constructed, and the variation of thrust in the detonation laser propulsion process is researched by simulating by utilizing FLUENT software.
In order to achieve the purpose, the invention provides the following scheme:
a detonation laser thruster analysis method based on FLUENT simulation comprises the following steps:
s1, determining a control equation according to a laser wave absorption stage and a subsequent shock wave attenuation stage in the working process of the detonation mode;
s2, constructing a user-defined function UDF for laser energy setting, dynamic boundary setting and plasma parameter adjustment in the FLUENT calculation process;
s3, establishing a parabolic model, modeling by using Gambit software, selecting a grid type to divide a grid, setting boundary conditions, and outputting a case file;
s4, opening FLUENT software, and importing the case file into the FLUENT software;
s5, detecting the grids, checking whether the divided grids are correct and setting grid units;
s6, importing a custom function UDF by FLUENT software;
s7, selecting a correct radiation model, an energy model and a turbulence model, and setting material properties, flow field pressure and temperature;
s8, in FLUENT, performing dynamic grid calculation, selecting a dynamic grid calculation module, and setting the overall calculation parameters and modes of the dynamic grid;
s9, selecting a solving method, and determining solving precision;
s10, initializing a boundary, setting iteration parameters, and calculating by using a FLUENT solver;
s11, look at residual map, velocity, pressure and temperature map.
Optionally, the control equation in step S1 specifically includes:
Figure BDA0002876193800000021
Figure BDA0002876193800000022
Figure BDA0002876193800000023
Figure BDA0002876193800000024
the corresponding heat flow is also simplified to:
Figure BDA0002876193800000025
in the formula, U is a conservation variable vector, E, F is convection flux in x and r directions under a cylindrical coordinate system, Ev and Fv are respectively convection flux in x and r directions, G, Gv is a source item generated by adopting the cylindrical coordinate system for an unbonded part and a bonded part, W is a mass and energy source item caused by a physicochemical process, and q isx、qrTotal and vibrational heat flow, q, respectivelynCombining the total heat conduction and diffusion heat flows after each internal energy mode for the heat flow q, wherein rho is the total density of the mixed gas, p is the pressure of the mixed gas, and u and v are two coordinates respectivelyThe velocity component in the direction, H is the enthalpy per unit mass of the mixed gas, and τ is the stress tensor.
Optionally, a custom function UDF is constructed in step S2, and is used for laser energy setting, dynamic boundary setting, and plasma parameter adjustment in the FLUENT calculation process, which specifically includes:
injecting laser energy: the instantaneous energy injection model is adopted, plasma ionization is complete and no energy loss exists, laser pulse energy is completely absorbed by the plasma with certain deposition efficiency, the shape of the plasma is a sphere, and if the radius is r, the formula is obtained:
Figure BDA0002876193800000031
calculating a moving boundary: when the flexible boundary calculation is carried out, a dynamic grid calculation module needs to be selected, and the DEFINE _ CG _ MOTION macro in the UDF is used for writing and realizing functions;
implementation of plasma equation: the method comprises the steps of modifying physical quantities of a flow field and calculating integration by utilizing a general solution macro DEFINE _ ADJUST in FLUENT, wherein DEFINE _ ADJUST mainly integrates the whole flow field, then the boundary is adjusted by analyzing a calculation result, and the macro defined by the function needs to be called for each iteration of each step or each solution of a transport equation.
Optionally, the grid unit in step S5 is mm, the flow field pressure in step S7 is 1atm, and the temperature is 298.15K.
Optionally, in step S8, the dynamic mesh calculation is performed in FLUENT, a dynamic mesh calculation module is selected, and the dynamic mesh overall calculation parameters and modes are set, which specifically includes:
processing the propeller as a dynamic boundary by using a dynamic mesh module provided by FLUENT, calling UDF to read flow field information when the calculation of each time step is finished in the calculation, and calculating the speed obtained by the propeller according to the force borne by the wall surface;
when the moving grid is selected for calculation and a boundary motion file is input, the required motion boundary type and motion boundary correspond to the motion boundary condition, so that a moving grid area is selected through a Setup-dynamic mesh-Create/Edit-dynamic mesh command;
for a rigid motion boundary, selecting a RigidBody option, wherein motion attributes and grid options need to be explained;
for the deformed boundary, a Deforming option is selected, and a geometric definition and a grid option need to be explained.
According to the specific steps provided by the invention, the invention discloses the following technical effects: according to the detonation laser thruster analysis method based on FLUENT simulation, the change reason of the turning point of the thrust in the laser propelling process is discussed by utilizing FLUENT simulation software under the traditional parabolic detonation model, and corresponding gas can be adopted in the rocket launching process according to the change of the thrust in the laser propelling process, so that the impulse can be effectively prevented from rapidly dropping along with the rising of the rocket height, and the speed stability in the rocket flying process is ensured.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a flow chart of a flow field calculation of a detonation laser thruster analysis method based on FLUENT simulation;
FIG. 2 is a flow field pressure diagram of a thrust turning point of a detonation laser thruster analysis method based on FLUENT simulation;
FIG. 3 is a pressure diagram of a thrust turning point two-flow field of a detonation laser thruster analysis method based on FLUENT simulation;
FIG. 4 is a three-flow field pressure diagram of a thrust turning point of the detonation laser thruster analysis method based on FLUENT simulation.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a detonation laser thruster analysis method based on FLUENT simulation, which is characterized in that dynamic mesh modeling is applied, a user-defined function UDF is constructed, and the variation of thrust in the detonation laser propulsion process is researched by simulating by utilizing FLUENT software.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
The invention provides a knocking laser thruster analysis method based on FLUENT simulation, which comprises the following specific steps of:
s1, determining a control equation according to a laser wave absorption stage and a subsequent shock wave attenuation stage in the working process of the detonation mode;
the detonation mode working process is a typical abnormal process and comprises a laser absorption wave stage and a subsequent shock wave attenuation stage; for the laser absorption wave stage, due to the existence of the high-speed developed laser supporting detonation wave or combustion wave, thermodynamic non-equilibrium and chemical non-equilibrium factors need to be considered, at the moment, the control equation system keeps a complete form, and each flux also keeps unchanged:
Figure BDA0002876193800000051
for the subsequent shock wave attenuation stage, after the laser pulse is finished, the laser absorbs an energy source item, simultaneously, thermodynamic imbalance and chemical imbalance are reduced to a secondary position, at the moment, the gas flow state is chemical equilibrium and thermodynamic equilibrium, a control equation set is an axial symmetry flow equation considering the characteristics of high-temperature real gas, and each vector becomes:
Figure BDA0002876193800000052
Figure BDA0002876193800000053
Figure BDA0002876193800000054
the corresponding heat flow is also simplified to:
Figure BDA0002876193800000055
in the formula, U is a conservation variable vector, E, F is convection flux in x and r directions under a cylindrical coordinate system, Ev and Fv are respectively convection flux in x and r directions, G, Gv is a source item generated by adopting the cylindrical coordinate system for an unbonded part and a bonded part, W is a mass and energy source item caused by a physicochemical process, and q isxQr total and vibrational heat flows, qnCombining total heat conduction and diffusion heat flows of all internal energy modes for a heat flow q, wherein rho is the total density of the mixed gas, p is the pressure of the mixed gas, u and v are velocity components in two coordinate directions respectively, H is the enthalpy of the mixed gas with unit mass, and tau is the stress tensor;
s2, constructing a custom function UDF for laser energy setting, dynamic boundary setting and plasma parameter adjustment in the FLUENT calculation process, wherein the custom function UDF specifically comprises the following steps: a laser energy injection program, a motion boundary setting program and a plasma equation realization program;
UDF (user defined functions) is written in C language, can be loaded into FLUENT in a compiling or interpreting mode, and can be used for controlling some model parameters or calculation processes in the FLUENT calculation process; the UDF can realize more functions, and the laser energy injection, the dynamic boundary calculation and the plasma equation realization are realized by mainly utilizing the UDF;
laser energy injection procedure: DEFINE _ SOURCE in UDF is used to DEFINE various SOURCE terms in FLUENT's solution for transport equations, including: continuity, momentum equation, k, e, energy and component mass fractions, etc.; the invention mainly compiles the function to realize the injection of laser energy, and the specific format is as follows: DEFINE _ SOURCE (name, c, t, dS, eqn); the laser injection model mainly adopts an instantaneous energy injection model, namely, plasma is completely ionized and has no energy loss, laser pulse energy is completely absorbed by the plasma with certain deposition efficiency, the shape of the plasma is a sphere, and if the radius is r, a formula can be obtained:
Figure BDA0002876193800000061
here, the energy conversion efficiency is η 40%;
motion boundary setting program: when the floating grid computing is to be performed, a moving grid computing module needs to be selected, and a DEFINE _ CG _ MOTION macro in the UDF is used for writing and realizing functions, wherein the specific format is as follows: DEFINE _ CG _ MOTION (name, dt, vel, omega, time, dtime);
implementation procedure of plasma equation: the method comprises the following steps of modifying physical quantities of a flow field and calculating products by utilizing a general solution macro DEFINE _ ADJUST in FLUENT, wherein DEFINE _ ADJUST mainly integrates the whole flow field, then the boundary is adjusted by analyzing a calculation result, the macro defined by the function needs to be called for each iteration of each step or each solution of a transport equation, and the macro format is as follows: DEFINE _ ADJUST (name, d);
s3, establishing a parabolic model, modeling by using Gambit software, selecting a grid type to divide a grid, setting boundary conditions, and outputting a case file;
s4, opening FLUENT software, and importing the case file into the FLUENT software;
s5, detecting the grids, checking whether the divided grids are correct and setting the grid unit as mm;
s6, importing a custom function UDF by FLUENT software;
s7, selecting a correct radiation model, an energy model and a turbulence model, setting material properties, setting the flow field pressure to be 1atm and the temperature to be 298.15K;
s8, in FLUENT, performing dynamic grid calculation, selecting a dynamic grid calculation module, and setting the overall calculation parameters and modes of the dynamic grid;
in FLUENT, a dynamic mesh model is generally used for simulating the motion situation of a basin shape changing along with time caused by the motion of a basin boundary; generally this flow can be divided into two cases: the linear velocity or the angular velocity of the object changing along with time can move in a determined form, or an undetermined motion (the linear velocity or the angular velocity is obtained according to the stress balance of the gravity center of the object), the motion condition of each step is calculated according to the motion condition of the previous step, the body grid needs to be updated in each step, and the motion is automatically completed by FLUENT according to the new position of the boundary condition of each step;
in FLUENT, a dynamic grid computing module needs to be selected for dynamic grid computing, and overall computing parameters and modes of the dynamic grid are set; processing the propeller as a moving boundary by using a dynamic mesh module provided by FLUENT; when the calculation of each time step in the calculation is finished, calling the UDF to read flow field information, calculating the speed obtained by the propeller according to the force borne by the wall surface, and selecting a motion grid calculation and inputting a boundary motion file, wherein the type and the motion boundary of the motion boundary required correspond to motion boundary conditions, so that a motion grid area is selected through a Setup-dynamic mesh-Create/Edit-dynamic mesh command; for a rigid motion boundary, selecting a RigidBlock option to explain motion attributes and grid options, and for a deformation boundary, selecting a Deforming option to explain geometric definitions and grid options;
s9, selecting a solving method, and determining solving precision;
s10, initializing a boundary, setting iteration parameters, and calculating by using a FLUENT solver;
s11, checking a residual error map, a speed map, a pressure map and a temperature map;
as shown in fig. 2 to 4, the main reasons for analyzing the existence of three turning points are: the turning point is that the high-pressure air mass acts on the wall surface to achieve the maximum thrust, and the maximum thrust is caused by the high pressure behind the reflected wave after the shock wave reaches the top end of the thruster and then is reflected; the second turning point marks that the high-pressure air mass moves because the flow field changes gradually, and the low-pressure area begins to leave the top end of the thruster and moves to the side wall; the third turning point is that the shock wave gradually leaves the thruster, part of high-pressure air mass does not act on the thruster, and the thrust is gradually zero.
According to the detonation laser thruster analysis method based on FLUENT simulation, the change reason of the turning point of the thrust in the laser propelling process is discussed by utilizing FLUENT simulation software under the traditional parabolic detonation model, and corresponding gas can be adopted in the rocket launching process according to the change of the thrust in the laser propelling process, so that the impulse can be effectively prevented from rapidly dropping along with the rising of the rocket height, and the speed stability in the rocket flying process is ensured.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (5)

1. A detonation laser thruster analysis method based on FLUENT simulation is characterized by comprising the following steps:
s1, determining a control equation according to a laser wave absorption stage and a subsequent shock wave attenuation stage in the working process of the detonation mode;
s2, constructing a user-defined function UDF for laser energy setting, dynamic boundary setting and plasma parameter adjustment in the FLUENT calculation process;
s3, establishing a parabolic model, modeling by using Gambit software, selecting a grid type to divide a grid, setting boundary conditions, and outputting a case file;
s4, opening FLUENT software, and importing the case file into the FLUENT software;
s5, detecting the grids, checking whether the divided grids are correct and setting grid units;
s6, importing a custom function UDF by FLUENT software;
s7, selecting a correct radiation model, an energy model and a turbulence model, and setting material properties, flow field pressure and temperature;
s8, in FLUENT, performing dynamic grid calculation, selecting a dynamic grid calculation module, and setting the overall calculation parameters and modes of the dynamic grid;
s9, selecting a solving method, and determining solving precision;
s10, initializing a boundary, setting iteration parameters, and calculating by using a FLUENT solver;
s11, look at residual map, velocity, pressure and temperature map.
2. The method for analyzing the knocking laser thruster based on FLUENT simulation as claimed in claim 1, wherein the control equation in the step S1 specifically comprises:
Figure FDA0002876193790000011
Figure FDA0002876193790000012
Figure FDA0002876193790000013
Figure FDA0002876193790000014
the corresponding heat flow is also simplified to:
Figure FDA0002876193790000021
in the formula, U is a conservation variable vector, E, F is convection flux in x and r directions under a cylindrical coordinate system, Ev and Fv are respectively convection flux in x and r directions, G, Gv is a source item generated by adopting the cylindrical coordinate system for an unbonded part and a bonded part, W is a mass and energy source item caused by a physicochemical process, and q isx、qrTotal and vibrational heat flow, q, respectivelynAnd combining the total heat conduction and diffusion heat flow of each internal energy mode for the heat flow q, wherein rho is the total density of the mixed gas, p is the pressure of the mixed gas, u and v are velocity components in two coordinate directions respectively, H is the enthalpy of the mixed gas with unit mass, and tau is the stress tensor.
3. The method for analyzing a knocking laser thruster based on FLUENT simulation of claim 1, wherein the step S2 is to construct a custom function UDF for laser energy setting, dynamic boundary setting and plasma parameter adjustment in FLUENT calculation, which specifically comprises:
injecting laser energy: the instantaneous energy injection model is adopted, plasma ionization is complete and no energy loss exists, laser pulse energy is completely absorbed by the plasma with certain deposition efficiency, the shape of the plasma is a sphere, and if the radius is r, the formula is obtained:
Figure FDA0002876193790000022
wherein, the energy conversion efficiency is eta which is 40 percent;
calculating a moving boundary: when the flexible boundary calculation is carried out, a dynamic grid calculation module needs to be selected, and the DEFINE _ CG _ MOTION macro in the UDF is used for writing and realizing functions;
implementation of plasma equation: the method comprises the steps of modifying physical quantities of a flow field and calculating integration by utilizing a general solution macro DEFINE _ ADJUST in FLUENT, wherein DEFINE _ ADJUST mainly integrates the whole flow field, then the boundary is adjusted by analyzing a calculation result, and the macro defined by the function needs to be called for each iteration of each step or each solution of a transport equation.
4. The FLUENT simulation-based knocking laser thruster analysis method according to claim 1, wherein the grid unit in the step S5 is mm, the flow field pressure in the step S7 is 1atm, and the temperature is 298.15K.
5. The method for analyzing a knocking laser thruster based on FLUENT simulation of claim 1, wherein the step S8 is to perform dynamic grid calculation in FLUENT, select a dynamic grid calculation module, and set the overall calculation parameters and modes of the dynamic grid, specifically comprising:
processing the propeller as a Dynamic boundary by using a Dynamic Mesh module provided by FLUENT, calling UDF to read flow field information when the calculation of each time step is finished in the calculation, and calculating the speed obtained by the propeller according to the force borne by the wall surface;
when the moving grid is selected for calculation and a boundary motion file is input, the required motion boundary type and motion boundary correspond to the motion boundary condition, so that a moving grid area is selected through a Setup-Dynamic Mesh-Create/Edit-Dynamic Mesh command;
for a Rigid motion boundary, selecting a Rigid Body option, wherein motion attributes and grid options need to be explained;
for the deformed boundary, a Deforming option is selected, and a geometric definition and a grid option need to be explained.
CN202011621232.0A 2020-12-31 2020-12-31 Knock laser thruster analysis method based on FLUENT simulation Active CN112651083B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011621232.0A CN112651083B (en) 2020-12-31 2020-12-31 Knock laser thruster analysis method based on FLUENT simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011621232.0A CN112651083B (en) 2020-12-31 2020-12-31 Knock laser thruster analysis method based on FLUENT simulation

Publications (2)

Publication Number Publication Date
CN112651083A true CN112651083A (en) 2021-04-13
CN112651083B CN112651083B (en) 2023-08-11

Family

ID=75367367

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011621232.0A Active CN112651083B (en) 2020-12-31 2020-12-31 Knock laser thruster analysis method based on FLUENT simulation

Country Status (1)

Country Link
CN (1) CN112651083B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113468832A (en) * 2021-08-02 2021-10-01 中国空气动力研究与发展中心超高速空气动力研究所 Two-dimensional simulation method for shock tube diaphragm rupture process based on overlapped moving grids
CN113792498A (en) * 2021-08-26 2021-12-14 华南理工大学 FLUENT-based self-defined k-e turbulence model and initialization method
CN114036809A (en) * 2021-09-30 2022-02-11 东北电力大学 Method for predicting fly ash particle deposition based on dynamic grid and random function
CN117932794A (en) * 2024-03-25 2024-04-26 山东科技大学 Rotary tail structure optimization method based on disposable ocean current profiler probe

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107672824A (en) * 2017-09-30 2018-02-09 北华航天工业学院 Optical-fiber laser self-propelled universe bareboat control system
CN109858192A (en) * 2019-03-15 2019-06-07 中国人民解放军陆军装甲兵学院 Numerical simulation calculation method based on FLUENT MHD module
CN111199099A (en) * 2019-12-26 2020-05-26 兰州空间技术物理研究所 Method for evaluating operation life of ion thruster based on grid corrosion

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107672824A (en) * 2017-09-30 2018-02-09 北华航天工业学院 Optical-fiber laser self-propelled universe bareboat control system
CN109858192A (en) * 2019-03-15 2019-06-07 中国人民解放军陆军装甲兵学院 Numerical simulation calculation method based on FLUENT MHD module
CN111199099A (en) * 2019-12-26 2020-05-26 兰州空间技术物理研究所 Method for evaluating operation life of ion thruster based on grid corrosion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HIROSHI KATSURYAMA等: "Numerical analyses on pressure wave propagation in repetitive pulse laser propulsion", 《AIAA-2001-3665》, pages 1 - 9 *
许仁萍: "激光推进热力冲击破坏机理和防护研究", 《中国博士学位论文全文数据库 (工程科技Ⅱ辑)》, no. 6, pages 031 - 8 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113468832A (en) * 2021-08-02 2021-10-01 中国空气动力研究与发展中心超高速空气动力研究所 Two-dimensional simulation method for shock tube diaphragm rupture process based on overlapped moving grids
CN113468832B (en) * 2021-08-02 2023-02-24 中国空气动力研究与发展中心超高速空气动力研究所 Two-dimensional simulation method for shock tube diaphragm rupture process based on overlapped moving grids
CN113792498A (en) * 2021-08-26 2021-12-14 华南理工大学 FLUENT-based self-defined k-e turbulence model and initialization method
CN113792498B (en) * 2021-08-26 2024-03-29 华南理工大学 Method for customizing and initializing k-e turbulence model based on FLUENT
CN114036809A (en) * 2021-09-30 2022-02-11 东北电力大学 Method for predicting fly ash particle deposition based on dynamic grid and random function
CN117932794A (en) * 2024-03-25 2024-04-26 山东科技大学 Rotary tail structure optimization method based on disposable ocean current profiler probe

Also Published As

Publication number Publication date
CN112651083B (en) 2023-08-11

Similar Documents

Publication Publication Date Title
CN112651083A (en) Detonation laser thruster analysis method based on FLUENT simulation
Pulliam et al. Implicit finite-difference simulations of three-dimensional compressible flow
Yi et al. A three-dimensional numerical study of rotational detonation in an annular chamber
Sun et al. Influences of design parameters on a double serpentine convergent nozzle
Heo et al. Numerical study of the dynamic characteristics of pintle nozzles for variable thrust
Mo et al. Design of an asymmetric scramjet nozzle with circular to rectangular shape transition
CN112632709B (en) Continuous laser thruster working medium analysis method based on FLUENT simulation
Sreenath et al. Design and analysis of contour bell nozzle and comparison with dual bell nozzle
Schwane Numerical prediction and experimental validation of unsteady loads on ARIANE5 and VEGA
Haider et al. Mathematical analysis of flow passing through a rectangular nozzle
Cai et al. Detonation simulations in supersonic combustible mixtures with nonuniform species
Yumuşak Analysis and design optimization of solid rocket motors in viscous flows
Haselbacher et al. Open-ended shock tube flows: influence of pressure ratio and diaphragm position
Zhuang et al. Effectiveness of reaction control system in hypersonic rarefied reactive flow
Cummings et al. Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies
Li et al. Coupled Navier–Stokes/direct simulation Monte Carlo simulation of multicomponent mixture plume flows
Wu et al. Self-starting simulation of a hypersonic inlet with variable free stream condition
Folusiak et al. Numerical Modeling of the RDE
Haws et al. Computational investigation of a method to compress air fluidically in supersonic inlets
Jana et al. A numerical study of the flow field driven by a submerged, high-speed, gaseous jet
Xie et al. Fluidic nozzle throats in solid rocket motors
Zhu et al. Research Progress on Active Secondary Jet Technology in Supersonic Flow Field of Aerospace Propulsion Systems
Avital et al. Flow design and simulation of a gas compression system for hydrogen fusion energy production
CHEN et al. Numerical analysis of complex internal and external viscous flows with a second-order pressure-based method
Anand et al. 3D CFD analysis in an afterburner using NUMECA

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