CFD中文网

这是编译截图

dell@dell-optiplex-7070:~/OpenFOAM/dell-v2012/run/myinterIsoFoam/myoverInterDyMFoam/myoverInterIsoDyMFoam1$
wmake | tee log.wmake2

Making dependency list for source file myoverInterIsoDyMFoam1.C

g++ -std=c++11 -m64 -pthread -DOPENFOAM=2012 -DWM_DP
-DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -Wno-unknown-pragmas -O3 -DNoRepository -ftemplate-depth-100 -I../myoverVoF -I../myoverinterFoam -I../myoverInterIsoFoam -I
省略
-lm -o /home/dell/OpenFOAM/dell-v2012/platforms/linux64GccDPInt32Opt/bin/myoverInterIsoDyMFoam1

--------------------------------------------------problrm---------------------------------------------这是具体的报错信息

/usr/bin/ld: Make/linux64GccDPInt32Opt/myoverInterIsoDyMFoam1.o: in
function Foam::isoAdvection::type() const': myoverInterIsoDyMFoam1.C:(.text._ZNK4Foam12isoAdvection4typeEv[_ZNK4Foam12isoAdvection4typeEv]+0x7): undefined reference to Foam::isoAdvection::typeName' /usr/bin/ld:
Make/linux64GccDPInt32Opt/myoverInterIsoDyMFoam1.o: in function
main': myoverInterIsoDyMFoam1.C:(.text.startup+0x238d): undefined reference to Foam::isoAdvection::isoAdvection(Foam::GeometricField<double,
Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double,
Foam::fvsPatchField, Foam::surfaceMesh> const&,
Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,
Foam::volMesh> const&)'
---------------------------------------------omitted--------------------
/usr/bin/ld: myoverInterIsoDyMFoam1.C:(.text.ZN4Foam12isoAdvection6advectINS_9zeroFieldES2_EEvRKT_RKT0[ZN4Foam12isoAdvection6advectINS_9zeroFieldES2_EEvRKT_RKT0]+0x671): undefined reference to `Foam::isoAdvection::debug'
collect2: error: ld returned 1 exit status
make: *** [/home/dell/OpenFOAM/OpenFOAM-v2012/wmake/makefiles/general:150：/home/dell/OpenFOAM/dell-v2012/platforms/linux64GccDPInt32Opt/bin/myoverInterIsoDyMFoam1] 错误 1

---

想利用头文件#include “isoAdvection.H”替换MULES，并耦合到重叠网格中间区，扩展openfoam的求解器
isoAdvection的实现文件和MULES位置类似，源代码位于~/src/transportModels/geometricVoF/

请问这个报错可能是那个方向的问题?感谢！

@卡洛 在 使用olaFlow与ofv1906自带造波应用重叠网格模拟浮体运动 中说：

@yuno 应该不是反射的问题，我在出口添加消波了

您好，我最到遇到反射的问题，我可以请教一下出口添加消波是怎么处理的，谢谢！

@一朵苍凉 非常谢谢你！

@一朵苍凉 好的，我认同你的解释，

@一朵苍凉 惯性矩这个我再确认一下，谢谢！上图显示的不是浮体表面网格，是背景网格和重叠区域网格的叠加一起。 目前的centerofmass以及orientation，在我的理解中似乎是一样的物理意义，质心？

@一朵苍凉 dynamicMeshDict 这时我的6自由度文件

@一朵苍凉 这个算例我在商软starccm+里面已经算过，目前在openfoamv2012进行同样的计算。网格仿造DTChull算例，利用snappyHexMesh生成，我的求解文件输出是这样的，运动输出的6Dof结果感觉比较大，所以我认为可能是我的dynamicsMeshdict文件可能有问题，但不知道那儿有问题，您方便的话，再帮我看以下，指点一二，感谢。
Create time

Create mesh for time = 0

Selecting dynamicFvMesh dynamicOversetFvMesh
Selecting motion solver: sixDoFRigidBodyMotion
Applying solid body motion to entire mesh
Selecting sixDoFSolver Newmark
Translational constraint tensor (0 0 0 0 0 0 0 0 1)
Rotational constraint tensor (1 0 0 0 1 0 0 0 1)

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Using LTS
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Creating cellMask field to block out hole cells

Creating interpolatedCells field

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 9.39433050916e-07, No Iterations 467
Reading/calculating face velocity Uf

Starting time loop

Flow time scale min/max = 8.6866267787e-07, 0.0291922459613
Smoothed flow time scale min/max = 8.6866267787e-07, 0.0291922459613
Time = 0.001

PIMPLE: iteration 1
forces forces:
rho: rho
Not including porosity effects

Restraint translationDamper: force (-0 -0 -0)
Restraint rotationDamper: moment (-0 -0 -0)
6-DoF rigid body motion
Centre of rotation: (0 0 -0.461009082997)
Centre of mass: (0 0 -0.461009082997)
Orientation: (0.999999999122 -4.18404783667e-05 2.19816547713e-06 4.18405388615e-05 0.999999998746 -2.75277355352e-05 -2.19701370075e-06 2.75278274834e-05 0.999999999619)
Linear velocity: (0 0 0.0218340063689)
Angular velocity: (0.0550532334843 0.00439510067628 0.0836815815307)
cellVolumeWeight : detected 2 mesh regions

zone:0 nCells:1428000
zone:1 nCells:1108410

Marking patch-cells on zone 0
Marking patch-cells on zone 1

After patch analysis : nCells : 2536410
other : 2371976
patch : 141190
overset: 23244

Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight
meshToMesh: Using AABBTree method
Overlap volume: 2.68136988135
cellVolumeWeight : detected 3 mesh regions after overset

Front : fraction:1 size:23970
Overset analysis : nCells : 2536410
calculated : 2512836
interpolated : 23432
hole : 142

Execution time for mesh.update() = 56.91 s
smoothSolver: Solving for alpha.water, Initial residual = 1.32501360637e-07, Final residual = 5.40621940608e-09, No Iterations 1
Phase-1 volume fraction = 0.80497110631 Min(alpha.water) = 0 Max(alpha.water) = 1.00094392533
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.804971106201 Min(alpha.water) = -5.22912457509e-16 Max(alpha.water) = 1.00094392533
smoothSolver: Solving for alpha.water, Initial residual = 1.32287726888e-07, Final residual = 5.33117835751e-09, No Iterations 1
Phase-1 volume fraction = 0.80497110685 Min(alpha.water) = -3.28734455684e-21 Max(alpha.water) = 1.00148843719
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.804971106742 Min(alpha.water) = -3.32478347046e-14 Max(alpha.water) = 1.00148843719
DILUPBiCGStab: Solving for p_rgh, Initial residual = 1, Final residual = 0.00904252579238, No Iterations 4
time step continuity errors : sum local = 4.98432915358e-06, global = -4.35030739554e-09, cumulative = -4.35030739554e-09
DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.000104613428998, Final residual = 9.48129717208e-07, No Iterations 95
time step continuity errors : sum local = 5.17819342099e-08, global = 9.75203425207e-09, cumulative = 5.40172685653e-09
smoothSolver: Solving for epsilon, Initial residual = 0.00318321052028, Final residual = 5.00754349231e-05, No Iterations 2
bounding epsilon, min: -0.00891808144257 max: 106.248958364 average: 1.03357845974
smoothSolver: Solving for k, Initial residual = 0.999999999979, Final residual = 0.0989395295117, No Iterations 1
PIMPLE: iteration 2
smoothSolver: Solving for alpha.water, Initial residual = 1.5141031309e-05, Final residual = 5.25145730347e-09, No Iterations 2
Phase-1 volume fraction = 0.804973015496 Min(alpha.water) = -7.90828268587e-26 Max(alpha.water) = 1.00074698501
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.804972981578 Min(alpha.water) = -0.000375075266952 Max(alpha.water) = 1.00074698501
smoothSolver: Solving for alpha.water, Initial residual = 1.54146759159e-05, Final residual = 5.35120277128e-09, No Iterations 2
Phase-1 volume fraction = 0.804974867545 Min(alpha.water) = -4.8534415386e-10 Max(alpha.water) = 1.0014871078
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.80497478029 Min(alpha.water) = -0.000141219315799 Max(alpha.water) = 1.0014871078
DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.91238817862e-05, Final residual = 9.99981139229e-07, No Iterations 15
time step continuity errors : sum local = 1.97312896592e-06, global = 8.73490744761e-07, cumulative = 8.78892471617e-07
DILUPBiCGStab: Solving for p_rgh, Initial residual = 6.95825944208e-06, Final residual = 8.05766273076e-07, No Iterations 19
time step continuity errors : sum local = 1.37353052288e-06, global = -1.53311028217e-07, cumulative = 7.25581443401e-07
smoothSolver: Solving for epsilon, Initial residual = 0.0476145055724, Final residual = 0.00249278364183, No Iterations 6
bounding epsilon, min: -1.06527371456 max: 3100.56090246 average: 2.84781772564
smoothSolver: Solving for k, Initial residual = 0.79259262557, Final residual = 0.0311989522415, No Iterations 4
ExecutionTime = 112.02 s ClockTime = 112 s

Flow time scale min/max = 1.85015907037e-06, 1e+15
Smoothed flow time scale min/max = 1.85015907037e-06, 0.0296574788498
Time = 0.002

PIMPLE: iteration 1
forces forces:
rho: rho
Not including porosity effects

Restraint translationDamper: force (-0 -0 -187.685118747)
Restraint rotationDamper: moment (-637.846763149 -50.9216364354 -969.534803614)
6-DoF rigid body motion
Centre of rotation: (0 0 -0.459430003293)
Centre of mass: (0 0 -0.459430003293)
Orientation: (0.998240479028 -0.0584060732538 0.0102311601265 0.058508037056 0.998237177305 -0.00996731887657 -0.00963097244877 0.0105483862658 0.999897982755)
Linear velocity: (0 0 3.13632540091)
Angular velocity: (6.05745573944 18.6026037453 117.298982883)
cellVolumeWeight : detected 2 mesh regions

zone:0 nCells:1428000
zone:1 nCells:1108410

Marking patch-cells on zone 0
Marking patch-cells on zone 1

After patch analysis : nCells : 2536410
other : 2371976
patch : 141190
overset: 23244

Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight
meshToMesh: Using AABBTree method
Overlap volume: 2.67649933669
cellVolumeWeight : detected 4 mesh regions after overset

Front : fraction:1 size:24277
Overset analysis : nCells : 2536410
calculated : 2512742
interpolated : 23425
hole : 243

Execution time for mesh.update() = 54.91 s
[3] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
[3] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> 一堆错误提示省略

Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.

---

mpirun noticed that process rank 3 with PID 0 on node dell-optiplex-7070 exited on signal 8 (Floating point exception).

@一朵苍凉 不知道您能不能帮我看下我的算例，看看那儿可以优化以下？我的算例是重叠网格进行自由运动计算。运算几步以后就崩了，运动不对，但不知道那儿存在问题，可以的话，我上传详细的细节给您。

@迟铁 你好，我想找您请教一下6dof的设置问题，期待得到您的回复！感谢

@dy大世界 你好，我最近也在做类似的一些学习，可以一起交流吗？不是是否方便加你qq？我的QQ1332169833

@enthusiasm 在 如何设置初速度 中说：

本人目前在做一个物体斜射入水的问题，需要计算物体入水后的弹道和速度变化，使用的是重叠网格。
我的疑惑是如何定义入水时的初速度，众所周知，有些方法可以设置恒定的速度或角速度，可这与我的问题模型不符，因此诚恳请教各位如何设置初速度，这个初速度是不是应该在 6DOF 模型里设置呢

请问一下师兄怎么设置？我最近再算一个高速入水，之前用的商业软件。同样的困惑。如有解答十分感谢！