- •TABLE OF CONTENTS
- •Chapter 1 INTRODUCTION
- •The es-ice Environment
- •es-ice Meshing Capabilities
- •Tutorial Structure
- •Trimming Tutorial Overview
- •Required Files
- •Trimming Tutorial files
- •Automatic 2D Tutorial files
- •Wall Temperature Tutorial files
- •Mesh Replacement Tutorial files
- •Multiple Cylinder Tutorial files
- •Closed-Cycle Tutorial files
- •Sector Tutorial files
- •Two-Stroke Tutorial files
- •Mapping Tutorial files
- •ELSA Tutorial files
- •Chapter 2 SURFACE PREPARATION IN STAR-CCM+
- •Importing and Scaling the Geometry
- •Creating Features
- •Defining Surfaces
- •Remeshing and Exporting the Geometry
- •Chapter 3 GEOMETRY IMPORT AND VALVE WORK
- •Importing the Surfaces
- •Modelling the Valves
- •Saving the Model
- •Chapter 4 MESHING WITH THE TRIMMING METHOD
- •Modifying Special Cell Sets in the Geometry
- •Defining Flow Boundaries
- •Creating the 2D Base Template
- •Creating the 3D Template
- •Trimming the 3D Template to the Geometry
- •Improving cell connectivity
- •Assembling the Trimmed Template
- •Running Star Setup
- •Saving the Model
- •Chapter 5 CREATING AND CHECKING THE MESH
- •Chapter 6 STAR SET-UP in es-ice
- •Load Model
- •Analysis Set-up
- •Valve Lifts
- •Assembly
- •Combustion
- •Initialization
- •Cylinder
- •Port 1 and Port 2
- •Boundary Conditions
- •Cylinder
- •Port and Valve 1
- •Port and Valve 2
- •Global settings
- •Post Set-up
- •Cylinder
- •Port 1 and Port 2
- •Global settings
- •Time Step Control
- •Write Data
- •Saving the Model
- •Chapter 7 STAR SET-UP in pro-STAR
- •Using the es-ice Panel
- •Setting Solution and Output Controls
- •File Writing
- •Chapter 8 RUNNING THE STAR SOLVER
- •Running in Serial Mode
- •Running in Parallel Mode
- •Running in Parallel on Multiple Nodes
- •Running in Batch
- •Restarting the Analysis
- •Chapter 9 POST-PROCESSING: GENERAL TECHNIQUES
- •Creating Plots with the es-ice Graph Tool
- •Calculating Apparent Heat Release
- •Plotting an Indicator Diagram
- •Calculating Global Engine Quantities
- •Creating a Velocity Vector Display
- •Creating an Animation of Fuel Concentration
- •Creating an Animation of Temperature Isosurfaces
- •Chapter 10 USING THE AUTOMATIC 2D TEMPLATE
- •Importing the Geometry Surface
- •Defining Special Cell Sets in the Geometry
- •Modelling the Valves
- •Creating the Automatic 2D Template
- •Refining the 2D Template Around the Injector
- •Adding Features to the Automatic 2D Template
- •Using Detailed Automatic 2D Template Parameters
- •Saving the es-ice Model File
- •Chapter 11 MULTIPLE-CYCLE ANALYSIS
- •Setting Up Multiple Cycles in es-ice
- •Setting Up Multiple Cycles in pro-STAR
- •Chapter 12 HEAT TRANSFER ANALYSIS
- •Resuming the es-ice Model File
- •Mapping Wall Temperature
- •Exporting Wall Heat Transfer Data
- •Saving the es-ice Model File
- •Cycle-averaging Wall Heat Transfer Data
- •Post-processing Wall Heat Transfer Data in pro-STAR
- •Plotting average wall boundary temperatures
- •Plotting average heat transfer coefficients
- •Plotting average near-wall gas temperature at Y-plus=100
- •Mapping Heat Transfer Data to an Abaqus Model via STAR-CCM+
- •Chapter 13 MESH REPLACEMENT
- •Preparing the File Structure
- •Rebuilding the Dense Mesh
- •Creating Ahead Files for the Dense Mesh
- •Defining Mesh Replacements
- •Setting Up Mesh Replacement in pro-STAR
- •Setting up the coarse model
- •Setting up the dense model
- •Chapter 14 MULTIPLE CYLINDERS
- •Resuming the es-ice Model File
- •Making, Cutting and Assembling the Template
- •Setting Up Multiple Cylinders
- •Checking the Computational Mesh
- •STAR Set-Up in es-ice
- •Analysis set-up
- •Assembly
- •Combustion
- •Initialization
- •Boundary Conditions
- •Post Setup
- •Time Step Control
- •Write Data
- •Saving the es-ice Model File
- •Importing the Geometry
- •Generating the Closed-Cycle Polyhedral Mesh
- •Assigning shells to geometry cell sets
- •Specifying General, Events and Cylinder parameters
- •Creating a spray-optimised mesh zone
- •Importing a user intermediate surface
- •Checking the spray-optimised zone
- •Creating the closed-cycle polyhedral mesh
- •Running Star Setup
- •Creating and checking the computational mesh
- •Saving the Model File
- •Chapter 16 DIESEL ENGINE: SECTOR MODEL
- •Importing the Bowl Geometry
- •Defining the Bowl Shape
- •Defining the Fuel Injector
- •Creating the 2D Template
- •Creating the Sector Mesh
- •Creating and Checking the Mesh
- •Saving the Model
- •Chapter 17 DIESEL ENGINE: STAR SET-UP IN es-ice and pro-STAR
- •STAR Set-up in es-ice
- •Load model
- •Analysis setup
- •Assembly
- •Combustion
- •Initialization
- •Boundary conditions
- •Post setup
- •Time step control
- •Write data
- •Saving the Model File
- •STAR Set-up in pro-STAR
- •Using the es-ice Panel
- •Selecting Lagrangian and Liquid Film Modelling
- •Setting up the Fuel Injection Model
- •Setting up the Liquid Film Model
- •Setting up Analysis Controls
- •Writing the Geometry and Problem Files and Saving the Model
- •Chapter 18 DIESEL ENGINE: POST-PROCESSING
- •Creating a Scatter Plot
- •Creating a Spray Droplet Animation
- •Chapter 19 TWO-STROKE ENGINES
- •Importing the Geometry
- •Meshing with the Trimming Method
- •Assigning shells to geometry cell sets
- •Creating the 2D template
- •Creating the 3D template
- •Trimming the 3D template to the geometry
- •Assembling the trimmed template
- •Running Star Setup
- •Checking the mesh
- •STAR Set-up in es-ice
- •Analysis setup
- •Assembly
- •Combustion
- •Initialization
- •Boundary conditions
- •Post setup
- •Time step control
- •Write data
- •Saving the es-ice Model File
- •Chapter 20 MESHING WITH THE MAPPING METHOD
- •Creating the Stub Surface in the Geometry
- •Creating the 2D Base Template
- •Creating the 3D Template
- •General Notes About Edges and Splines
- •Creating Edges and Splines Near the Valve Seat
- •Creating the Remaining Edges and Splines
- •Creating Patches
- •The Mapping Process
- •Chapter 21 IMPROVING THE MAPPED MESH QUALITY
- •Creating Plastered Cells
- •Chapter 22 PISTON MODELING
- •Meshing the Piston with the Shape Piston Method
- •Chapter 23 ELSA SPRAY MODELLING
- •Importing the Bowl Geometry
- •Defining the Bowl Shape
- •Setting the Events and Cylinder Parameters
- •Creating the Spray Zone
- •Creating the Sector Mesh
- •STAR Set-up in es-ice
- •Load model
- •Analysis setup
- •Assembly
- •Combustion
- •Initialization
- •Boundary Conditions
- •Time step control
- •Write data
- •Saving the Model File
- •STAR Set-up in pro-STAR
- •Using the es-ice panel
- •Activating the Lagrangian model
- •Defining the ELSA scalars
- •Setting up the Lagrangian droplets
- •Defining boundary regions and boundary conditions
- •Setting up analysis controls
- •Adding extended data for the ELSA model
- •Writing the Geometry and Problem Files and Saving the Model
ELSA SPRAY MODELLING |
Chapter 23 |
Saving the Model File |
|
|
|
Saving the Model File
Having completed the es-ice part of the CFD model set-up, save your work up to this point in an es-ice model file.
•In the Select panel, click Write data
•In the Write Tool, enter save_es-ice.ELSA and click Save to save the model
•Close es-ice
STAR Set-up in pro-STAR
This section covers the required settings in pro-STAR, which is where the ELSA spray model, boundary conditions and STAR analysis controls are specified. Before attempting this part of the tutorial, it is important that you familiarize yourself with the pro-STAR interface by completing the example in Chapter 7 of this volume. In the following sections, most panel settings are only presented in summary form but information specific to ELSA and ECFM-CLEH is given in more detail.
The required steps are as follows:
1.Start up pro-STAR, open panel es-ice.PNL and use it to import the data created in es-ice
2.Activate the Lagrangian model for particle tracking
3.Define scalars required by the ELSA model
4.Set up the Lagrangian droplet physical models and controls
5.Define the injector boundary conditions
6.Set up the numerical analysis controls
7.Add an Extended Data segment to the problem file, as currently required by the ELSA model
Using the es-ice panel
Use the es-ice panel to import the mesh and physics settings created in es-ice and saved via the Write data operation in Star Controls. Also increase the pro-STAR memory allocation and define moving mesh events.
•Launch pro-STAR in the usual manner
•Select Panels > .es-ice from the menu bar to open the es-ice panel (see Figure 23-34)
•Click the Resize, Model and Events buttons in sequence
•Close the es-ice panel
23-24 |
Version 4.20 |
Chapter 23 |
ELSA SPRAY MODELLING |
|
STAR Set-up in pro-STAR |
|
|
Figure 23-34 es-ice panel
Activating the Lagrangian model
To set up an ELSA run, you must first activate the Lagrangian model in pro-STAR:
•In the pro-STAR Model Guide, select Analysis Features (see Figure 23-35)
•Set the Multi-Phase Treatment option to Lagrangian
•Click Apply
Figure 23-35 ELSA Analysis Features panel
Defining the ELSA scalars
Three ELSA-specific scalars are required for this case:
•LIQM_ELSA — The liquid phase of the fuel injected through the nozzle
•LSFD_ELSA — The liquid surface area density, used to determine the mean size of the liquid droplets (i.e. the droplet diameters)
•LIQV_ELSA — The vapour phase of the fuel taking part in combustion
Version 4.20 |
23-25 |
ELSA SPRAY MODELLING |
Chapter 23 |
STAR Set-up in pro-STAR |
|
|
|
The LIQM_ELSA and LSFD_ELSA scalars were defined as passive scalars (tracers) in the es-ice Star Controls panel. In this section, you must redefine LIQM_ELSA as an active scalar with properties that match the liquid fuel. The LSFD_ELSA scalar remains passive. You must also change the settings defining the mass diffusivity constant for both scalars.
A fuel vapour scalar (Scalar 1) has already been defined automatically by the ECFM-CLEH model, but you must change its name here to LIQV_ELSA so that it is compatible with ELSA’s naming conventions. Changing the fuel vapour scalar name also requires redefinition of the reaction scheme.
First, redefine the LIQM_ELSA scalar:
•In the pro-STAR Model Guide, select Thermophysical Models and Properties > Additional Scalars > Molecular Properties (Scalars) (see Figure 23-36)
•Set the Scalar # to 48
•Set Influence to Active
•Set the panel parameters as follows:
•Density (kg/m3) — 824.6
•Molecular Weight (kg/kmol) — 208.2
•Expansion Coefficient — 0
•Molecular Viscosity (kg/ms) — 0.00137563
•Specific Heat (J/kgK) — 2203.67
•Conductivity (W/mK) — 0.134905
•Heat of Formation (J/kg) — -1.69896e+06
•Temperature of Formation (K) — 298.15
•Click Apply and then click OK to close the warning message
23-26 |
Version 4.20 |
Chapter 23 |
ELSA SPRAY MODELLING |
|
STAR Set-up in pro-STAR |
|
|
Figure 23-36 LIQM_ELSA properties
•Click Define Polynomials
•In the Polynomial Function Definition panel (see Figure 23-37), set the panel parameters as follows:
•Tmin — 200
•Tmax — 750
•C1 — 6.8912
•C2 — 0.0895922
•C3 — 0.000294643
•C4 — -1.05188e-06
•C5 — 1.04676e-09
•Enthalpy — -42247.5
•Entropy — -48.5848
•Click Apply Coefficients followed by Close
Version 4.20 |
23-27 |
ELSA SPRAY MODELLING |
Chapter 23 |
STAR Set-up in pro-STAR |
|
|
|
Figure 23-37 LIQM_ELSA polynomial
Check the LSFD_ELSA scalar:
•Set the Scalar # to 49 (see Figure 23-38)
•Make sure that the Influence setting is Passive
•Click Apply
Figure 23-38 LSFD_ELSA properties
23-28 |
Version 4.20 |
Chapter 23 |
ELSA SPRAY MODELLING |
|
STAR Set-up in pro-STAR |
|
|
Define the LIQV_ELSA scalar:
•Enter the following commands to rename Scalar 1 (fuel vapour) as LIQV_ELSA
SCMODIFY, 1, NAME LIQV_ELSA
•Set the Scalar # to 1. Note that the name displayed on the panel has changed to LIQV_ELSA (see Figure 23-39)
Figure 23-39 LIQV_ELSA properties
Redefine the reaction scheme:
•Enter the following commands to define the leading reactant (i.e. LIQV_ELSA)
CRTYPE, 1 LREACT, 1
LIQV_ELSA
•Enter the following commands to define the reaction scheme
REACTION, 1 1 14.9 0 10.3 9.2
Version 4.20 |
23-29 |
ELSA SPRAY MODELLING |
Chapter 23 |
STAR Set-up in pro-STAR |
|
|
|
where:
1 — the number of LIQV_ELSA kilomoles in the reaction 14.9 — the number of Oxygen kilomoles in the reaction 0 — the number of Oxygen kilomoles in the products 10.3 — the number of CO2 kilomoles in the products 9.2 — the number of H2O kilomoles in the products
• Enter the following command for a stoichiometric check of the chemical reaction scheme
CRSCALAR, MAP, 1, DBASE
Define the binary properties of the ELSA scalars:
•In the pro-STAR Model Guide, select Additional Scalars > Binary Properties (see Figure 23-40)
•Select the Diffusion Velocity Correction check box
•In Scalar Mass Diffusivity, select Scalar Number 48 and set the Constant (m2/s) to 1e-10
•Click Apply
•Select Scalar Number 49 and set the Constant (m2/s) to 1e-10
•Click Apply
23-30 |
Version 4.20 |