CN104036095A - 基于区域分解的耦合高精度复杂外形流场快速算法 - Google Patents
基于区域分解的耦合高精度复杂外形流场快速算法 Download PDFInfo
- Publication number
- CN104036095A CN104036095A CN201410300877.2A CN201410300877A CN104036095A CN 104036095 A CN104036095 A CN 104036095A CN 201410300877 A CN201410300877 A CN 201410300877A CN 104036095 A CN104036095 A CN 104036095A
- Authority
- CN
- China
- Prior art keywords
- weno
- coupled interface
- coupling
- numerical
- region
- 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
Landscapes
- Complex Calculations (AREA)
Abstract
本发明提供了一种基于区域分解的耦合高精度DG和WENO方法求解双曲守恒律方程和Euler方程组的快速计算方法,该方法首先对原问题进行区域分解,即在物理边界附近区域采用结构或非结构的间断有限元方法(DG),在其余规则区域使用结构网格下的有限体差分型加权本质无振荡(WENO)格式。在处理区域耦合界面过程中,有两种处理方法,一种是守恒的耦合处理方法;另一种是非守恒的耦合处理方法。在实际计算过程中,我们使用坏单元指示器来判断界面附近解是否充分光滑,若解充分光滑,则界面处采用非守恒的耦合方法,否则,采用守恒的耦合方法。
Description
技术领域:
本发明涉及计算流体力学数值方法领域,特别是涉及一种求解双曲守恒律方程的高精度数值方法。
背景技术:
飞行器三维复杂流动的数值模拟和相关多目标优化问题是目前计算流体力学中的前沿热点问题,同时也是一个面向工程实际需求的应用问题。然而在目前的计算机规模和求解能力的条件下,目前主流的数值方法在计算效率上还不能满足这一工程实际应用问题的需要,解决这个问题的关键之一是提高流场解算器的效率。
目前流行的高精度数值方法主要包括间断有限元方法(DG),高精度有限体积方法,如(k-exact)有限体积方法以及高精度有限差分型方法,如有限差分型的加权本质无振荡格式(WENO)。间断有限元方法具有高精度和易处理复杂边界这些优点,但其计算量大,计算效率低,不能满足工程实际需求;高精度有限体积方法有传统的有限体积法推广而来,具有处理复杂边界的能力,但其重构模板一般不是紧致的,这给该方法应用到实际三维问题上带来了一定的困难,另外此方法的计算量也比较大,计算效率较低;有限差分型方法具有高精度和计算量小,计算效率高的优点,但有限差分方法一般只能在结构网格上应用,很难处理复杂物形和边界。故针对三维飞行器复杂流动数值模拟和相关多目标优化这一问题,目前缺少一种高精度、高效且易处理复杂边界的数值方法。
发明内容:
为了克服上述现有技术的不足,本发明提出一种基于区域分解的耦合DG和WENO方法,即将间断有限元(DG)方法和有限差分型的加权本质无振荡(WENO)格式以区域分解的方式耦合在一起,在复杂物形边界附近区域使用结构或非结构网格下的间断有限元方法处理计算区域的边界,同时在流场其余规则区域使用有限差分型WENO方法以大规模的提高计算效率。多区域耦合DG和WENO方法具有高阶精度,易处理复杂边界且计算量小,计算效率高等优势。
本发明的技术方案是:
对求解的实际问题首先进行区域分解,将整体计算区域划分为物形边界附近区域和其余规则计算区域;对物形边界附近区域使用结构或非结构网格进行区域剖分,对其余规则计算区域采用结构网格进行区域剖分;
对物形边界附近区域使用结构或非结构网格下的间断有限元方法进行初始化,对其余规则计算区域采用结构网格下的有限差分型WENO进行初始化;
在耦合界面处分别构造两类数值通量,一类是耦合界面处DG数值通量,一类是耦合界面处有限差分型WENO数值通量;
使用坏单元指示子判断耦合界面两侧邻居单元是否为坏单元,若存在坏单元,说明解在界面附近可能存在间断,故在界面处使用是守恒的耦合方式,若不存在坏单元,说明解在界面附近充分光滑,故在界面处使用非守恒的耦合方式;
确定界面处耦合方式之后,可对计算区域内各个子区域完成相应的空间离散,得到半离散个是的方程(组),此半离散方程(组)可使用三阶TVD Runge-Kutta方法求解。
本发明的有益效果是:
本发明结合了目前主流高精度算法的优点,采用区域分解的方式,耦合两种不同种类的高精度算法,从而达到易于适应和处理各种复杂边界的同时,大幅度提高计算效率的目的。耦合的DG和WENO方法可以方便的在混合网格上使用,在处理复杂物形边界附近区域时使用非结构网格的DG方法,在处理远场规则计算区域时,使用结构网格的有限差分型WENO方法。相比于传统的DG方法,耦合方法可以幅度降低计算量,提高计算效率;相比于传统的有限差分型格式,耦合方法在处理复杂边界时更加灵活,从而满足工程实际的需要。耦合方法作为计算流体力学中求解双曲守恒律的一类数值方法,在工程实际应用中有广阔的前景和应用价值。
附图说明:
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面对实施例的附图做简要介绍。
图1:NACA0012翼型示意图
图2:耦合算法针对计算NACA0012机翼绕流问题的区域分解示意图
图3:耦合算法混合网格下耦合界面局部示意图
图4:构造耦合界面处WENO-FD数值通量示意图
图5:构造耦合界面处DG数值通量示意图
图6:耦合算法计算亚音速NACA0012机翼绕流计算结果密度等值线示意图
图7:耦合算法整体流程示意图
具体实施方式:
下面将结合本发明实施例中的附图,对本发明实施例中的方法进行清楚、完整的描述。显然,所描述的实例仅仅是本发明的一个应用实例。基于本发明中的实例,本领域技术人员在没有做出创造性工作前提下所获得的所有其他实例,都属于本发明保护的范围。
针对图一物形,确定计算区域。在此问题中,需要我们进行数值模拟求解的问题时亚音速NACA0012机翼绕流问题,流动初始条件为马赫数Ma=0.4,攻角AoA=5.0°。我们选定的计算区域为一规则的矩形区域[-15.0,15.0]×[-15.0,15.0];
针对求解问题物形特征,进行区域分解。对于本问题,由于处理物形的非规则区域集中于机翼附近区域,因此我们将整体计算区域分为两大区域,区域一为机翼附近的非结构网格区域(图二中红色区域),该区域范围为[-0.4,1.4]×[-0.4,0.4],该区域我们使用非结构网格上的DG方法计算;区域二为规则的远场结构网格区域(图二中绿色区域),该区域为规则区域,可以使用结构网格的有限差分型WENO-FD方法计算。
确定网格尺度,进行网格剖分。针对本问题的特点和NACA0012机翼的物形特征长度,我们的网格剖分参数如下:网格尺寸h=0.05,非结构区域三角形网格单元数量N1=824,结构区域四边形网格单元数量N2=57456,计算区域总体网格数量为N=58520。
初始化流场区域,构造耦合界面处数值通量。在耦合界面处我们需要构造两类数值通量,分别是耦合界面处WENO-FD的数值通量和DG的数值通量。
构造耦合界面处WENO-FD数值通量步骤如下:
寻找和获得耦合界面构造WENO-FD数值通量所需的虚拟节点的位置及节点上的函数值Uh,如图4所示,设耦合界面的位置为在构造单元II+1,J处的WENO-FD数值通量时,需要虚拟节点II,J,II-1,J,II-2,J,这三个节点上的函数值Uh由该虚拟节点对应的三角单元上的DG解函数多项式提供:
其中u(l)(t)为单元上DG自由度,v(l)(x,y)相应的基函数;
在由虚拟节点II,J,II-1,J,II-2,J和WENO-FD区域内节点II+1,J,II+2,J做成的重构模板上,使用WENO-FD数值通量构造方法构造单元II+1,J位于耦合界面处的WENO-FD数值通量
构造耦合界面处DG数值通量的步骤如下:
构造耦合界面处的DG数值通量需要提供界面左右两侧高斯积分节点上的和如图5所示,其中可以通过DG求解区域内相应单元上的解函数多项式获得,对于则需要使用重构或插值手段获得,这里我们使用基于WENO的直接插值方法获得位于耦合界面处的对于二维情况下,基于WENO的直接插值,我们逐维度进行插值的方法,即先沿y-轴方向进行WENO插值,得到x-轴方向WENO插值所需节点值,再进行x-轴方向的WENO插值,从而得到耦合界面处位于高斯积分节点上的
以耦合界面上点为例,给出沿x-轴方向进行WENO型插值的方法如下:
选取x-轴方向WENO插值的插值模板S={II-2,J,II-1,J,II,J,II+1,J,II+2,J},将此模板分割为三个小模板S1={II-2,J,II-1,J,II,J},S2={II-1,JII,J,II+1,J},S3={II,J,II+1,J,II+2,J};
在每个小模板上构造拉格朗日插值多项式Pl(x),l=1,2,3,对于每个模板Sl内的拉格朗日多项式,需满足Pl(xi,yj)=Ui,j,Ii,j∈Sl;
计算每个小模板上的拉格朗日多项式的线性权重dl,l=1,2,3和光滑因子βl,l=1,2,3,得到每个小模板对应的非线性权,对于本例三个小模板Sl的线性权分别为:
每个小模板上的光滑因子βl的计算方法如下:
其中N为小模板上拉格朗日插值多项式的次数。于是,可以得到每个小模板对应的非线性权如下:
其中ε=1.0e-6;
计算耦合界面处点的函数值,该值由上述模板中各个多项式的值与相应的非线性权重加权组合得到
在完成耦合界面处WENO-FD数值通量和DG数值通量的构造后,在耦合界面位置使用坏单元指示子,判断耦合界面两侧单元是否为坏单元:若耦合界面两侧单元存在坏单元,则耦合界面处使用守恒的耦合方式,否则耦合界面处使用非守恒的耦合方式;对于守恒的耦合方式,即在耦合界面处采用唯一的数值通量,这里我们可以选取WENO-FD的数值通量,或选取DG数值通量;对于非守恒耦合方式,耦合界面处不同的子区域选择使用不同的数值通量,例如对于DG子区域,耦合界面处则使用DG数值通量,对于WENO-FD子区域,耦合界面处使用WENO-FD数值通量。
最后,对于不同的子区域,使用相应的重构的数值通量完成空间离散,得到半离散形式的常微分方程(组),可以使用三阶TVD Runge-Kutta方法求解,本实例的计算结果如图6所示。
Claims (5)
1.基于区域分解的耦合高精度复杂外形流场快速算法,其特征在于所述方法具体包括以下步骤:
(1)对求解的实际问题首先进行区域分解,将整体计算区域划分为物形边界附近区域和其余规则计算区域;对物形边界附近区域使用结构或非结构网格进行区域剖分,对其余规则计算区域采用结构网格进行区域剖分;
(2)对物形边界附近子区域使用结构或非结构网格下的间断有限元(DG)方法进行初始化和空间离散,对其余规则子区域采用结构网格下的有限差分型加权本质无振荡(WENO)格式进行初始化和空间离散;
(3)在耦合界面处的单元边界上分别构造两种数值通量,一种是耦合界面处单元边界上的DG数值通量,一类是耦合界面处单元边界上的有限差分型WENO数值通量;
(4)使用坏单元指示子判断耦合界面两侧邻居单元是否为坏单元:若存在坏单元,说明解在界面附近可能存在间断,此时在界面处使用是守恒的耦合方式;若不存在坏单元,说明解在界面附近充分光滑,此时在界面处使用非守恒的耦合方式;
(5)确定界面处耦合方式之后,可对计算区域内各个子区域完成相应的空间离散,得到半离散个是的方程(组)。
2.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法,其特征在于:在步骤(3)中构造有限差分型WENO数值通量过程中使用DG计算区域提供的虚拟节点,通过相应节点所属单元上DG解函数多项式可以得到该虚拟节点的值,然后使用有限差分型WENO格式方法构造耦合界面处单元的WENO数值通量。
3.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法,其特征在于:在步骤(3)中构造DG数值通量过程中需要使用基于WENO思想的拉格朗日插值得到计算DG数值通量所需的耦合界面处高斯积分节点的解函数近似值,从而构造耦合界面处单元DG的数值通量。
4.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法,其特征在于:在步骤(4)中守恒的耦合方式是在耦合界面处使用唯一的数值通量对耦合界面两侧的不同子区域进行空间离散,此数值通量可以为步骤(3)中在耦合界面处构造WENO数值通量,也可以为步骤(3)中在耦合界面处构造的DG数值通量。
5.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法,其特征在于:在步骤(4)中非守恒的耦合方式是在耦合界面两侧针对不同的子区域使用不同数值通量对耦合界面两侧的子区域进行空间离散,此时对有限差分型WENO方法处理的子区域使用步骤(3)中在耦合界面处构造的WENO数值通量,对DG方法处理的子区域使用步骤(3)中在耦合界面处构造的DG数值通量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300877.2A CN104036095B (zh) | 2014-06-27 | 2014-06-27 | 基于区域分解的耦合高精度复杂外形流场快速算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410300877.2A CN104036095B (zh) | 2014-06-27 | 2014-06-27 | 基于区域分解的耦合高精度复杂外形流场快速算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104036095A true CN104036095A (zh) | 2014-09-10 |
CN104036095B CN104036095B (zh) | 2018-01-02 |
Family
ID=51466866
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410300877.2A Active CN104036095B (zh) | 2014-06-27 | 2014-06-27 | 基于区域分解的耦合高精度复杂外形流场快速算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036095B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105975645A (zh) * | 2016-02-26 | 2016-09-28 | 西北工业大学 | 一种基于多步的含激波区域飞行器流场快速计算方法 |
CN106295088A (zh) * | 2016-11-03 | 2017-01-04 | 上海卫星工程研究所 | 航天器空间机热耦合载荷高精度自适应传递的通用方法 |
CN107422387A (zh) * | 2017-06-27 | 2017-12-01 | 吉林大学 | 一种虚拟时域有限差分的瞬变电磁发射源加载方法 |
CN107967382A (zh) * | 2017-11-20 | 2018-04-27 | 中国运载火箭技术研究院 | 一种吸气式飞行器耦合影响程度评价方法 |
CN108197367A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
CN108229083A (zh) * | 2018-04-11 | 2018-06-29 | 南京航空航天大学 | 一种基于改进的有限差分格式的水流数值模拟方法 |
CN111159853A (zh) * | 2019-12-10 | 2020-05-15 | 北京航空航天大学 | 一种针对高雷诺数粘性流动问题的边界层高精度处理方法 |
CN111563314A (zh) * | 2020-03-23 | 2020-08-21 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
CN114001654A (zh) * | 2021-11-01 | 2022-02-01 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
CN114254572A (zh) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及*** |
CN116702571A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070011644A1 (en) * | 2003-04-06 | 2007-01-11 | Daniel Abrams | Optimized photomasks for photolithography |
CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
-
2014
- 2014-06-27 CN CN201410300877.2A patent/CN104036095B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070011644A1 (en) * | 2003-04-06 | 2007-01-11 | Daniel Abrams | Optimized photomasks for photolithography |
CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
Non-Patent Citations (8)
Title |
---|
HONG LUO等: "A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids", 《45TH AIAA AEROSPACE SCIENCES MEETING AND EXHIBIT》 * |
JIAN CHENG等: "Multidomain hybrid RKDG and weno methods for hyperbolic conservation laws", 《SIAM J. SCI. COMPUT.》 * |
ZHEN-HUA JIANG 等: "Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids", 《ACTA MECH. SIN.》 * |
卢耀文: "求解双曲守恒律的非守恒型RKDG和WENO耦合方法", 《科技风》 * |
张来平等: "基于静动态混合重构的DG/FV混合格式", 《力学学报》 * |
朱洪强: "简化的WENO重构在h自适应RKDG方法中的应用", 《科学技术与工程》 * |
王坤等: "间断有限元方法处理边界的WEND-FD格式", 《第十一六届全国流体力学数值方法研讨会2013论文集》 * |
程剑等: "基于间断有限元方法与加权本质非振荡格式的多区域藕合算法", 《第十一六届全国流体力学数值方法研讨会2013论文集》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105975645A (zh) * | 2016-02-26 | 2016-09-28 | 西北工业大学 | 一种基于多步的含激波区域飞行器流场快速计算方法 |
CN105975645B (zh) * | 2016-02-26 | 2019-01-04 | 西北工业大学 | 一种基于多步的含激波区域飞行器流场快速计算方法 |
CN106295088A (zh) * | 2016-11-03 | 2017-01-04 | 上海卫星工程研究所 | 航天器空间机热耦合载荷高精度自适应传递的通用方法 |
CN106295088B (zh) * | 2016-11-03 | 2020-03-10 | 上海卫星工程研究所 | 航天器空间机热耦合载荷高精度自适应传递的通用方法 |
CN107422387A (zh) * | 2017-06-27 | 2017-12-01 | 吉林大学 | 一种虚拟时域有限差分的瞬变电磁发射源加载方法 |
CN107967382B (zh) * | 2017-11-20 | 2020-12-18 | 中国运载火箭技术研究院 | 一种吸气式飞行器耦合影响程度评价方法 |
CN107967382A (zh) * | 2017-11-20 | 2018-04-27 | 中国运载火箭技术研究院 | 一种吸气式飞行器耦合影响程度评价方法 |
CN108197367A (zh) * | 2017-12-27 | 2018-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种基于流场通量阶跃的高精度间断Galerkin人工粘性激波捕捉方法 |
CN108197367B (zh) * | 2017-12-27 | 2021-07-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高精度间断Galerkin人工粘性激波捕捉方法 |
CN108229083A (zh) * | 2018-04-11 | 2018-06-29 | 南京航空航天大学 | 一种基于改进的有限差分格式的水流数值模拟方法 |
CN111159853A (zh) * | 2019-12-10 | 2020-05-15 | 北京航空航天大学 | 一种针对高雷诺数粘性流动问题的边界层高精度处理方法 |
CN111159853B (zh) * | 2019-12-10 | 2022-07-15 | 北京航空航天大学 | 一种针对高雷诺数粘性流动问题的边界层高精度处理方法 |
CN111563314A (zh) * | 2020-03-23 | 2020-08-21 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
CN111563314B (zh) * | 2020-03-23 | 2023-09-12 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
CN114001654A (zh) * | 2021-11-01 | 2022-02-01 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
CN114001654B (zh) * | 2021-11-01 | 2024-03-26 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
CN114254572A (zh) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及*** |
CN114254572B (zh) * | 2021-12-16 | 2024-01-02 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及*** |
CN116702571A (zh) * | 2023-08-07 | 2023-09-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
CN116702571B (zh) * | 2023-08-07 | 2023-10-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104036095B (zh) | 2018-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN103150424B (zh) | 一种获取反应堆堆芯三维中子通量密度精细分布的方法 | |
CN106446432B (zh) | 一种求解材料大变形的最优输运无网格方法 | |
CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
CN103778298B (zh) | 改进的模拟多孔介质中二维水流运动的多尺度有限元方法 | |
CN110188423A (zh) | 一种基于有限元网格划分的线性工程结构快速bim建模方法 | |
CN109948214B (zh) | 城市多尺度风环境数值模拟方法 | |
CN103838913A (zh) | 曲线箱梁弯桥的有限单元法 | |
CN114943167B (zh) | 一种结构网格壁面距离的计算方法、***、介质和设备 | |
CN103914879A (zh) | 一种在抛物线方程中由三角面元数据生成立方网格数据的方法 | |
CN110765695A (zh) | 一种基于高阶有限元法获取混凝土重力坝裂纹扩展路径的模拟计算方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及***及设备 | |
CN105302964A (zh) | 一种用于芯片结构的热分布分析方法 | |
CN105404751B (zh) | 基于热-力-电磁场网络统一的实现方法 | |
CN104615812B (zh) | 一种高超声速返回舱外形气动热数值计算网格生成方法 | |
CN105677995A (zh) | 一种基于全网格配点理论的模糊稳态热传导问题数值求解方法 | |
CN112100890B (zh) | 一种压裂水平缝起裂扩展全三维数学模型的计算方法 | |
CN107808021B (zh) | 基于cfd的流体设备的阻力计算方法 | |
CN111797447B (zh) | 一种复杂地基混凝土坝孔道配筋方法 | |
CN111339688A (zh) | 基于大数据并行算法求解火箭仿真模型时域方程的方法 | |
CN102930143A (zh) | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 | |
CN111488704A (zh) | 一种排管敷设电缆外部热阻计算方法及*** | |
CN104809297A (zh) | 一种磁场-结构场耦合计算中的异型网格间的电磁力密度传递方法 | |
CN102073769B (zh) | 一种快速预测铸件热节的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210618 Address after: Room 501, 5 / F, building 4, 1418-66 Moganshan Road, Hangzhou, Zhejiang 310000 Patentee after: Sanduo (Hangzhou) Technology Co.,Ltd. Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road Patentee before: BEIHANG University |