CN104036095B - 基于区域分解的耦合高精度复杂外形流场快速算法 - Google Patents
基于区域分解的耦合高精度复杂外形流场快速算法 Download PDFInfo
- Publication number
- CN104036095B CN104036095B CN201410300877.2A CN201410300877A CN104036095B CN 104036095 B CN104036095 B CN 104036095B CN 201410300877 A CN201410300877 A CN 201410300877A CN 104036095 B CN104036095 B CN 104036095B
- Authority
- CN
- China
- Prior art keywords
- interface
- coupled
- region
- numerical
- conservation
- 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.)
- Active
Links
- 238000010168 coupling process Methods 0.000 title claims abstract description 28
- 230000008878 coupling Effects 0.000 title claims abstract description 24
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 title claims description 12
- 238000000034 method Methods 0.000 claims abstract description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000004907 flux Effects 0.000 claims description 40
- 238000013316 zoning Methods 0.000 claims description 13
- 230000007480 spreading Effects 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000004804 winding Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 230000001788 irregular Effects 0.000 claims description 2
- 239000012141 concentrate Substances 0.000 claims 1
- 230000003534 oscillatory effect Effects 0.000 claims 1
- 238000003672 processing method Methods 0.000 abstract 3
- 238000010586 diagram Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
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)对求解的作为实际问题的亚音速NACA0012机翼绕流问题首先进行区域分解，将整体计算区域划分为物形边界附近区域和其余规则计算区域；对物形边界附近区域使用结构或非结构网格进行区域剖分，对其余规则计算区域采用结构网格进行区域剖分；由于处理物形的非规则区域集中于机翼附近区域，因此将整体计算区域分为两大区域，区域一为机翼附近的非结构网格区域，该区域范围为[-0.4,1.4]×[-0.4,0.4]，该区域使用非结构网格上的间断有限元方法计算；区域二为规则的远场结构网格区域，该区域为规则区域，使用结构网格下的有限差分型加权本质无振荡格式进行初始化和空间离散；
(2)对物形边界附近子区域使用非结构网格下的间断有限元方法进行初始化和空间离散，对其余规则子区域采用结构网格下的有限差分型加权本质无振荡格式进行初始化和空间离散；
(3)在耦合界面处的单元边界上分别构造两种数值通量，一种是耦合界面处单元边界上的间断有限元数值通量，一类是耦合界面处单元边界上的有限差分型加权本质无振荡数值通量；
(4)使用坏单元指示子判断耦合界面两侧邻居单元是否为坏单元；若存在坏单元，说明解在界面附近可能存在间断，此时在界面处使用是守恒的耦合方式；若不存在坏单元，说明解在界面附近充分光滑，此时在界面处使用非守恒的耦合方式；
(5)确定界面处耦合方式之后，对计算区域内各个子区域完成相应的空间离散，得到半离散形式的方程。
2.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法，其特征在于：在步骤(3)中构造有限差分型加权本质无振荡数值通量过程中使用间断有限元计算区域提供的虚拟节点，通过相应节点所属单元上间断有限元解函数多项式得到该虚拟节点的值，然后使用有限差分型加权本质无振荡格式方法构造耦合界面处单元的加权本质无振荡数值通量。
3.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法，其特征在于：在步骤(3)中构造间断有限元数值通量过程中需要使用基于加权本质无振荡思想的拉格朗日插值得到计算间断有限元数值通量所需的耦合界面处高斯积分节点的解函数近似值，从而构造耦合界面处单元间断有限元的数值通量。
4.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法，其特征在于：在步骤(4)中守恒的耦合方式是在耦合界面处使用唯一的数值通量对耦合界面两侧的不同子区域进行空间离散，此数值通量为步骤(3)中在耦合界面处构造加权本质无振荡数值通量，或为步骤(3)中在耦合界面处构造的间断有限元数值通量。
5.根据权利要求1所述的基于区域分解的耦合高精度复杂外形流场快速算法，其特征在于：在步骤(4)中非守恒的耦合方式是在耦合界面两侧针对不同的子区域使用不同数值通量对耦合界面两侧的子区域进行空间离散，此时对有限差分型加权本质无振荡方法处理的子区域使用步骤(3)中在耦合界面处构造的加权本质无振荡数值通量，对间断有限元方法处理的子区域使用步骤(3)中在耦合界面处构造的间断有限元数值通量。
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 CN104036095A (zh) | 2014-09-10 |
| CN104036095B true 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) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105975645B (zh) * | 2016-02-26 | 2019-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 | 中国运载火箭技术研究院 | 一种吸气式飞行器耦合影响程度评价方法 |
| CN108197367B (zh) * | 2017-12-27 | 2021-07-27 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种高精度间断Galerkin人工粘性激波捕捉方法 |
| CN108229083A (zh) * | 2018-04-11 | 2018-06-29 | 南京航空航天大学 | 一种基于改进的有限差分格式的水流数值模拟方法 |
| CN111159853B (zh) * | 2019-12-10 | 2022-07-15 | 北京航空航天大学 | 一种针对高雷诺数粘性流动问题的边界层高精度处理方法 |
| CN111563314B (zh) * | 2020-03-23 | 2023-09-12 | 空气动力学国家重点实验室 | 一种七点weno格式的构造方法 |
| CN114001654B (zh) * | 2021-11-01 | 2024-03-26 | 北京卫星制造厂有限公司 | 工件端面位姿评价方法 |
| CN114254572B (zh) * | 2021-12-16 | 2024-01-02 | 西北工业大学太仓长三角研究院 | 考虑污染物沉积的航发压气机流场性能预测方法及系统 |
| CN116702571B (zh) * | 2023-08-07 | 2023-10-20 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于多重光滑度量因子的数值模拟方法及装置 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7480889B2 (en) * | 2003-04-06 | 2009-01-20 | Luminescent Technologies, Inc. | Optimized photomasks for photolithography |
-
2014
- 2014-06-27 CN CN201410300877.2A patent/CN104036095B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102830121A (zh) * | 2012-08-17 | 2012-12-19 | 浙江工业大学 | 一种软性磨粒流磨粒群实时检测方法 |
Non-Patent Citations (8)
| Title |
|---|
| A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids;Hong Luo等;《45th AIAA Aerospace Sciences Meeting and Exhibit》;20070510;全文 * |
| Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids;Zhen-Hua Jiang 等;《Acta Mech. Sin.》;20121231;全文 * |
| Multidomain hybrid RKDG and weno methods for hyperbolic conservation laws;JIAN CHENG等;《SIAM J. SCI. COMPUT.》;20130411;第35卷(第2期);第1-6节 * |
| 基于间断有限元方法与加权本质非振荡格式的多区域藕合算法;程剑等;《第十一六届全国流体力学数值方法研讨会2013论文集》;20130823;全文 * |
| 基于静动态混合重构的DG/FV混合格式;张来平等;《力学学报》;20101130;第42卷(第6期);全文 * |
| 求解双曲守恒律的非守恒型RKDG和WENO耦合方法;卢耀文;《科技风》;20110131;全文 * |
| 简化的WENO重构在h自适应RKDG方法中的应用;朱洪强;《科学技术与工程》;20130831;第13卷(第22期);全文 * |
| 间断有限元方法处理边界的WEND-FD格式;王坤等;《第十一六届全国流体力学数值方法研讨会2013论文集》;20130823;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104036095A (zh) | 2014-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104036095B (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
| CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
| CN110069800B (zh) | 具有光滑边界表达的三维结构拓扑优化设计方法及设备 | |
| CN106446432B (zh) | 一种求解材料大变形的最优输运无网格方法 | |
| CN103604729B (zh) | 一种颗粒随机分布的复合材料宏观有效性质的预测方法 | |
| CN104268322B (zh) | Weno差分方法的一种边界处理技术 | |
| CN103969627B (zh) | 基于fdtd的探地雷达大规模三维正演模拟方法 | |
| CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
| Han et al. | An adaptive GRP scheme for compressible fluid flows | |
| CN106682262A (zh) | 一种获取飞行器流场的数值模拟方法 | |
| CN106557638A (zh) | 焊接混合双向过渡单元网格模型的建立方法 | |
| CN109918821A (zh) | 一种迎风守恒型河道漫溢出流数值模拟方法 | |
| CN103838913A (zh) | 曲线箱梁弯桥的有限单元法 | |
| CN111125963A (zh) | 基于拉格朗日积分点有限元的数值仿真系统及方法 | |
| CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
| CN102646288B (zh) | 一种间隔基本均匀的三维流线放置方法 | |
| CN102063577B (zh) | 一种基于等势面的地下水渗流量计算方法 | |
| CN106844963A (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
| Zhang et al. | A simple technique to improve computational efficiency of meshless methods | |
| CN104809297A (zh) | 一种磁场-结构场耦合计算中的异型网格间的电磁力密度传递方法 | |
| CN104535040A (zh) | 用于叶片的有限元单元划分方法和叶片的检测方法 | |
| CN111488704A (zh) | 一种排管敷设电缆外部热阻计算方法及系统 | |
| CN104992046A (zh) | 流体力学计算系统及方法 | |
| Cai et al. | A robust algorithm for the generation of integration cells in Numerical Manifold Method | |
| CN107609274B (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 |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240819 Address after: 200000 Room 1903, 19th Floor, No. 17, Lane 166, Guchuan Road, Putuo District, Shanghai Patentee after: Haifang (Shanghai) Technology Co.,Ltd. Country or region after: China Address before: Room 501, 5 / F, building 4, 1418-66 Moganshan Road, Hangzhou, Zhejiang 310000 Patentee before: Sanduo (Hangzhou) Technology Co.,Ltd. Country or region before: China |