JP5397634B2 - Thermal analysis method - Google Patents
Thermal analysis method Download PDFInfo
- Publication number
- JP5397634B2 JP5397634B2 JP2010215974A JP2010215974A JP5397634B2 JP 5397634 B2 JP5397634 B2 JP 5397634B2 JP 2010215974 A JP2010215974 A JP 2010215974A JP 2010215974 A JP2010215974 A JP 2010215974A JP 5397634 B2 JP5397634 B2 JP 5397634B2
- Authority
- JP
- Japan
- Prior art keywords
- exhaust pipe
- heat shield
- exhaust
- thermal analysis
- temperature
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- 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
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Description
本発明は、熱解析方法に関し、特に、自動車の排気管の温度分布を解析する方法に関する。 The present invention relates to a thermal analysis method, and more particularly to a method of analyzing a temperature distribution of an exhaust pipe of an automobile.
従来、自動車の開発においては、CAE(Computer Aided Engineering)システムが用いられている(例えば、特許文献1参照)。このようなシステムは、自動車の熱の挙動を計算する機能を有しており、その計算結果に基づいて、熱の挙動を考慮した自動車の設計を行うことができる。 Conventionally, in the development of automobiles, a CAE (Computer Aided Engineering) system has been used (see, for example, Patent Document 1). Such a system has a function of calculating the thermal behavior of the automobile, and based on the calculation result, the automobile can be designed in consideration of the thermal behavior.
特に、高温となる排気管の温度分布は、その周辺部分の遮熱対策などを設計するために重要な情報であるので、排気管の温度分布を精度よく計算することは、重要な課題となっている。 In particular, the temperature distribution of the exhaust pipe, which is at a high temperature, is important information for designing measures for heat shielding in the surrounding area. Therefore, it is important to accurately calculate the temperature distribution of the exhaust pipe. ing.
排気管の温度を推定する方法としては、コンピュータを用いた熱流体解析が知られている。熱流体解析では、計算対象の排気管をモデル化し、所定の形状、大きさの計算格子に分割する。そして、計算格子に対して、排気管に流れる排気ガスの温度や粘度や密度などの流体物性、排気管の材質の熱伝導率等の物性値などを設定し、コンピュータで解析することで、計算格子における温度を求めることができる。 As a method for estimating the temperature of the exhaust pipe, thermal fluid analysis using a computer is known. In thermal fluid analysis, an exhaust pipe to be calculated is modeled and divided into calculation grids of a predetermined shape and size. Then, set the physical properties such as the temperature, viscosity and density of the exhaust gas flowing in the exhaust pipe, the physical properties such as the thermal conductivity of the exhaust pipe material, etc. to the calculation grid, and analyze it with a computer. The temperature at the grid can be determined.
排気管には、その周囲に遮熱板が設けられる場合があり、この遮熱板を考慮して排気管の温度分布を計算しなければならない。すなわち、対流や輻射より車体の熱が遮熱板に伝わり、さらに、遮熱板からその内側の排気管に対流や輻射により熱が伝わるので、このことを考慮した計算が必要となる。したがって、排気管と遮熱板との間の空間にも計算格子を設定し、その空間での対流や輻射を考慮した解析を行う必要がある。 A heat shield plate may be provided around the exhaust pipe, and the temperature distribution of the exhaust pipe must be calculated in consideration of the heat shield plate. That is, the heat of the vehicle body is transmitted to the heat shield plate by convection and radiation, and further, the heat is transmitted from the heat shield plate to the exhaust pipe inside thereof by convection and radiation. Therefore, it is necessary to set a calculation grid also in the space between the exhaust pipe and the heat shield, and perform analysis in consideration of convection and radiation in the space.
しかし、排気管と遮熱板との間は、狭いところで数ミリ程度と狭く、この空間を計算格子に分割することは困難である。また、その空間に多量の計算格子を設定して解析すると、解析に要する時間が膨大となり、収束しない虞がある。 However, the space between the exhaust pipe and the heat shield is as narrow as several millimeters in a narrow space, and it is difficult to divide this space into calculation grids. Also, if a large amount of calculation grid is set in the space and analyzed, the time required for the analysis becomes enormous and may not converge.
本発明は、上記事情に鑑み、遮熱板が配設された排気管の温度分布を高速かつ高精度に解析することができる熱解析方法を提供することを目的とする。 In view of the above circumstances, an object of the present invention is to provide a thermal analysis method capable of analyzing a temperature distribution of an exhaust pipe provided with a heat shield plate at high speed and with high accuracy.
上記目的を達成するための本発明の第1の態様は、排気ガスが流通する排気管及び該排気管の周囲に空間を隔てて配設された遮熱板を備える排気系のモデルを用いて前記排気管の温度分布を解析する熱解析方法であって、前記排気管及びその内部の流路には計算格子を設定すると共に、前記排気管と前記遮熱板との間の空間部には計算格子を設定せずに前記排気系のモデルに計算格子を設定するステップと、前記遮熱板から前記排気管への輻射と、前記空間部の対流伝熱とにより前記排気管に伝熱することを表す境界条件を設定し、前記排気管の温度分布の熱解析を行うステップとをコンピュータにより実行させることを特徴とする熱解析方法にある。 In order to achieve the above object, a first aspect of the present invention uses an exhaust system model including an exhaust pipe through which exhaust gas flows and a heat shield disposed around the exhaust pipe with a space therebetween. A thermal analysis method for analyzing a temperature distribution of the exhaust pipe, wherein a calculation grid is set for the exhaust pipe and a flow path inside the exhaust pipe, and in a space portion between the exhaust pipe and the heat shield plate Heat is transferred to the exhaust pipe by setting a calculation grid in the exhaust system model without setting a calculation grid, radiation from the heat shield plate to the exhaust pipe, and convective heat transfer in the space. And a step of performing a thermal analysis of the temperature distribution of the exhaust pipe by a computer .
かかる第1の態様では、遮熱板と排気管との間の空間部に計算格子を設けずに熱解析を行うため、計算格子を設ける時間が短縮され、かつ熱解析に要する計算時間が短縮される。さらに、遮熱板から排気管への輻射の影響と、空間部での対流の影響とを考慮した境界条件を用いるので、実測の温度に近似した排気管及び遮熱板の温度分布を得ることができる。 In the first aspect, since the thermal analysis is performed without providing the calculation grid in the space between the heat shield and the exhaust pipe, the time for providing the calculation grid is reduced and the calculation time required for the thermal analysis is also reduced. Is done. In addition, since boundary conditions that take into account the effects of radiation from the heat shield to the exhaust pipe and the effects of convection in the space are used, the temperature distribution of the exhaust pipe and heat shield that approximates the measured temperature can be obtained. Can do.
本発明の第2の態様は、第1の態様に記載する熱解析方法において、前記排気系を備える車体のモデルを用い、前記車体及び当該車体から前記遮熱板の外表面までの空間に計算格子を設定して熱解析を行うことで、前記車体から前記遮熱板への入射輻射量、前記遮熱板の外表面側の流体温度及び熱伝達係数を解析するステップと、これらの入射輻射量、流体温度及び熱伝達係数を含む前記境界条件に基づいて熱解析を行うステップとをコンピュータにより実行させることを特徴とする熱解析方法にある。 According to a second aspect of the present invention, in the thermal analysis method described in the first aspect, calculation is performed on a space from the vehicle body and the vehicle body to the outer surface of the heat shield plate using a vehicle body model including the exhaust system. set the grid by performing thermal analysis, the incident amount of radiation to the heat shield plate from the vehicle body, and a step of analyzing the fluid temperature and heat transfer coefficient of the outer surface of the heat shield plate, these incoming radiation And a step of performing a thermal analysis based on the boundary conditions including a quantity, a fluid temperature and a heat transfer coefficient by a computer .
かかる第2の態様では、車体から遮熱板及び排気管に伝わる熱の影響を解析する際に、計算格子を設定する範囲を遮熱板の外表面までとしたので、複雑な遮熱板の内側に計算格子を設定する必要がない。このため、車体のモデルに計算格子を作成する時間と、解析に要する時間を大幅に短縮することができる。 In the second aspect, when analyzing the influence of the heat transmitted from the vehicle body to the heat shield and the exhaust pipe, the range for setting the calculation grid is extended to the outer surface of the heat shield. There is no need to set a calculation grid inside. For this reason, it is possible to greatly reduce the time required to create a calculation grid in the vehicle body model and the time required for analysis.
本発明の第3の態様は、第2の態様に記載する熱解析方法において、前記遮熱板が前記排気管を覆っていない部分については、前記入射輻射量、前記流体温度、前記熱伝達係数及び前記排気管の内表面の温度から構成された境界条件に基づいて熱解析を行うステップをコンピュータにより実行させることを特徴とする熱解析方法にある。 According to a third aspect of the present invention, in the thermal analysis method according to the second aspect, the incident radiation amount, the fluid temperature, and the heat transfer coefficient for a portion where the heat shield plate does not cover the exhaust pipe. And a step of performing a thermal analysis on the basis of a boundary condition constituted by the temperature of the inner surface of the exhaust pipe.
かかる第3の態様では、排気管の全長に亘って遮熱板が覆っていない場合においても、排気管の温度分布を解析することができる。
In the third aspect, the temperature distribution of the exhaust pipe can be analyzed even when the heat shield does not cover the entire length of the exhaust pipe.
本発明によれば、遮熱板が配設された排気管の温度分布を高速かつ高精度に解析することができる熱解析方法が提供される。 ADVANTAGE OF THE INVENTION According to this invention, the thermal analysis method which can analyze the temperature distribution of the exhaust pipe in which the heat shield was arrange | positioned at high speed and with high precision is provided.
本実施形態に係る熱解析方法の対象となる排気管及び遮熱板について説明する。
図1は、自動車の排気系を示す概略構成図であり、図2は、排気管及び遮熱板の一部断面を拡大した図である。
An exhaust pipe and a heat shield that are targets of the thermal analysis method according to the present embodiment will be described.
FIG. 1 is a schematic configuration diagram showing an exhaust system of an automobile, and FIG. 2 is an enlarged view of a partial cross section of an exhaust pipe and a heat shield plate.
図1に示すように、自動車の排気系1は、エンジン(図示せず)からの排気ガスの流路となる排気管2を備えている。排気管2には、排気ガスを浄化する排気浄化装置4が設けられ、さらに、排気浄化装置4の下流には、マフラー5が設けられている。エンジンから排出された排気ガスは、排気管2を介して排気浄化装置4に送られ、浄化される。さらに浄化された排気ガスは、排気管2を介してマフラー5に送られ、減音され外部に放出される。
As shown in FIG. 1, an
図2に示すように、排気管2は、円筒状に形成されており、その内表面2a側に排気ガスが流れるようになっている。また、排気管2の周囲には、空間を隔てて遮熱板3が配設されている。具体的には、遮熱板3は、排気管2よりも径の大きい円筒状又は円筒状の一部に形成されており、その内側(内表面3a側)に排気管2が配設されている。また、排気管2と遮熱板3との間には、空間部6があり、空間部6には流体(空気)が存在している。
As shown in FIG. 2, the
本実施形態では、遮熱板3は、排気管2の全長に亘って設けられてはおらず、排気管2の外表面2bの一部が外部に露出している(この露出した部分を排気管露出部7と称する)。すなわち、排気系1としては、排気管露出部7と遮熱板3の外表面3bとが外部に露出している。
In the present embodiment, the
図3は、車体及び遮熱板から排気管に伝わる熱の挙動を説明する図である。図示するように、排気管2の表面の温度は、その内部の流路8を流れる排気ガスに影響を受ける。また、遮熱板3には、車体9で生じた熱が輻射及び対流により伝わり、排気管2には、その遮熱板3からの熱が輻射及び空間部6での対流により伝わる。すなわち、排気管2の温度分布を解析するためには、排気管2を流通する排気ガスだけでなく、遮熱板3からの熱を考慮する必要がある。
FIG. 3 is a diagram for explaining the behavior of heat transmitted from the vehicle body and the heat shield to the exhaust pipe. As shown in the figure, the temperature of the surface of the
ここで、上述した構造の排気管2を対象とし、本実施形態に係る熱解析方法を詳細に説明する。図4は、本実施形態に係る熱解析方法の処理手順を示すフロー図である。各ステップは、コンピュータによって実行される。コンピュータは、演算を行うCPU、演算中のデータなどの各種データやプログラムを記憶させておくRAMやROMなどの記憶装置、キーボードやマウスなどの入力装置、演算結果などを表示するモニタなどの出力装置を備えている。
Here, the thermal analysis method according to the present embodiment will be described in detail with the
まず、車体9の熱解析をすることで、車体9から遮熱板3への入射輻射量と、遮熱板3表面の流体の温度及び熱伝達係数を計算する(ステップS1)。これらの入射輻射量、流体の温度及び熱伝達係数は、排気管2の温度分布を解析する際に用いる境界条件に必要とされるデータである。
First, by performing thermal analysis of the
これらのデータは、具体的には、排気系1を含む車体9全体のモデルを形成し、モデル化された車体9を計算格子に分割し、各種の解析条件を設定して、熱解析を行うことにより得られる。
Specifically, these data form a model of the
排気系1を含む車体9のモデルは、車体9の形状を表した数値データであり、3次元CADシステムなどにより作成されたものである。この車体9のモデルを計算格子に分割する。
The model of the
図5に、車体のモデルに設定した計算格子を例示する。図示するように、車体9の熱解析においては、車体9及び遮熱板3の外表面3bまでの間の空間に計算格子10を設定する。なお、特に図示しないが、排気管露出部7(図2参照)については、車体9から排気管2の外表面2bの露出した部分までの空間に計算格子10を設定する。
FIG. 5 shows an example of a calculation grid set in the vehicle body model. As shown in the figure, in the thermal analysis of the
そして、排気系1及び車体9が計算範囲に入るように境界条件を設定し、車体9の表面から遮熱板3の外表面3bへの輻射、及び車体9と遮熱板3の外表面3bとの間での対流とによる伝熱を解析する。なお、排気管露出部7に関しては、車体9の表面からの排気管2の外表面2bへの輻射、及び車体9と排気管2との間での対流による伝熱を解析する。
Then, boundary conditions are set so that the
この解析により、車体9から遮熱板3(及び排気管露出部7における排気管2の外表面2b)への入射輻射量と、車体9と遮熱板3(及び排気管露出部7における排気管2の外表面2b)との間における流体の温度及び熱伝達係数が得られる。
By this analysis, the amount of incident radiation from the
車体9の熱解析においては、計算格子10を設定する範囲を遮熱板3の外表面までとしたので、複雑な遮熱板3の内側に計算格子10を設定する必要がない。このため、車体9のモデルに計算格子10を作成する時間と、解析に要する時間を大幅に短縮することができる。
In the thermal analysis of the
次に、排気管2の外表面2bの温度を計算する。具体的には、図1及び図2に示した排気系1のモデルを形成する(図4ステップS2)。このモデルは、排気系1の形状を表した数値データであり、3次元CADシステムなどにより作成されたものである。
Next, the temperature of the
次に、モデル化された排気系1に計算格子10を設定する(図4ステップS3)。図6に、排気系に設定した計算格子10を例示する。図示するように、排気系1の熱解析においては、排気管2及びその内部の流路8に計算格子10を設定し、排気管2の外表面2bから遮熱板3の内表面3aの間の空間部6には計算格子10を設定しない。
Next, the
このように計算格子10を設定することで、狭い空間部6に多数の計算格子10を設定することに要する時間を短縮することができる。また、空間部6に計算格子10を設定しないので排気系1の熱解析に要する時間を短縮することができ、また、解析が収束しないということを回避することができる。
By setting the
次に、計算格子10を設定した排気系1に、境界条件を設定して熱解析を実行する(図4ステップS4)。排気管2の遮熱板3に覆われている部分については、遮熱板3から排気管2の外表面2bへの輻射により伝熱し、かつ遮熱板3と排気管2との間の空間部6における対流により排気管2に伝熱することを表す境界条件を設定する。
Next, boundary conditions are set in the
具体的には、排気管2の遮熱板3に覆われている部分については、次の式(1)〜(6)を境界条件として用いる。
Specifically, the following formulas (1) to (6) are used as boundary conditions for the portion of the
ただし、ε1:排気管2外表面2bの輻射放射率、ε2:遮熱板3の輻射放射率、σ:ボルツマン定数、G1:排気管2の外表面2bの射度、G2:遮熱板3の射度、Tw0:排気管2の外表面2bの温度、Tw:遮熱板3の温度(厚さは無視している)、Ti:排気管2の内表面2aの温度、λ:排気管2の熱伝導率、d:排気管2の壁面の厚さ、h1:空間部6の熱伝達係数、Ta:空間部6の温度、Fx,x:輻射形態係数、I:車体9から遮熱板3への入射輻射量、Tf:遮熱板3の外表面3b側の流体の温度、h2:遮熱板3の外表面3b側の流体の熱伝達係数、α、β:係数である。
Where ε 1 : radiation emissivity of the
また、排気管2の遮熱板3に覆われていない排気管露出部7については、次の式(7)を境界条件として用いる。
For the exhaust pipe exposed
これらの境界条件に加えて、入口境界として、排気管2の上流のエンジン(排気マニホールド)出口の排気ガスの流量、密度、温度などの物性値を用い、出口境界として、マフラー5から外部に排出される排気ガスの圧力を用いる。本実施形態では、エキゾーストマニホールドについて熱解析を行うことにより排気ガスの流量、密度、温度の物性値を計算し、これらを入口境界として用いた。また、入口境界は、エキゾーストマニホールドの出口であり、排気ガスの非定常性の影響が弱くなることから、定常計算で排気管2の表面温度を解析することとした。
In addition to these boundary conditions, exhaust gas flow rate, density, temperature, and other physical property values at the outlet of the engine (exhaust manifold) upstream of the
そして、上述した各種条件を用いると共に、ステップS1で算出した入射輻射量I、遮熱板3の外表面3b側の流体の温度Ta、熱伝達係数h2を読み込み、かつ、排気管2に流れる排気ガスの熱解析をして排気管の内表面2aの温度Tiを計算し、上記式(1)〜(6)の連立方程式を解いて排気管2の外表面2bの温度Tw0と遮熱板3の温度Twを計算する。なお、式(1)〜(6)の解法の一例は、本明細書の最後に示す。
And while using the various conditions mentioned above, the incident radiation amount I calculated in step S1, the temperature T a of the fluid on the
排気管露出部7については、同様にして、入射輻射量Iを読み込み、かつ、排気管2に流れる排気ガスの熱解析をして排気管の内表面2aの温度Tiを計算し、上記式(7)を解いて排気管2の外表面2bの温度Tw0と遮熱板3の温度Twを計算する。
For the exhaust pipe exposed
なお、本実施形態では、排気浄化装置4についてもモデル化している。排気浄化装置4は、排気ガスから窒素酸化物等を除去する触媒が設けられており、触媒の発熱を解析してある。
In the present embodiment, the
このようにして計算格子10ごとに得られた排気管2の外表面2bの温度Tw0と、遮熱板3の温度Twを車体9の熱解析の排気系境界に与え、車体各部品の温度を算出することができる。これにより、各部品の温度を考慮した自動車の設計を行うことができる。
And thus calculated temperature of the
以上に説明したように、本実施形態に係る熱解析方法では、遮熱板3と排気管2との間の空間部6に計算格子10を設けずに熱解析を行うため、計算格子10を設ける時間が短縮され、かつ熱解析に要する計算時間が短縮される。
As described above, in the thermal analysis method according to the present embodiment, since the thermal analysis is performed without providing the
さらに、排気管2の温度を解析する際に、遮熱板3から排気管2への輻射の影響と、空間部6での対流の影響とを考慮した境界条件を用いたので、実測の温度に近似した排気管2及び遮熱板3の温度分布を得ることができる。
Furthermore, when analyzing the temperature of the
また、排気管2には、遮熱板3が外側に配設され、また排気管2の一部は遮熱板3に覆われずに外部に露出した部分がある排気系1について熱解析を行ったが、これに限らず、排気管2は全長に亘って遮熱板3に覆われた排気系1であってもよい。また、排気系1は、排気管2、遮熱板3、排気浄化装置4、及びマフラー5から構成されていたが、これに限らず、少なくとも排気管2(排気管2のみの場合は式(7)を利用する)から構成されていればよい。
Further, the
さらに、上述の実施形態では、自動車の排気系1を対象としたが、これに限らず、自動車以外の排気系についても本発明の熱解析方法を適用することができる。
Furthermore, in the above-described embodiment, the
本実施形態に係る熱解析方法で得られた遮熱板3の温度及び排気管2の外表面2bの温度と、実測の遮熱板3の温度及び排気管2の外表面2bの温度とを比較した。
The temperature of the
解析対象のモデルは、図1に示した排気系1を対象とし、空間部6を除いた排気系1全体を約53万の計算格子に分割した。
The model to be analyzed is the
解析条件の一例として次のものを用いた。
解析タイプ :定常解析
乱流モデル :k−ε/RNG乱流モデル
入口境界 :排気菅の入口における排気ガスの流量、密度、温度
出口境界 :排気菅の出口における排気ガスの圧力
境界条件 :壁面熱流速条件
排気ガス物性 :圧縮性がある熱空気
密度 :温度と圧力の関数
粘性 :定数、4.61E-05kg/ms
比熱 :定数、995J/kgK
熱伝導率 :温度の関数
Con = 1.5207E-11 × T3 - 4.8574E-08 ×T2
+1.0184E-04 ×T - 3.9333E-04
分子量 :26.68kg/kmol
触媒 :ポーラス、体積発熱
発熱量 :Q=L×Cp×ΔT=0.0481×995×38=1.82KW
Lはガス流量(kg/s)、Wave04にて算出
Cpは定圧比熱(J/kgK)
ΔTは出入口の排気ガス温度差(K)
The following was used as an example of analysis conditions.
Analysis type: Steady analysis Turbulence model: k-ε / RNG turbulence model Inlet boundary: Exhaust gas flow rate, density and temperature at the exhaust soot Outlet boundary: Exhaust gas pressure at the exhaust soot Boundary condition: Wall heat Flow rate condition Exhaust gas physical property: Hot air with compressibility Density: Function of temperature and pressure Viscosity: Constant, 4.61E-05kg / ms
Specific heat: Constant, 995J / kgK
Thermal conductivity: function of temperature
Con = 1.5207E-11 × T3-4.8574E-08 × T2
+ 1.0184E-04 × T-3.9333E-04
Molecular weight: 26.68kg / kmol
Catalyst: Porous, Volumetric heat generation Calorific value: Q = L × Cp × ΔT = 0.0481 × 995 × 38 = 1.82KW
L is gas flow rate (kg / s), calculated by Wave04
Cp is constant pressure specific heat (J / kgK)
ΔT is the exhaust gas temperature difference at the inlet and outlet (K)
図7に、解析結果を示すグラフを示す。横軸は、排気管2の上流から順に適宜選択した部位の番号であり、縦軸は温度である。同図には、各部位の排気ガスの温度と、遮熱板3及び排気管2の外表面2bの温度の折れ線グラフが示されている。折れ線グラフの各部位にある四角印は、本実施形態に係る熱解析方法で得られた遮熱板3及び排気管2の外表面2bの温度(計算値)を表し、×印は、実測の遮熱板3及び排気管2の外表面2bの温度(実測値)を表している。
FIG. 7 shows a graph showing the analysis results. The horizontal axis is the number of the part appropriately selected from the upstream of the
排気管2の上流(横軸左側)から下流(横軸右側)までの全体に亘り、計算値は実測値との誤差が小さく、非常に精度が高いことが示されている。
From the upstream (left side of the horizontal axis) to the downstream (right side of the horizontal axis) of the
[式(1)〜(6)の解法例]
上述した境界条件方程式(1)〜(6)の解法を説明する。
ここで、輻射形態係数Fx,xは、排気管2を「1」、遮熱板3を「2」とすると、F1,1=0、F1,2=1、F2,1=A2/A1、F2,2=(A2−A1)/A2である。遮熱板3と排気管2の距離が近くA2≒A1とすると、F2,1=1、F2,2=2となる。
[Example of solving equations (1) to (6)]
A method for solving the above-described boundary condition equations (1) to (6) will be described.
Here, when the
式(1)〜(6)のG1とG2とを消去すると、次の式(8)、(9)が得られる。
式(8)、(9)に示されるように、未知の変数はTw0とTwの2つである。ここで、 Equation (8), (9) as shown in, the unknown variables are two of T w0 and T w. here,
式(11)について、ニュートン反復法を用いてTwを計算し、そのTwを式(10)に代入する。式(10)の左辺>右辺であるならば、Tw0を大きくしてTwを再計算し、左辺<右辺であるならば、Tw0を小さくしてTwを再計算する。左辺と右辺の差が所定の精度以下になるまで繰り返し計算することで、Tw0とTwが算出できる。 For formula (11), to calculate the T w using Newton iterative method, and substitutes the T w in equation (10). If the left side of the equation (10)> the right side, T w0 is increased and T w is recalculated. If the left side <the right side, T w0 is reduced and T w is recalculated. T w0 and T w can be calculated by repeatedly calculating until the difference between the left side and the right side becomes a predetermined accuracy or less.
本発明の熱解析方法は、自動車の排気系の表面温度、特に排気管の表面温度の解析の他、各種交通機関や産業機械の部品であって内部に排気ガスが流れる排気系部品の表面温度を解析する場合に適用できる。 The thermal analysis method of the present invention is applicable to the analysis of the surface temperature of the exhaust system of an automobile, particularly the surface temperature of the exhaust pipe, as well as the surface temperature of exhaust system parts through which exhaust gas flows, which are parts of various transportation and industrial machines. Applicable when analyzing
1 排気系
2 排気管
3 遮熱板
4 排気浄化装置
5 マフラー
6 空間部
7 排気管露出部
8 流路
9 車体
DESCRIPTION OF
Claims (3)
前記排気管及びその内部の流路には計算格子を設定すると共に、前記排気管と前記遮熱板との間の空間部には計算格子を設定せずに前記排気系のモデルに計算格子を設定するステップと、
前記遮熱板から前記排気管への輻射と、前記空間部の対流伝熱とにより前記排気管に伝熱することを表す境界条件を設定し、前記排気管の温度分布の熱解析を行うステップとをコンピュータにより実行させることを特徴とする熱解析方法。 A thermal analysis method for analyzing a temperature distribution of the exhaust pipe using an exhaust system model including an exhaust pipe through which exhaust gas flows and a heat shield disposed around the exhaust pipe with a space therebetween,
A calculation grid is set for the exhaust pipe and the flow path inside the exhaust pipe, and a calculation grid is set in the exhaust system model without setting a calculation grid in the space between the exhaust pipe and the heat shield plate. Steps to set,
Step of setting the radiation from the heat shield plate to the exhaust pipe, a boundary condition indicating that heat transfer to the exhaust pipe by the convective heat transfer of the space, for thermal analysis of the temperature distribution of the exhaust pipe And a computer-executable thermal analysis method.
前記排気系を備える車体のモデルを用い、
前記車体及び当該車体から前記遮熱板の外表面までの空間に計算格子を設定して熱解析を行うことで、前記車体から前記遮熱板への入射輻射量、前記遮熱板の外表面側の流体温度及び熱伝達係数を解析するステップと、
これらの入射輻射量、流体温度及び熱伝達係数を含む前記境界条件に基づいて熱解析を行うステップとをコンピュータにより実行させることを特徴とする熱解析方法。 In the thermal analysis method according to claim 1 ,
Using a model of the car body with the exhaust system,
The amount of incident radiation from the vehicle body to the heat shield plate and the outer surface of the heat shield plate are determined by setting a calculation grid in the vehicle body and a space from the vehicle body to the outer surface of the heat shield plate. a step of analyzing the fluid temperature and heat transfer coefficient of the side,
And a step of performing a thermal analysis based on the boundary condition including the incident radiation amount, the fluid temperature, and the heat transfer coefficient by a computer .
前記遮熱板が前記排気管を覆っていない部分については、
前記入射輻射量、前記流体温度、前記熱伝達係数及び前記排気管の内表面の温度から構成された境界条件に基づいて熱解析を行うステップをコンピュータにより実行させることを特徴とする熱解析方法。 In the thermal analysis method according to claim 2 ,
For the part where the heat shield plate does not cover the exhaust pipe,
A thermal analysis method, comprising: causing a computer to execute a step of performing a thermal analysis based on a boundary condition configured from the incident radiation amount, the fluid temperature, the heat transfer coefficient, and the temperature of the inner surface of the exhaust pipe.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010215974A JP5397634B2 (en) | 2010-09-27 | 2010-09-27 | Thermal analysis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010215974A JP5397634B2 (en) | 2010-09-27 | 2010-09-27 | Thermal analysis method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012073670A JP2012073670A (en) | 2012-04-12 |
JP5397634B2 true JP5397634B2 (en) | 2014-01-22 |
Family
ID=46169806
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010215974A Expired - Fee Related JP5397634B2 (en) | 2010-09-27 | 2010-09-27 | Thermal analysis method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5397634B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109145388B (en) * | 2018-07-25 | 2023-02-28 | 中国航发沈阳发动机研究所 | Thermal analysis method for aircraft engine component |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002032419A (en) * | 2000-07-19 | 2002-01-31 | Mitsubishi Heavy Ind Ltd | Management system for piping stress |
JP2002183223A (en) * | 2000-12-15 | 2002-06-28 | Sharp Corp | Device and method for processing mesh and computer- readable recording medium with mesh processing program recorded thereon |
JP4389843B2 (en) * | 2005-06-14 | 2009-12-24 | パナソニック株式会社 | Analysis mesh model generation device, analysis mesh model generation method, analysis mesh model generation program, and recording medium |
JP4967752B2 (en) * | 2007-03-28 | 2012-07-04 | マツダ株式会社 | Vehicle design support system |
-
2010
- 2010-09-27 JP JP2010215974A patent/JP5397634B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2012073670A (en) | 2012-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lange et al. | Impact of manufacturing variability and nonaxisymmetry on high-pressure compressor stage performance | |
Maffulli et al. | Impact of wall temperature on heat transfer coefficient and aerodynamics for three-dimensional turbine blade passage | |
Zou et al. | Numerical modeling of steady state errors for shielded thermocouples based on conjugate heat transfer analysis | |
Montomoli et al. | Uncertainty quantification and conjugate heat transfer: a stochastic analysis | |
Sharma et al. | Evaluation of modelling parameters for computing flow-induced noise in a small high-speed centrifugal compressor | |
Kumar et al. | Underhood thermal simulation of a small passenger vehicle with rear engine compartment to evaluate and enhance radiator performance | |
Shahane et al. | Finite volume simulation framework for die casting with uncertainty quantification | |
Musch et al. | Optimization strategy for a coupled design of the last stage and the successive diffuser in a low pressure steam turbine | |
Winchler et al. | Conjugate heat transfer methodology for thermal design and verification of gas turbine cooled components | |
Duda et al. | Solution of an inverse transient heat conduction problem in a part of a complex domain | |
Suzuki et al. | Experimental study of graphite ablation in nitrogen flow, part II: Further numerical analysis | |
Lavagnoli et al. | Analysis of the heat transfer driving parameters in tight rotor blade tip clearances | |
JP2006284214A (en) | Thermal analyzing method and its program | |
Cozzi et al. | Capturing radial mixing in axial compressors with computational fluid dynamics | |
Rahaim et al. | Coupled dual reciprocity boundary element/finite volume method for transient conjugate heat transfer | |
El-Sharkawy et al. | Transient modelling of vehicle exhaust surface temperature | |
JP5397634B2 (en) | Thermal analysis method | |
Da Soghe et al. | Development of experimental and numerical methods for the analysis of active clearance control systems | |
VanderVeer et al. | Optimization of an inverse convection solution strategy | |
VanderVeer et al. | Solution of an inverse convection problem by a predictor–corrector approach | |
Mick et al. | Study on relevant effects concerning heat transfer of a convection cooled gas turbine blade under realistic engine temperature conditions | |
Mathison et al. | Temperature Predictions and Comparison With Measurements for the Blade Leading Edge and Platform of a 1 1/2 Stage Transonic HP Turbine | |
Haehndel et al. | The development of exhaust surface temperature models for 3D CFD vehicle thermal management simulations Part 1-general exhaust configurations | |
Fu et al. | Uncertainty propagation through a simulation of industrial high pressure die casting | |
Choi et al. | A metamodeling approach for uncertainty analysis of nondeterministic systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120824 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130501 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130515 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130625 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20130925 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20131008 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5397634 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |