|
[Sponsors] |
November 12, 2012, 16:07 |
Initial conditions for uniform flow
|
#1 |
New Member
Join Date: Oct 2012
Posts: 14
Rep Power: 14 |
Hello,
I'm new in CFD and would like to know how to write the initial conditions for an incompressible fluid between two plane surfaces. I used channelFlow for the mesh, but have to change the 0 folder. My initial conditions are: p_inlet = p_outlet V_initial = 5 m/s As you imagine, the fluid will slow down untill it stops since there is no pressure gradient. My initial pressure conditions are: Code:
dimensions [0 2 -2 0 0 0 0]; internalField uniform 1; boundaryField { inlet { type fixedValue; value uniform 1; } outlet { type fixedValue; value uniform 1; } upAndDown { type zeroGradient; } frontAndBack { type empty; } } Code:
dimensions [0 1 -1 0 0 0 0]; internalField uniform (5 0 0); boundaryField { inlet { type cyclic; } outlet { type cyclic; } TopAndBottom { type fixedValue; value uniform (0 0 0); } frontAndBack { type empty; } } Thanks for your help, Andreas Last edited by andreas; November 12, 2012 at 19:39. |
|
November 15, 2012, 07:22 |
|
#2 |
Member
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 14 |
Maybe it's because you have called the same patches with 2 different names in the 2 files?
upAndDown & TopAndBottom |
|
November 15, 2012, 14:05 |
|
#3 |
New Member
Join Date: Oct 2012
Posts: 14
Rep Power: 14 |
It was not the problem, paraFoam still doesn't work.
And it still closes automaticaly when the "apply" button is pushed. This is the response on the terminal when blockMesh is generated: Code:
ubuntu@ubuntu-VirtualBox:~/OpenFOAM/Simu/U0$ blockMesh /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.1.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.1.1-221db2718bbb Exec : blockMesh Date : Nov 15 2012 Time : 18:51:18 Host : "ubuntu-VirtualBox" PID : 2587 Case : /home/ubuntu/OpenFOAM/Simu/U0 nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Creating block mesh from "/home/ubuntu/OpenFOAM/Simu/U0/constant/polyMesh/blockMeshDict" Creating curved edges Creating topology blocks Creating topology patches Creating block mesh topology Check topology Basic statistics Number of internal faces : 0 Number of boundary faces : 6 Number of defined boundary faces : 6 Number of undefined boundary faces : 0 Checking patch -> block consistency Creating block offsets Creating merge list . Creating polyMesh from blockMesh Creating patches Creating cells Creating points with scale 0.1 Writing polyMesh ---------------- Mesh Information ---------------- boundingBox: (0 0 0) (0.1 0.05 0.01) nPoints: 1722 nCells: 800 nFaces: 3260 nInternalFaces: 1540 ---------------- Patches ---------------- patch 0 (start: 1540 size: 20) name: inlet patch 1 (start: 1560 size: 20) name: outlet patch 2 (start: 1580 size: 80) name: upAndDown patch 3 (start: 1660 size: 1600) name: frontAndBack End Code:
ubuntu@ubuntu-VirtualBox:~/OpenFOAM/Simu/U0$ paraFoam created temporary 'U0.OpenFOAM' OpenGL Warning: Failed to connect to host. Make sure 3D acceleration is enabled for this VM. --> FOAM FATAL ERROR: Attempt to cast type wall to type lduInterface From function refCast<To>(From&) in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::error::abort() addr2line failed #2 Foam::lduInterface const& Foam::refCast<Foam::lduInterface const, Foam::fvPatch const>(Foam::fvPatch const&) addr2line failed #3 Foam::cyclicFvPatchField<Foam::Vector<double> >::cyclicFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) addr2line failed #4 Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::cyclicFvPatchField<Foam::Vector<double> > >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) addr2line failed #5 Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) addr2line failed #6 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricBoundaryField(Foam::fvBoundaryMesh const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libvtkPV3Foam.so" #7 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::dictionary const&) addr2line failed #8 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) addr2line failed #9 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) addr2line failed #10 void Foam::vtkPV3Foam::convertVolFields<Foam::Vector<double> >(Foam::fvMesh const&, Foam::PtrList<Foam::PrimitivePatchInterpolation<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > > > const&, Foam::IOobjectList const&, bool, vtkMultiBlockDataSet*) addr2line failed #11 Foam::vtkPV3Foam::convertVolFields(vtkMultiBlockDataSet*) addr2line failed #12 Foam::vtkPV3Foam::Update(vtkMultiBlockDataSet*, vtkMultiBlockDataSet*) addr2line failed #13 vtkPV3FoamReader::RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #14 vtkMultiBlockDataSetAlgorithm::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #15 vtkExecutive::CallAlgorithm(vtkInformation*, int, vtkInformationVector**, vtkInformationVector*) addr2line failed #16 vtkDemandDrivenPipeline::ExecuteData(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #17 vtkCompositeDataPipeline::ExecuteData(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #18 vtkDemandDrivenPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #19 vtkStreamingDemandDrivenPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #20 vtkCompositeDataPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #21 vtkCompositeDataPipeline::ForwardUpstream(vtkInformation*) addr2line failed #22 vtkDemandDrivenPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #23 vtkStreamingDemandDrivenPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #24 vtkCompositeDataPipeline::ProcessRequest(vtkInformation*, vtkInformationVector**, vtkInformationVector*) addr2line failed #25 vtkDemandDrivenPipeline::UpdateData(int) addr2line failed #26 vtkStreamingDemandDrivenPipeline::Update(int) addr2line failed #27 vtkSISourceProxy::UpdatePipeline(int, double, bool) addr2line failed #28 vtkSISourceProxyCommand(vtkClientServerInterpreter*, vtkObjectBase*, char const*, vtkClientServerStream const&, vtkClientServerStream&) addr2line failed #29 vtkClientServerInterpreter::ProcessCommandInvoke(vtkClientServerStream const&, int) addr2line failed #30 vtkClientServerInterpreter::ProcessOneMessage(vtkClientServerStream const&, int) addr2line failed #31 vtkClientServerInterpreter::ProcessStream(vtkClientServerStream const&) addr2line failed #32 vtkPVSessionCore::ExecuteStreamInternal(vtkClientServerStream const&, bool) addr2line failed #33 vtkPVSessionCore::ExecuteStream(unsigned int, vtkClientServerStream const&, bool) addr2line failed #34 vtkPVSessionBase::ExecuteStream(unsigned int, vtkClientServerStream const&, bool) addr2line failed #35 vtkSMProxy::ExecuteStream(vtkClientServerStream const&, bool, unsigned int) addr2line failed #36 vtkSMOutputPort::UpdatePipelineInternal(double, bool) addr2line failed #37 vtkSMOutputPort::UpdatePipeline(double) addr2line failed #38 vtkSMSourceProxy::UpdatePipeline(double) addr2line failed #39 pqPipelineSource::updatePipeline() addr2line failed #40 pqDisplayPolicy::getPreferredViewType(pqOutputPort*, bool) const addr2line failed #41 pqDisplayPolicy::getPreferredView(pqOutputPort*, pqView*) const addr2line failed #42 pqDisplayPolicy::createPreferredRepresentation(pqOutputPort*, pqView*, bool) const addr2line failed #43 pqObjectInspectorWidget::show(pqPipelineSource*) addr2line failed #44 pqObjectInspectorWidget::accept() addr2line failed #45 at moc_pqObjectInspectorWidget.cxx:0 #46 QMetaObject::activate(QObject*, QMetaObject const*, int, void**) addr2line failed #47 QAbstractButton::clicked(bool) addr2line failed #48 addr2line failed #49 addr2line failed #50 QAbstractButton::mouseReleaseEvent(QMouseEvent*) addr2line failed #51 QWidget::event(QEvent*) addr2line failed #52 QAbstractButton::event(QEvent*) addr2line failed #53 QPushButton::event(QEvent*) addr2line failed #54 QApplicationPrivate::notify_helper(QObject*, QEvent*) addr2line failed #55 QApplication::notify(QObject*, QEvent*) addr2line failed #56 QCoreApplication::notifyInternal(QObject*, QEvent*) addr2line failed #57 QApplicationPrivate::sendMouseEvent(QWidget*, QMouseEvent*, QWidget*, QWidget*, QWidget**, QPointer<QWidget>&, bool) addr2line failed #58 addr2line failed #59 QApplication::x11ProcessEvent(_XEvent*) addr2line failed #60 addr2line failed #61 g_main_context_dispatch addr2line failed #62 addr2line failed #63 g_main_context_iteration addr2line failed #64 QEventDispatcherGlib::processEvents(QFlags<QEventLoop::ProcessEventsFlag>) addr2line failed #65 addr2line failed #66 QEventLoop::processEvents(QFlags<QEventLoop::ProcessEventsFlag>) addr2line failed #67 QEventLoop::exec(QFlags<QEventLoop::ProcessEventsFlag>) addr2line failed #68 QCoreApplication::exec() addr2line failed #69 QApplication::exec() addr2line failed #70 main addr2line failed #71 __libc_start_main addr2line failed #72 addr2line failed Aborted (core dumped) |
|
November 15, 2012, 14:12 |
|
#4 | |
Member
Simone
Join Date: Sep 2012
Posts: 95
Rep Power: 14 |
Maybe there's something wrong in the patch wall
Quote:
|
||
November 16, 2012, 11:20 |
|
#5 |
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 23 |
Hi,
you defined inlet and outlet as cyclic (as in U), but you give them a fixedValue. you should use cyclic instead. mad |
|
November 16, 2012, 16:00 |
|
#6 | |
Senior Member
Lieven
Join Date: Dec 2011
Location: Leuven, Belgium
Posts: 299
Rep Power: 23 |
I'm somewhat confused why you would say
Quote:
So your flow will certainly not stop. The fact that this pressure gradient is imposed independently from p means however, that on average this field will be equal to 0. But if you want to have the total pressure in the channel, you need to account for the pressure gradient in a post-processing step. Regards, L |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Interfoam blows on parallel run | danvica | OpenFOAM Running, Solving & CFD | 16 | December 22, 2012 03:09 |
Orifice Plate with a fully developed flow - Problems with convergence | jonmec | OpenFOAM Running, Solving & CFD | 3 | July 28, 2011 06:24 |
MRFSimpleFOAM goes divergenced! | renyun0511 | OpenFOAM Running, Solving & CFD | 0 | November 19, 2009 03:11 |
Differences between serial and parallel runs | carsten | OpenFOAM Bugs | 11 | September 12, 2008 12:16 |
Convergence moving mesh | lr103476 | OpenFOAM Running, Solving & CFD | 30 | November 19, 2007 15:09 |