燃燒仿真軟件CFD++的前處理技術(shù)教程

燃燒仿真軟件CFD++的前處理技術(shù)教程1燃燒仿真的基礎(chǔ)理論1.1燃燒反應(yīng)機理燃燒是一種復(fù)雜的化學(xué)反應(yīng)過程，涉及到燃料與氧化劑（通常是空氣中的氧氣）的快速氧化反應(yīng)，產(chǎn)生熱能和光能。在燃燒仿真中，理解燃燒反應(yīng)機理是至關(guān)重要的，因為它直接影響到燃燒模型的選擇和仿真結(jié)果的準(zhǔn)確性。1.1.1燃燒反應(yīng)類型燃燒反應(yīng)可以分為幾種類型，包括：均相燃燒：燃料和氧化劑在分子水平上混合，如氣體燃燒。非均相燃燒：燃料和氧化劑在不同相態(tài)下反應(yīng)，如液體燃料的燃燒。1.1.2燃燒反應(yīng)方程燃燒反應(yīng)方程描述了燃料和氧化劑反應(yīng)生成產(chǎn)物的化學(xué)過程。例如，甲烷（CH4）在氧氣（O2）中的燃燒反應(yīng)方程可以表示為：CH4+2O2->CO2+2H2O1.1.3燃燒反應(yīng)速率燃燒反應(yīng)速率受多種因素影響，包括溫度、壓力、反應(yīng)物濃度和催化劑的存在。在CFD仿真中，通常使用Arrhenius方程來描述燃燒反應(yīng)速率：r=A*exp(-Ea/RT)*[fuel]^m*[O2]^n其中，r是反應(yīng)速率，A是頻率因子，Ea是活化能，R是氣體常數(shù)，T是溫度，[fuel]和[O2]分別是燃料和氧氣的濃度，m和n是反應(yīng)級數(shù)。1.2湍流模型與燃燒模型在燃燒仿真中，湍流模型和燃燒模型是兩個關(guān)鍵的組成部分，它們幫助我們理解和預(yù)測燃燒過程中的流體動力學(xué)和化學(xué)反應(yīng)行為。1.2.1湍流模型湍流模型用于描述流體中的湍流行為，常見的湍流模型包括：k-ε模型：基于湍流動能（k）和湍流耗散率（ε）的模型。k-ω模型：基于湍流動能（k）和渦旋頻率（ω）的模型。1.2.2燃燒模型燃燒模型用于描述燃燒過程中的化學(xué)反應(yīng)，常見的燃燒模型包括：層流火焰模型：適用于層流燃燒，假設(shè)火焰?zhèn)鞑ニ俣仁浅?shù)。PDF模型：概率密度函數(shù)模型，適用于湍流燃燒，考慮了湍流對燃燒反應(yīng)的影響。1.3網(wǎng)格生成原理網(wǎng)格生成是CFD仿真的基礎(chǔ)，它將連續(xù)的物理空間離散化為一系列有限的單元，以便進行數(shù)值計算。1.3.1網(wǎng)格類型網(wǎng)格可以分為幾種類型，包括：結(jié)構(gòu)網(wǎng)格：網(wǎng)格單元按照規(guī)則排列，如矩形或六面體網(wǎng)格。非結(jié)構(gòu)網(wǎng)格：網(wǎng)格單元不規(guī)則排列，適用于復(fù)雜幾何形狀的仿真。1.3.2網(wǎng)格質(zhì)量網(wǎng)格質(zhì)量直接影響到仿真的準(zhǔn)確性和計算效率。網(wǎng)格質(zhì)量的評估通常包括：網(wǎng)格單元的形狀：單元應(yīng)盡可能接近正方形或正六面體。網(wǎng)格單元的大小：在高梯度區(qū)域（如火焰附近）應(yīng)使用更小的單元。1.3.3網(wǎng)格生成工具CFD++軟件通常支持多種網(wǎng)格生成工具，包括：Gambit：用于生成結(jié)構(gòu)和非結(jié)構(gòu)網(wǎng)格。TetGen：專門用于生成非結(jié)構(gòu)三維網(wǎng)格。1.3.4示例：使用Gambit生成非結(jié)構(gòu)網(wǎng)格#使用Gambit生成非結(jié)構(gòu)網(wǎng)格的示例命令

gambit-batch-inputinput.gmb-outputoutput.msh在上述命令中，input.gmb是Gambit的輸入文件，包含了幾何形狀和網(wǎng)格生成的參數(shù)；output.msh是輸出的網(wǎng)格文件，可以被CFD++軟件讀取并用于仿真。1.3.5示例：網(wǎng)格質(zhì)量評估#使用Python評估網(wǎng)格質(zhì)量的示例代碼

importnumpyasnp

#假設(shè)我們有網(wǎng)格單元的體積和面積數(shù)據(jù)

volumes=np.array([1.0,1.2,1.1,1.3,1.4])

areas=np.array([0.5,0.6,0.55,0.65,0.7])

#計算網(wǎng)格單元的形狀質(zhì)量

aspect_ratios=volumes/areas

#輸出網(wǎng)格單元的平均形狀質(zhì)量

print("平均形狀質(zhì)量:",np.mean(aspect_ratios))在這個示例中，我們使用了Python的numpy庫來計算網(wǎng)格單元的形狀質(zhì)量。volumes和areas數(shù)組分別包含了網(wǎng)格單元的體積和面積數(shù)據(jù)。通過計算體積與面積的比值，我們可以得到網(wǎng)格單元的形狀質(zhì)量，進而評估網(wǎng)格的整體質(zhì)量。通過以上內(nèi)容，我們深入了解了燃燒仿真的基礎(chǔ)理論，包括燃燒反應(yīng)機理、湍流模型與燃燒模型以及網(wǎng)格生成原理。這些理論知識是進行燃燒仿真前處理技術(shù)的關(guān)鍵，能夠幫助我們更準(zhǔn)確地設(shè)置仿真參數(shù)，提高仿真的可靠性和效率。2CFD++軟件介紹2.1軟件功能與應(yīng)用領(lǐng)域CFD++是一款高性能的計算流體動力學(xué)(CFD)軟件，廣泛應(yīng)用于航空航天、汽車、能源、化工等多個領(lǐng)域。其核心功能包括：流體動力學(xué)模擬：能夠處理復(fù)雜的流體流動問題，包括湍流、層流、可壓縮流和不可壓縮流。熱力學(xué)分析：支持熱傳導(dǎo)、對流和輻射的計算，適用于燃燒、傳熱等過程的仿真?；瘜W(xué)反應(yīng)模擬：能夠模擬多種化學(xué)反應(yīng)，特別適用于燃燒仿真，包括預(yù)混燃燒、擴散燃燒等。多物理場耦合：可以同時模擬流體、熱力學(xué)和化學(xué)反應(yīng)，實現(xiàn)多物理場的耦合分析。網(wǎng)格生成與處理：內(nèi)置強大的網(wǎng)格生成工具，支持結(jié)構(gòu)化、非結(jié)構(gòu)化和混合網(wǎng)格的生成與優(yōu)化。2.1.1應(yīng)用實例在汽車工業(yè)中，CFD++被用于優(yōu)化發(fā)動機燃燒過程，提高燃油效率和減少排放。例如，通過模擬燃燒室內(nèi)燃料的噴射、混合和燃燒過程，工程師可以設(shè)計更高效的燃燒室形狀和燃料噴射策略。2.2用戶界面與操作流程CFD++提供了一個直觀的用戶界面，使用戶能夠輕松地設(shè)置和運行仿真。操作流程大致如下：前處理：定義幾何模型，生成網(wǎng)格，設(shè)置邊界條件和物理模型。求解設(shè)置：選擇求解器，設(shè)置求解參數(shù)，如時間步長、迭代次數(shù)等。求解運行：運行仿真，軟件將根據(jù)設(shè)定的物理模型和邊界條件求解流體動力學(xué)方程。后處理：分析仿真結(jié)果，可視化流場、溫度分布、化學(xué)反應(yīng)等數(shù)據(jù)。2.2.1前處理技術(shù)詳解2.2.1.1幾何模型定義在前處理階段，首先需要定義幾何模型。CFD++支持導(dǎo)入多種格式的幾何模型，包括STL、IGES、STEP等。對于復(fù)雜的幾何結(jié)構(gòu)，軟件提供了強大的幾何編輯工具，允許用戶進行切割、合并、旋轉(zhuǎn)等操作。2.2.1.2網(wǎng)格生成網(wǎng)格生成是CFD仿真中的關(guān)鍵步驟。CFD++提供了多種網(wǎng)格生成技術(shù)，包括：結(jié)構(gòu)化網(wǎng)格：適用于規(guī)則幾何形狀，網(wǎng)格單元為正方形或立方體。非結(jié)構(gòu)化網(wǎng)格：適用于復(fù)雜幾何形狀，網(wǎng)格單元可以是三角形、四邊形、四面體或六面體?；旌暇W(wǎng)格：結(jié)合結(jié)構(gòu)化和非結(jié)構(gòu)化網(wǎng)格的優(yōu)點，適用于既有規(guī)則部分又有復(fù)雜部分的幾何模型。2.2.1.3設(shè)置邊界條件邊界條件的設(shè)置直接影響仿真的準(zhǔn)確性和可靠性。CFD++支持多種邊界條件，包括：入口邊界條件：可以設(shè)置速度、壓力、溫度和化學(xué)組分等。出口邊界條件：通常設(shè)置為自由出口或壓力出口。壁面邊界條件：可以設(shè)置為絕熱壁面、恒溫壁面或熱流壁面。2.2.1.4物理模型選擇根據(jù)仿真需求，用戶需要選擇合適的物理模型。例如，在燃燒仿真中，可能需要選擇：湍流模型：如k-ε模型、k-ω模型或雷諾應(yīng)力模型。燃燒模型：如預(yù)混燃燒模型、擴散燃燒模型或PDF模型。化學(xué)反應(yīng)模型：定義化學(xué)反應(yīng)方程式和反應(yīng)速率。2.2.2示例：網(wǎng)格生成與邊界條件設(shè)置#使用CFD++的前處理工具生成非結(jié)構(gòu)化網(wǎng)格

#假設(shè)我們有一個汽車發(fā)動機燃燒室的STL模型

#下面的命令將生成一個非結(jié)構(gòu)化網(wǎng)格，并保存為CGNS格式

grid_generatorinput.stl--nonstructured--output=engine.cgns

#設(shè)置邊界條件

#假設(shè)我們有三個邊界：入口、出口和壁面

#入口設(shè)置為速度入口，出口為壓力出口，壁面為絕熱壁面

boundary_conditionengine.cgns--inlet=velocity--outlet=pressure--wall=adiabatic在上述示例中，我們首先使用grid_generator命令生成了一個非結(jié)構(gòu)化網(wǎng)格，并將其保存為CGNS格式。然后，我們使用boundary_condition命令設(shè)置了邊界條件，入口設(shè)置為速度入口，出口設(shè)置為壓力出口，壁面設(shè)置為絕熱壁面。通過這些步驟，我們可以為CFD++準(zhǔn)備一個完整的前處理文件，用于后續(xù)的流體動力學(xué)和燃燒仿真。以上內(nèi)容詳細介紹了CFD++軟件的功能、應(yīng)用領(lǐng)域以及用戶界面和操作流程，特別是前處理技術(shù)中的幾何模型定義、網(wǎng)格生成、邊界條件設(shè)置和物理模型選擇。通過一個具體的示例，展示了如何使用CFD++的前處理工具生成非結(jié)構(gòu)化網(wǎng)格，并設(shè)置邊界條件，為燃燒仿真準(zhǔn)備必要的前處理文件。3前處理技術(shù)詳解3.1幾何模型的創(chuàng)建與導(dǎo)入在進行燃燒仿真之前，首先需要創(chuàng)建或?qū)霂缀文Ｐ?。幾何模型是燃燒仿真中的基礎(chǔ)，它定義了燃燒室的形狀、尺寸以及內(nèi)部結(jié)構(gòu)。CFD++軟件支持多種幾何模型的導(dǎo)入格式，包括但不限于STL、IGES、STEP等。3.1.1創(chuàng)建幾何模型幾何模型的創(chuàng)建通常在CAD軟件中完成，如SolidWorks、AutoCAD或CATIA。在創(chuàng)建模型時，需要確保模型的準(zhǔn)確性和細節(jié)，以反映真實的燃燒室結(jié)構(gòu)。3.1.2導(dǎo)入幾何模型一旦幾何模型創(chuàng)建完成，接下來的步驟是將其導(dǎo)入到CFD++中。這通常通過軟件的導(dǎo)入功能完成，選擇相應(yīng)的文件格式，確保模型正確無誤地導(dǎo)入。3.2網(wǎng)格劃分技術(shù)網(wǎng)格劃分是將幾何模型分割成一系列小的、離散的單元，以便進行數(shù)值計算。網(wǎng)格的質(zhì)量直接影響仿真的準(zhǔn)確性和計算效率。3.2.1結(jié)構(gòu)化網(wǎng)格生成結(jié)構(gòu)化網(wǎng)格生成通常用于形狀規(guī)則的區(qū)域，如圓柱、矩形等。這種網(wǎng)格的特點是每個單元的形狀和大小相對一致，網(wǎng)格點在空間中呈有規(guī)律的分布。3.2.1.1示例代碼#使用Python的OpenFOAM工具包生成結(jié)構(gòu)化網(wǎng)格

fromfoamToolsimportstructuredMesh

#定義幾何參數(shù)

length=1.0

width=0.5

height=0.2

nCellsX=100

nCellsY=50

nCellsZ=20

#生成結(jié)構(gòu)化網(wǎng)格

mesh=structuredMesh(length,width,height,nCellsX,nCellsY,nCellsZ)

mesh.writeMesh("structuredMesh")在CFD++中，雖然直接使用Python代碼生成網(wǎng)格可能不適用，但上述示例展示了結(jié)構(gòu)化網(wǎng)格生成的基本思路。3.2.2非結(jié)構(gòu)化網(wǎng)格生成非結(jié)構(gòu)化網(wǎng)格生成適用于復(fù)雜幾何形狀，如燃燒室內(nèi)部的不規(guī)則結(jié)構(gòu)。這種網(wǎng)格的單元形狀和大小可以變化，以適應(yīng)幾何的復(fù)雜性。3.2.2.1示例代碼#使用Gmsh生成非結(jié)構(gòu)化網(wǎng)格

importgmsh

#初始化Gmsh

gmsh.initialize()

#創(chuàng)建模型

model=gmsh.model

model.add("nonStructuredMesh")

#導(dǎo)入幾何模型

model.occ.importShapes("complexGeometry.stl")

#設(shè)置網(wǎng)格參數(shù)

model.mesh.setSize(model.occ.getEntities(0),0.01)

#生成網(wǎng)格

model.mesh.generate(3)

#導(dǎo)出網(wǎng)格

gmsh.write("nonStructuredMesh.msh")

#關(guān)閉Gmsh

gmsh.finalize()Gmsh生成的非結(jié)構(gòu)化網(wǎng)格可以轉(zhuǎn)換為CFD++可讀的格式，如CGNS或VTK，以便在軟件中使用。3.3邊界條件設(shè)置邊界條件是燃燒仿真中不可或缺的一部分，它定義了仿真域與外部環(huán)境的交互。3.3.1入口邊界條件入口邊界條件通常包括速度、溫度和化學(xué)組分的初始值。這些條件對于模擬燃燒過程至關(guān)重要。3.3.1.1示例在入口邊界，設(shè)置速度為10m/s，溫度為300K，氧氣濃度為21%。3.3.2出口邊界條件出口邊界條件通常設(shè)定為大氣壓力或自由出流條件，以模擬燃燒產(chǎn)物的排放。3.3.2.1示例在出口邊界，設(shè)定壓力為1atm，確保燃燒產(chǎn)物可以自由流出。3.3.3壁面邊界條件壁面邊界條件用于模擬燃燒室壁面的熱傳遞和化學(xué)反應(yīng)。常見的壁面條件包括絕熱壁面和指定溫度的壁面。3.3.3.1示例在燃燒室壁面，設(shè)定為絕熱壁面，以模擬壁面無熱傳遞的情況。以上內(nèi)容詳細介紹了在CFD++軟件中進行燃燒仿真前處理技術(shù)的關(guān)鍵步驟，包括幾何模型的創(chuàng)建與導(dǎo)入、網(wǎng)格劃分技術(shù)以及邊界條件的設(shè)置。通過這些步驟，可以為后續(xù)的燃燒仿真提供準(zhǔn)確的幾何和物理條件，確保仿真的可靠性和準(zhǔn)確性。4燃燒模型的設(shè)置4.1化學(xué)反應(yīng)模型化學(xué)反應(yīng)模型在燃燒仿真中至關(guān)重要，它描述了燃料與氧化劑之間的化學(xué)反應(yīng)過程。在CFD++軟件中，化學(xué)反應(yīng)模型的設(shè)置通常涉及反應(yīng)機理的選擇、反應(yīng)速率的計算以及化學(xué)物種的輸運。4.1.1選擇反應(yīng)機理反應(yīng)機理的選擇依賴于所研究的燃料類型和燃燒條件。例如，對于甲烷燃燒，可以使用GRI3.0或GRI1.2機理，這些機理包含了詳細的化學(xué)反應(yīng)路徑和物種。4.1.2反應(yīng)速率計算反應(yīng)速率的計算基于Arrhenius定律，其公式為：r其中，r是反應(yīng)速率，A是頻率因子，Ea是活化能，R是氣體常數(shù)，T4.1.3化學(xué)物種輸運在CFD++中，化學(xué)物種的輸運通過擴散方程來描述，這需要設(shè)置擴散系數(shù)和混合模型。4.2湍流燃燒模型湍流燃燒模型考慮了湍流對燃燒過程的影響，這對于理解實際燃燒環(huán)境中的火焰?zhèn)鞑ズ腿紵手陵P(guān)重要。4.2.1湍流模型選擇CFD++提供了多種湍流模型，如k-ε模型、k-ω模型和雷諾應(yīng)力模型。在燃燒仿真中，通常選擇k-ε模型或k-ω模型，因為它們在計算成本和準(zhǔn)確性之間提供了良好的平衡。4.2.2火焰?zhèn)鞑ツＰ突鹧鎮(zhèn)鞑ツＰ兔枋隽送牧魅绾斡绊懟鹧娴膫鞑ニ俣?。常見的模型包括EDC（EddyDissipationConcept）和PDF（ProbabilityDensityFunction）模型。EDC模型假設(shè)湍流渦旋迅速耗散化學(xué)反應(yīng)，而PDF模型則考慮了化學(xué)反應(yīng)和湍流的統(tǒng)計分布。4.3多相燃燒模型多相燃燒模型適用于涉及固體、液體和氣體相的燃燒過程，如煤燃燒或噴霧燃燒。4.3.1相間相互作用在CFD++中，相間相互作用通過顆粒軌跡模型和液滴破碎模型來描述。顆粒軌跡模型考慮了固體顆粒在流場中的運動，而液滴破碎模型則描述了液滴在湍流中的破碎和蒸發(fā)。4.3.2燃燒區(qū)域劃分多相燃燒模型需要對燃燒區(qū)域進行劃分，以分別處理不同相的燃燒過程。這通常通過定義燃燒區(qū)域的邊界條件和相的初始條件來實現(xiàn)。4.3.3示例：設(shè)置多相燃燒模型以下是一個在CFD++中設(shè)置多相燃燒模型的示例代碼：[PHASES]

-GAS

-LIQUID

-SOLID

[PHASE_INTERACTIONS]

GAS-LIQUID:

-PARTICLE_TRAJECTORY_MODEL

-DROPLET_BREAKUP_MODEL

GAS-SOLID:

-PARTICLE_TRAJECTORY_MODEL

[COMBUSTION_REGIONS]

-GAS_REGION

BOUNDARY_CONDITIONS:

-INLET

-WALL

-OUTLET

SPECIES:

-O2

-CO2

-H2O

-N2

-CO

-CH4

-C2H4

-C2H6

-H2

-NO

-NO2

-N2O

-OH

-H

-O

-N

-NH

-CH3

-CH2

-CH

-C

-HCO

-HO2

-CH2O

-CH3O

-CH3OH

-CH2OH

-CHOH

-CH2O2

-CH3O2

-CH3OOH

-CH3O2H

-CH3O2

-CH3OH2

-CH3OH2+

-CH3O

-CH3O+

-CH3+

-CH2+

-CH+

-C+

-H2O2

-H2O+

-H2+

-H+

-e

-Ar

-Ar+

-Ar++

-Ar+++

-Ar++++

-Ar+++++

-Ar++++++

-Ar+++++++

-Ar++++++++

-Ar+++++++++

-Ar++++++++++

-Ar+++++++++++

-Ar++++++++++++

-Ar+++++++++++++

-Ar++++++++++++++

-Ar+++++++++++++++

-Ar++++++++++++++++

-Ar+++++++++++++++++

-Ar++++++++++++++++++

-Ar+++++++++++++++++++

-Ar++++++++++++++++++++

-Ar+++++++++++++++++++++

-Ar++++++++++++++++++++++

-Ar+++++++++++++++++++++++

-Ar++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

-Ar+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++