CN106886666A - The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio under a kind of different angles of attack - Google Patents
The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio under a kind of different angles of attack Download PDFInfo
- Publication number
- CN106886666A CN106886666A CN201710221791.4A CN201710221791A CN106886666A CN 106886666 A CN106886666 A CN 106886666A CN 201710221791 A CN201710221791 A CN 201710221791A CN 106886666 A CN106886666 A CN 106886666A
- Authority
- CN
- China
- Prior art keywords
- underwater glider
- attack
- lift
- drag ratio
- underwater
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention proposes a kind of body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio under different angles of attack, purpose is to determine that the optimal of underwater glider slides angle, and resistance and lift under the different angles of attack can be obtained, for the selection of the control system and power source of follow-up underwater glider provides theory support.The accuracy of body fitted anisotropic mesh numerical simulation finally by experimental verification, the identical situation that obtains is preferable.This analysis method be can be generalized under other Series of Water in exercise machine, and calculating cycle is saved compared with traditional complicated formulas are derived, and computational efficiency can be improved more than twice, the problems such as more effectively predicting and solve the efficiency of underwater.Step A, underwater glider physical model simplify.Step B, the outer domain model of establishment underwater glider fluid.Step C, establishment underwater glider boundary-layer fence lattice model.Step D, the numerical solution calculating for carrying out underwater glider lift-drag ratio under the different angles of attack.Step E, carry out data result post processing.The present invention is applied to effective solution of underwater glider lift-drag ratio.
Description
Technical field
The invention belongs to target property and identification technology field, underwater glider lift-drag ratio under a kind of different angles of attack is relate to
Body fitted anisotropic mesh method for solving.
Background technology
Underwater glider is a kind of new under water without cable robot, and its power source is the net buoyancy of seawater, by net
Buoyancy underwater glider can be slided slowly along jagged track in the seawater.Because it does not have push system, tool
There is noise small, continuation of the journey is long, the features such as coasting distance is long.And various sensors can be installed thereon as auxiliary, underwater gliding
Machine plays vital effect at aspects such as the collection of ocean big data, military investigation, marine environment and monitoring resources.In recent years
Come, underwater glider numerical simulation is mainly the calculating of lift-drag ratio, bigger lift, and smaller resistance and bigger lift-drag ratio can
Aerodone is brought preferably continuation of the journey, improve its hydrodynamic efficiency.Therefore, the present invention proposes to be slided under water under a kind of different angles of attack
The method of value solving of Xiang machine lift-drag ratio, it is therefore an objective to determine that the optimal of underwater glider slides angle, and the different angles of attack can be obtained
Under resistance and lift, for the selection of the control system and power source of follow-up underwater glider provides theory support.
The content of the invention
It is poor to the accuracy of underwater glider glide efficiency problem solving the present invention is to solve existing computational methods
Present situation, a kind of body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio under different angles of attack is provided for this.
The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio is realized according to the following steps under a kind of different angles of attack：
Step A, underwater glider physical model simplify：
First have to simplify model in ICEM, carried out repairing in simplifying on tail vane, eliminated on fluid emulation influence
Hook and communication system before less release, as shown in figure.
Step B, establishment underwater glider fluid foreign lands are as follows：
Cylindrical fluid domain is set up to underwater glider using ICEM mesh generations software, it is underwater gliding head that fluid domain is anterior
The similar semi-circular structure in portion, as shown in Figure 2, and carries out mesh generation to it, destructuring net of the generation with body fitted anisotropic mesh
Lattice, carry out mesh refinement treatment, as shown in Figure 3 to boundary layer.
Step C, the numerical simulation condition for determining underwater glider：
Underwater glider grid is imported into Fluent, the turbulence model being derived from according to early stage, and underwater gliding
The real work situation of machine, is respectively provided with boundary condition for entrance boundary condition, export boundary condition and wall perimeter strip
Part, sets seawater fluid medium, and carrying out conventional -10 ° to 10 ° of angles of attack to underwater glider respectively carries out Numerical Simulation.
And the interative computation of 6000 steps is carried out, and monitor is set, analog simulation obtains flow field data.
Step D, the numerical solution calculating for carrying out underwater glider lift-drag ratio under the different angles of attack.
Step E, data result post processing is carried out, the lift that numerical analysis is obtained is collected with drag consequence, is obtained
Lift-drag ratio, fits curve of the lift-drag ratio with angle of attack variation, draws contrast curve, and its optimal glide angle is drawn by contrast
Degree.
Invention effect：
The present invention obtains the optimal glide efficiency of underwater glider by the method for liquid numerical simulation, has drawn liter under the different angles of attack
Resistance ratio, is contrasted by the analogue value and experiment value, and situation of coincideing preferably, demonstrates the accuracy of method of value solving.This analysis side
Method can be generalized under other Series of Water in exercise machine, be derived compared to saving calculating cycle with traditional complicated formulas, meter
Calculating efficiency can be improved more than twice, and manpower and materials are saved compared with traditional experiment, more economically, more can effectively be predicted and be asked
The problems such as solving the efficiency of underwater.
Brief description of the drawings
Fig. 1 is underwater glider rough schematic view；
Fig. 2 is underwater glider fluid domain schematic diagram；
Fig. 3 is body fitted anisotropic mesh schematic diagram；
Fig. 4 is pressure field distribution.
Specific embodiment
The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio is realized according to the following steps under a kind of different angles of attack：
The present invention is main research meanses with Computer Simulation, is broadly divided into numerical modeling and CFD emulation, wherein CFD emulation
It is divided into mesh generation again（Pre-treatment）, calculate and solve, and interpretation of result（Post processing）；
Step A, underwater glider model simplification：
First have to simplify model in ICEM, the simplification of model can not only improve mesh quality, improve the essence for calculating
Degree, is an essential link in simulation analysis while also greatly reduce the time of calculating, and tail vane is entered in simplifying
Gone repairing, eliminate release before hook and communication system.
Step B, establishment underwater glider fluid foreign lands are as follows：
Because underwater glider does relative motion with Flow Field outside, it is assumed that underwater glider is motionless, fluid is with relative underwater glider
The speed flowing slided.Cylindrical fluid domain is set up to underwater glider using ICEM mesh generations software, fluid domain front portion is
The similar semi-circular structure of underwater gliding chain-drive section, and mesh generation is carried out to it, destructuring of the generation with body fitted anisotropic mesh
Grid, mesh refinement treatment is carried out to boundary layer, as shown in Figure 3, number of grid 2,100,000, nodes 450,000.
Step C, the numerical simulation condition for determining underwater glider：
It is assumed that condition is：Seawater is incompressible fluid, and the influence of wave is not considered.Density of sea water and temperature keep constant, i.e.,
At 15 DEG C, density isρ=1025.0kg/m3, dynamic viscosity coefficient isμ=1.219×10-3kg/m·s-1.Do not consider further that the earth
Surface curvature and earth rotation, seawater is considered as plane.Omit seawater pressure and make the elasticity of underwater glider housing generation
Deformation, is considered as underwater glider the rigid body of constancy of volume.
The turbulence model being derived from according to early stage, and underwater glider real work situation, be respectively provided with border
Condition is entrance boundary condition, export boundary condition and wall boundary condition, seawater fluid medium is set, respectively under water
Aerodone carries out conventional -10 ° to 10 ° of angles of attack and carries out Numerical Simulation.
Step D, the numerical solution calculating for carrying out underwater glider lift-drag ratio under the different angles of attack.
Step D1, opening fluid calculation software, grid model file is imported by the Import in File
Shuixiahuaxiangji. mesh, and mesh quality inspection is carried out by Grid；
Step D2, the parameter setting that solver Solver is carried out by Define, are flowed using Three Dimensional Steady fluid and set；
Step D3, material properties Materials settings are carried out by Define, mainly provide the viscosity and density of studied seawater
Parameter etc.；
Step D4, model M odels settings are carried out by Define, selection is carriedRNGk-εEquations of Turbulence, Wall-function for mark
Quasi- Wall-function, remaining keeps acquiescence；
Step D5, the setting of Define-Boundary Conditions boundary conditions is carried out, selection speed entrance velocity-
Inlet, when the angle of attack is 0 °, and speed is 0.5m/s, it is as shown in the table for specific setting.Difference according to the angle of attack is come appropriate regulation Y
The speed of axle, in gliding, the angle of attack is usually no more than 10 ° to aerodone, its Y-axis component is calculated, while bringing Y- into
Velocity（m/s）, exit selection free discharge outflow, wall is quiet for fricton-tight wall wall ROVs surface is
Only without sliding wall wall, and surface roughness is 0；
Step D6, carry out monitor plane setting；
Step D7, residual error are set to 10-5Initialization flow field, is iterated calculating, and flow field numerical value meter is pointed out after the step number of iteration 6000
Convergence is calculated, result of calculation is obtained.
Step E, data result post processing is carried out, the lift that numerical analysis is obtained is collected with drag consequence, is obtained
Lift-drag ratio, fits curve of the lift-drag ratio with angle of attack variation, draws contrast curve, and its optimal glide angle is drawn by contrast
Degree, it is as shown in the table, can draw this type underwater glider angle of attack be 8 ° when, lift-drag ratio for maximum, now there is highest to slide effect
Rate.
Present embodiment effect：
Analyzed by numerical result of the invention, there is highest lift-drag ratio when the angle of attack is 8 °, now underwater gliding
Machine glide efficiency highest, and the accuracy of body fitted anisotropic mesh numerical simulation by experimental verification, the identical situation that obtains are preferable.
This analysis method be can be generalized under other Series of Water in exercise machine, and calculating is saved compared with traditional complicated formulas are derived
Cycle, computational efficiency can be improved more than twice, and manpower and materials have been saved compared with traditional experiment, more can effectively be predicted and be solved
The problems such as efficiency of underwater.
Claims (2)
1. under a kind of different angles of attack underwater glider lift-drag ratio body fitted anisotropic mesh method for solving, it is characterised in that by boundary layer
The method that grid is solved obtains underwater glider Flow Field outside distribution situation, has drawn the liter resistance of the underwater glider under the different angles of attack
Than, contrasted by the analogue value and test value, verify the accuracy of method of value solving；
This analysis method can be generalized under other Series of Water in exercise machine, be saved compared with traditional complicated formulas are derived
Calculating cycle, computational efficiency can be improved more than twice, the problems such as more effectively predicting and solve underwater glider performance；
It is poor to the accuracy and efficiency of underwater glider lift-drag ratio problem solving the present invention is to solve existing computational methods
Present situation, for this provides underwater glider lift-drag ratio body fitted anisotropic mesh method for solving under a kind of different angles of attack；
The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio is realized according to the following steps under a kind of different angles of attack：
Step A, slide under water machine physical model simplification：
First have to simplify model, carried out repairing in simplifying on tail vane, eliminated release less on fluid emulation influence
Hook and communication system before putting；
Step B, establishment underwater glider boundary-layer fence lattice model：
Cylindrical fluid domain is set up to underwater glider, it is the similar semi-circular structure of underwater gliding chain-drive section that fluid domain is anterior, and
Mesh generation is carried out to it, unstructured grid of the generation with body fitted anisotropic mesh carries out mesh refinement treatment to boundary layer；
Step C, the numerical solution condition for determining underwater glider：
The turbulence model being derived from according to early stage, and underwater glider real work situation, be respectively provided with boundary condition,
Seawater fluid medium is set, and carrying out conventional -10 ° to 10 ° of angles of attack to underwater glider respectively carries out body fitted anisotropic mesh pattern number
Value is solved；
Step D, the 6000 step iterative computings for carrying out underwater glider lift-drag ratio under the different angles of attack；
Step E, data result post processing is carried out, fit body fitted anisotropic mesh and solve the lift-drag ratio parameter for obtaining with angle of attack variation
Curve, show that it most preferably slides angle.
2. under a kind of different angles of attack according to claim 1 underwater glider lift-drag ratio body fitted anisotropic mesh method for solving,
It is characterized in that with the increase of the angle of attack, lift coefficient and resistance coefficient increase therewith, when it is 8 ° to slide angle, resistance is now risen
Than being maximum, glide efficiency highest, this is consistent with the conclusion that experiment draws, demonstrates the standard of body fitted anisotropic mesh method of value solving
True property；
This analysis method can be generalized under other Series of Water in exercise machine, be saved compared with traditional complicated formulas are derived
Calculating cycle, computational efficiency can be improved more than twice, more can effectively be predicted and be solved hydrodynamic performance of underwater etc. and ask
Topic.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710221791.4A CN106886666B (en) | 2017-04-06 | 2017-04-06 | Boundary layer grid solving method for lift-drag ratio of underwater glider under different attack angles |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710221791.4A CN106886666B (en) | 2017-04-06 | 2017-04-06 | Boundary layer grid solving method for lift-drag ratio of underwater glider under different attack angles |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106886666A true CN106886666A (en) | 2017-06-23 |
| CN106886666B CN106886666B (en) | 2019-12-31 |
Family
ID=59181883
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710221791.4A Expired - Fee Related CN106886666B (en) | 2017-04-06 | 2017-04-06 | Boundary layer grid solving method for lift-drag ratio of underwater glider under different attack angles |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106886666B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109948269A (en) * | 2019-03-27 | 2019-06-28 | 湖南科技大学 | A kind of armouring umbilical cables model hexahedral mesh division methods |
| CN110276131A (en) * | 2019-06-24 | 2019-09-24 | 西北工业大学 | Blended wing-body underwater glider profile optimization method based on polynomial response surface model |
| CN110309571A (en) * | 2019-06-24 | 2019-10-08 | 西北工业大学 | Blended wing-body underwater glider profile optimization method based on radial basis function model |
| CN112947502A (en) * | 2021-03-04 | 2021-06-11 | 中国科学院自动化研究所 | Flexible bionic web underwater motion control method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060058985A1 (en) * | 2004-08-31 | 2006-03-16 | Supersonic Aerospace International, Llc | Adjoint-based design variable adaptation |
| CN101246515A (en) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | Digital ship model planar motion mechanism experimental method based on CFD software |
| CN102486633A (en) * | 2010-12-04 | 2012-06-06 | 中国科学院沈阳自动化研究所 | Motion parameter optimization method for optimal energy consumption of underwater glider |
| CN104102783A (en) * | 2014-07-18 | 2014-10-15 | 安徽工业大学 | Method for forecasting numerical value of cavitation of underwater propeller tip vortex |
| CN105912763A (en) * | 2016-04-06 | 2016-08-31 | 西安交通大学 | Heat conduction path topological optimization based underwater glider cruise path planning method |
-
2017
- 2017-04-06 CN CN201710221791.4A patent/CN106886666B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060058985A1 (en) * | 2004-08-31 | 2006-03-16 | Supersonic Aerospace International, Llc | Adjoint-based design variable adaptation |
| CN101246515A (en) * | 2008-03-03 | 2008-08-20 | 哈尔滨工程大学 | Digital ship model planar motion mechanism experimental method based on CFD software |
| CN102486633A (en) * | 2010-12-04 | 2012-06-06 | 中国科学院沈阳自动化研究所 | Motion parameter optimization method for optimal energy consumption of underwater glider |
| CN104102783A (en) * | 2014-07-18 | 2014-10-15 | 安徽工业大学 | Method for forecasting numerical value of cavitation of underwater propeller tip vortex |
| CN105912763A (en) * | 2016-04-06 | 2016-08-31 | 西安交通大学 | Heat conduction path topological optimization based underwater glider cruise path planning method |
Non-Patent Citations (1)
| Title |
|---|
| 刘国良: ""水下滑翔机水动力性能模拟及实验研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109948269A (en) * | 2019-03-27 | 2019-06-28 | 湖南科技大学 | A kind of armouring umbilical cables model hexahedral mesh division methods |
| CN109948269B (en) * | 2019-03-27 | 2020-11-03 | 湖南科技大学 | Hexahedral mesh division method for armored umbilical cable model |
| CN110276131A (en) * | 2019-06-24 | 2019-09-24 | 西北工业大学 | Blended wing-body underwater glider profile optimization method based on polynomial response surface model |
| CN110309571A (en) * | 2019-06-24 | 2019-10-08 | 西北工业大学 | Blended wing-body underwater glider profile optimization method based on radial basis function model |
| CN110309571B (en) * | 2019-06-24 | 2022-02-11 | 西北工业大学 | Wing body fusion underwater glider external shape optimization method based on radial basis function model |
| CN112947502A (en) * | 2021-03-04 | 2021-06-11 | 中国科学院自动化研究所 | Flexible bionic web underwater motion control method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106886666B (en) | 2019-12-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106886666A (en) | The body fitted anisotropic mesh method for solving of underwater glider lift-drag ratio under a kind of different angles of attack | |
| Batten et al. | Accuracy of the actuator disc-RANS approach for predicting the performance and wake of tidal turbines | |
| CN108319793B (en) | CFD-based analysis method for free movement of low-speed inclined water inlet | |
| CN104091085B (en) | The cavitation noise feature assessment method being calculated based on screw current pressure fluctuation | |
| US9399981B2 (en) | Methods and systems for comparing vertical axis turbine arrays and providing configurations thereof | |
| D'Alessandro et al. | Assessment of the dimples as passive boundary layer control technique for laminar airfoils operating at wind turbine blades root region typical Reynolds numbers | |
| CN101246515A (en) | Digital ship model planar motion mechanism experimental method based on CFD software | |
| CN109992889B (en) | Wind power plant model building method and system and wake value calculating method and system | |
| Hu et al. | Tip vortex prediction for contra-rotating propeller using large eddy simulation | |
| Liao et al. | A simulation-based actuator surface parameterization for large-eddy simulation of propeller wakes | |
| Alidadi et al. | A numerical method for calculation of power output from ducted vertical axis hydro-current turbines | |
| Peng et al. | Numerical study of cavitation on the surface of the guide vane in three gorges hydropower unit | |
| Shi et al. | Broadband force spectra of a pump-jet propulsor subjected to inflow turbulence: Comparison with ducted propeller and propeller | |
| CN113139306A (en) | Numerical prediction method for cavitation noise of composite propeller | |
| Li et al. | Propulsive properties of a flexible oscillating wing with time-varying camber deformation | |
| CN110489902A (en) | A kind of propeller cavitation wake flow fine-feature multivariate statistics modeling method | |
| CN115879231B (en) | Method and device for acquiring flow transition information of micro aircraft | |
| Rathod et al. | Evolving a Bio-Inspired Blade Shape of the Drag-Based Vertical-Axis Wind Rotor Derived From Orange Sea Pen (Ptilosarcus Gurneyi) | |
| Yamazaki | Experiment/simulation integrated shape optimization using variable fidelity Kriging model approach | |
| McNae | Unsteady hydrodynamics of tidal stream turbines | |
| Yeo et al. | Analysis of the effect of vortex reduction structures on submarine tonal noise via frequency-domain method employing thickness noise source | |
| Zhu et al. | Research Progress on Identification and Suppression Methods for Monitoring the Cavitation State of Centrifugal Pumps | |
| CN113392600B (en) | Relevance analysis method for propeller low-frequency line spectrum characteristics in non-uniform inflow | |
| Bruzzone et al. | Rudder-propeller interaction: analysis of different approximation techniques | |
| CN110032818B (en) | Borderless swimming pool water drop flow calculation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yu Xiaodong Inventor after: Quan Zhen Inventor after: Liu Guoliang Inventor after: Kong Xiangbin Inventor after: Zhang Yanqin Inventor after: Cheng Haikuo Inventor after: Guo Lili Inventor before: Zhang Yanqin Inventor before: Quan Zhen Inventor before: Liu Guoliang Inventor before: Kong Xiangbin Inventor before: Yu Xiaodong Inventor before: Cheng Haikuo Inventor before: Guo Lili |
|
| CB03 | Change of inventor or designer information | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191231 Termination date: 20210406 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |