CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Boudary conditions in dsmcfoamplus problem

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By arte

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 13, 2024, 19:39
Exclamation Boudary conditions in dsmcfoamplus problem
  #1
New Member
 
Artem
Join Date: Nov 2024
Posts: 1
Rep Power: 0
arte is on a distinguished road
Hello everyone!
I decided to calculate the free-molecular flow around bodies of rotation. I am using the dsmcfoamplus module in OpenFoam 2.4.0-MNF. As an example I use the calculation of the Orion capsule (https://github.com/Eliasfarah0/dsmcF...dule/tree/main). For my example I set the flow velocity to 7000 m/s, the temperature to 196 K and the same gas parameters as in the Orion example. As a result of the calculation, I want to obtain the pressure and temperature distribution over the surface of the streamlined body. BUT! In the example and in the articles, the boundariesDict initially sets the boundary condition DiffuseWall and the temperature at the body wall, although this is something that I do not know initially and want to find. It seems that the zeroGradient boundary condition should be used, but OpenFoam gives an error if I use it. Please help me how to correctly set the boundary condition on the body wall.
The boudariesdict code is below, the model image is https://ibb.co/hHxTBn1.
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.4.0-MNF                             |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      boundariesDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// **************** //
// PATCH BOUNDARIES //
// **************** //
dsmcPatchBoundaries
(
    boundary
    {
        patchBoundaryProperties
        {
            patchName		inlet;
        }
        boundaryModel		dsmcDeletionPatch;
        dsmcDeletionPatchProperties
        {
            allSpecies		yes;
        }
    }

    boundary
    {
        patchBoundaryProperties
        {
            patchName		outlet;
        }
        boundaryModel		dsmcDeletionPatch;
        dsmcDeletionPatchProperties
        {
            allSpecies		yes;
        }
    }

    boundary
    {
        patchBoundaryProperties
        {
            patchName		top;
        }
        boundaryModel		dsmcDeletionPatch;
        dsmcDeletionPatchProperties
        {
            allSpecies		yes;
        }
    }

    boundary
    {
        patchBoundaryProperties
        {
            patchName		bottom;
        }
        boundaryModel		dsmcDeletionPatch;
        dsmcDeletionPatchProperties
        {
            allSpecies		yes;
        }
    }

    boundary
    {
        patchBoundaryProperties
        {
            patchName		ocm;
        }
        boundaryModel		dsmcDiffuseWallPatch;
        dsmcDiffuseWallPatchProperties
        {
            velocity		(0 0 0);
            temperature		196.605; // from Orion Aerodynamics for Hypersonic Free Molecular to Continuum Conditions
        }
    }
    
    
);



// ****************** //
// GENERAL BOUNDARIES //
// ****************** //
dsmcGeneralBoundaries
(
    boundary
    {
        generalBoundaryProperties
        {
            patchName			inlet;
        }
        boundaryModel   dsmcFreeStreamInflowPatch;
        dsmcFreeStreamInflowPatchProperties
        {
            typeIds			(N2 O2 O);
            velocity			(7000 0 0);	
            translationalTemperature	196.605;	 
            rotationalTemperature	196.605;	 
            vibrationalTemperature	196.605;	 
            electronicTemperature	0;		
            numberDensities
            {
                N2			8.892e18; 	
                O2			2.106e18; 	
		O			7.020e17; 	
            };
        }
    }

    boundary
    {
        generalBoundaryProperties
        {
            patchName			outlet;
        }
        boundaryModel   dsmcFreeStreamInflowPatch;
        dsmcFreeStreamInflowPatchProperties
        {
            typeIds			(N2 O2 O);
            velocity			(7000 0 0);	
            translationalTemperature	196.605;	 	
            rotationalTemperature	196.605;	 	
            vibrationalTemperature	196.605;	 	
            electronicTemperature	0;		
            numberDensities
            {
                N2			8.892e18; 	
                O2			2.106e18; 	
		O			7.020e17; 	
            };
        }
    }

    boundary
    {
        generalBoundaryProperties
        {
            patchName			top;
        }
        boundaryModel   dsmcFreeStreamInflowPatch;
        dsmcFreeStreamInflowPatchProperties
        {
            typeIds			(N2 O2 O);
            velocity			(7000 0 0);	
            translationalTemperature	196.605;	 	
            rotationalTemperature	196.605;	 	
            vibrationalTemperature	196.605;	 	
            electronicTemperature	0;		
            numberDensities
            {
                N2			8.892e18; 	
                O2			2.106e18; 	
		O			7.020e17; 	
            };
        }
    }

    boundary
    {
        generalBoundaryProperties
        {
            patchName			bottom;
        }
        boundaryModel   dsmcFreeStreamInflowPatch;
        dsmcFreeStreamInflowPatchProperties
        {
            typeIds			(N2 O2 O);
            velocity			(7000 0 0);
            translationalTemperature	196.605;	 
            rotationalTemperature	196.605;	 
            vibrationalTemperature	196.605;	 
            electronicTemperature	0;		
            numberDensities
            {
                N2			8.892e18; 
                O2			2.106e18;
		O			7.020e17;
            };
        }
    }
);



// ***************** //
// CYCLIC BOUNDARIES //
// ***************** //
dsmcCyclicBoundaries
(
	// EMPTY
);
// ************************************************************************* //
R2P3 likes this.
arte is offline   Reply With Quote

Reply

Tags
boundary condition, dsmcfoam+


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boundary conditions problem pj6620 FLOW-3D 9 June 4, 2014 06:15
Problem Boundary Conditions, forceCoeff valerio galluzzi OpenFOAM Pre-Processing 9 May 7, 2014 12:16
Problem in Cell Zone conditions for Species Transport model CFDP FLUENT 3 March 5, 2013 05:14
Problem with periodic boundary conditions jm_ngs FLUENT 7 April 27, 2012 03:32
Fluent boundary conditions problem bobo FLUENT 2 July 3, 2009 07:28


All times are GMT -4. The time now is 06:52.