HFSS混合算法求解器研发历程与展望---Ansys

1©2016ANSYS,Inc.August11,2016ANSYSUGM2016HFSS混合算法求解器研发历程与展望赵克钟博士ANSYS首席研发专家2©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegionHFSS算法体系与混合算法求解器3©2016ANSYS,Inc.August11,2016ANSYSUGM2016FE-PML混合FE-BI分离FE-BIs混合FE-BI和IERegionN=7.7M,137GB,136minN=777K,18GB,29minN=5.0M,78GB,101minN=127K,5GB,6min混合FE-BI和PORegionN=127K,1.6GB,1.5minHFSS求解器演化历程4©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegion5©2016ANSYS,Inc.August11,2016ANSYSUGM2016FE-BI:FiniteElement–BoundaryIntegral有限元-边界积分无最近距离限制◦推荐最近/10无反射边界条件◦能理想吸收入射场，与辐射角度无关任意形状的边界◦能包含分离区域◦边界与辐射体共形，减小了空气体积FE-BI计算消耗少◦空气盒子比ABC或PML小，能显著减小运算量*注:ABC（AbsorbingBoundaryCondition）即常说的Radiation边界IE（积分方程法）求解表面FEM（有限元法）求解空间区域设置方法:•选择Objects•AssignHybridFE-BI11©2016ANSYS,Inc.August11,2016ANSYSUGM2016FE-BI与ABC对比FE-BIABCABCFE-BI12©2016ANSYS,Inc.August11,2016ANSYSUGM2016FE-BI与ABC对比•FE-BI的内存占用仅有共形ABC的60%左右◦共形ABC边界并不理想，因为对边界表面的斜入射波吸收不好◦PML边界尽管吸收良好，但空气盒子体积大很多，增加了仿真运算量MethodRAMElapsedTimeABC15GB70min共形ABC11GB40minFE-BI7GB30min共形ABCFE-BI50cmABC13©2016ANSYS,Inc.August11,2016ANSYSUGM2016PMLFE-BI与PML对比•FEM-IE混合法可显著减少计算资源需求◦FEM无需计算极耗资源的罩内空气区域若使用ABC或PML则需求解完整的空气区域◦使用两个相互分离的FEM-IE区域与天线罩完美共形与喇叭天线(26GHz)共形两个区域间的耦合自动计算26GHz占用内存消耗时间PML259G(DDM)840minFE-BI64G205min19143λ³使用FEBI边界后获得4.1倍的提速并节省75%的内存FE-BIFE-BI表面2860λ³14©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegion15©2016ANSYS,Inc.August11,2016ANSYSUGM2016❷IE-Region❶FE-BI❶FE-BI❶FE-BIPMLSetup10625λ³HybridFE-BI37λ³18.4meter求解器未知量内存时间DirectSolver13,023,494187GB10:47:08DDM并行(23domains)13,023,494116GB01:49:44混合算法528,70918GB00:33:4724cores改进的IE-Region16©2016ANSYS,Inc.August11,2016ANSYSUGM2016与FEM对比，IE区域混合法可减少81%的内存消耗140GRAM44GRAMBlack–FEM-ABCRed–HybridIE-RegionsFEMwithABC231GRAM22463λ³199λ³HybridIE-Region44GRAM❶FE-BI❷IE-Region50λContinuousCurrentacrossFEMandIE-Region改进的IE-Region17©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegion18©2016ANSYS,Inc.August11,2016ANSYSUGM2016混合算法(6.2GRAM)DielectricCavity(介质腔体）❶介质IE-Region•金属区域可由IE直接求解表面电流•无需在金属周围设置辐射边界•均匀介质也可以可利用IE区域求解•介质表面由IE求解金属介质εr=4有限元求解(20GRAM)设置方法:•选择Objects•AssignHybridDielectricCavity相比于ABC边界，使用IE区域可减少70%的内存❷金属IE-Region19©2016ANSYS,Inc.August11,2016ANSYSUGM2016DielectricCavity案例1:办公室边界类型网格量未知数总内存(GB)耗时(小时)ABC412,6425,861,2834411.5FE-BI277,9723,594,473303BIstatistic:N=395,598RAM=7.5GBTime=9minFE-BI收敛情况20©2016ANSYS,Inc.August11,2016ANSYSUGM2016DielectricCavity案例2:人头部+手机求解方法内存时间(小时)仅FEM6.2GB1FE-BI3.0GB0.5人头部材料特性:εr=79,σ=0.47simems/m仅手机FEM:手机+人头FE-BI:手机+人头仅FEMFE-BI1.8GHz21©2016ANSYS,Inc.August11,2016ANSYSUGM20161.FEM:人头和手机周围2.内部BI:介质人体3.IERegion:金属车体求解方式总内存耗时FEM+DDM并行160GB8hours混合算法求解11GB2.7hoursFEM-IE混合算法求解完全FEM求解应用案例:天线布局在车内人体附近22©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegion26©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBIMoMPO反射面与真实馈源求解器未知量内存迭代步数仿真时间IERegion127,9648.2GB1600:06:18PORegion127,9643.1GB1500:01:29蓝色实线:E-fieldPO红色实线:H-fieldPO黑色虚线:IE27©2016ANSYS,Inc.August11,2016ANSYSUGM2016反射面与真实馈源PORegionIERegion28©2016ANSYS,Inc.August11,2016ANSYSUGM2016PORegionvs.IERegion:碟形反射面@50GHzPO/IERegionFEBISolverUnknownsMemoryIterationsRuntimeIERegion317,13417.3GB8.000:11:29PORegion317,1343.80GB9.500:01:34IERegionPORegion29©2016ANSYS,Inc.August11,2016ANSYSUGM2016双向耦合vs.单向链接馈源与反射面双向场耦合馈源与反射面单向场链接30©2016ANSYS,Inc.August11,2016ANSYSUGM2016求解器未知量内存迭代步数仿真时间IERegion1,367,400129GB10160minPORegion1,367,40014GB109.5minPORegionIERegionIE/PORegionFEBIFEBI+IE/PORegion案例1:58碟形反射面31©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBI+IE/PORegion案例2:双反射面系统FEBIPOPO求解器未知量内存迭代步数仿真时间IERegion369,15418.5GB800:12:52PORegion369,1545.3GB800:02:32PORegionIERegion32©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBI+IE/PORegion案例3:波束波导系统POFEBIPOPO求解器未知量内存迭代步数仿真时间IERegion1,177,97758.1GB7.374minPORegion1,177,97728.6GB8.726min33©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBI+IERegion+PORegion案例:FSS卡塞格伦天线IERegionPORegionFE-BI工作频段：K波段馈源天线：圆极化波纹喇叭馈源主反射面：直径2700mm，200波长次反射面：直径530mm，40波长，FSS边框宽仅0.3mm34©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBI+IERegion+PORegion案例:FSS卡塞格伦天线FSS网格剖分细节仿真时间3h37min迭代步数3passes总网格量62,404Tetra+2,868,344Tri内存消耗144GB35©2016ANSYS,Inc.August11,2016ANSYSUGM2016r=3.9,losstangent=0.025FEBI+DPO案例1:喇叭天线与介质罩求解器未知量内存迭代步数仿真时间FEBI1,745,57124.9GB1000:20:31DPO146,1181.1GB900:01:1136©2016ANSYS,Inc.August11,2016ANSYSUGM2016FEBI+DPO案例2:喇叭天线与介质罩@26GHz求解器未知量内存迭代步数仿真时间FEBI3,861,85545.4GB1600:39:29DPO573,6406.4GB700:03:52r=3.4,losstangent=0.01237©2016ANSYS,Inc.August11,2016ANSYSUGM2016基于区域分解的混合算法求解器FE-BIIERegionPORegionDielectricPODielectricIERegionDielectricCavityLinkedRegion40©2016ANSYS,Inc.August11,2016ANSYSUGM2016r=3.9,losstangent=0.025FEBI+DIE案例1:喇叭天线与介质罩求解器未知量内存仿真时间FEBI1,745,57124.9GB11minDIE(OLD)233,43321.5GB23minDIE(NEW)427,6219.8GB5min4倍提速