CFD Predictions of Soot & CO Emissions Generated by a Partially-Fueled 9-Element Lean-Direct Injection Combustor

A study was undertaken to investigate the CO & soot emissions generated by a partially-fueled 9- element LDI (Lean-Direct Injection) combustor configuration operating in the idle range of jet engine conditions. In order to perform the CFD analysis, several existing soot/chemistry models were implemented into the OpenNCC (Open National Combustion Code). The calculations were based on a Reynolds-Averaged Navier Stokes (RANS) simulation with standard k-epsilon turbulence model, a 62- species jet-a/air chemistry, a 2-equation soot model, & a Lagrangian spray solver. A separate transport equation was solved for all individual species involved in jet-a/air combustion. In the test LDI configuration we examined, only five of the nine injectors were fueled with the major pilot injector operating at an equivalence ratio of near one and the other four main injectors operating at an equivalence ratio near 0.55. The calculations helped to identify several reasons behind the soot & CO formation in different regions of the combustor. The predicted results were compared with the reported experimental data on soot mass concentration (SMC) & emissions index of CO (EICO). The experimental results showed that an increase in either T3 and/or F/A ratio lead to a reduction in both EICO & SMC. The predicted results were found to be in reasonable agreement. However, the predicted EICO differed substantially in one test condition associated with higher F/A ratio

Stress-intensity factor equations for cracks in three-dimensional finite bodies subjected to tension and bending loads

Stress intensity factor equations are presented for an embedded elliptical crack, a semielliptical surface crack, a quarter elliptical corner crack, a semielliptical surface crack along the bore of a circular hole, and a quarter elliptical corner crack at the edge of a circular hole in finite plates. The plates were subjected to either remote tension or bending loads. The stress intensity factors used to develop these equations were obtained from previous three dimensional finite element analyses of these crack configurations. The equations give stress intensity factors as a function of parametric angle, crack depth, crack length, plate thickness, and, where applicable, hole radius. The ratio of crack depth to plate thickness ranged from 0 to 1, the ratio of crack depth to crack length ranged from 0.2 to 2, and the ratio of hole radius to plate thickness ranged from 0.5 to 2. The effects of plate width on stress intensity variation along the crack front were also included

Three dimensional finite-element analysis of finite-thickness fracture specimens

The stress-intensity factors for most of the commonly used fracture specimens (center-crack tension, single and double edge-crack tension, and compact), those that have a through-the-thickness crack, were calculated using a three dimensional finite-element elastic stress analysis. Three-dimensional singularity elements were used around the crack front. The stress intensity factors along the crack front were evaluated by using a force method, developed herein, that requires no prior assumption of either plane stress or plane strain. The calculated stress-intensity factors from the present analysis were compared with those from the literature whenever possible and were generally found to be in good agreement. The stress-intensity factors at the midplane for all specimens analyzed were within 3 percent of the two dimensional plane strain values. The stress intensity factors at the specimen surfaces were considerably lower than at the midplanes. For the center-crack tension specimens with large thickness to crack-length ratios, the stress-intensity factor reached a maximum near the surface of the specimen. In all other specimens considered the maximum stress intensity occurred at the midplane

Improved stress-intensity factors for semi-elliptical surface cracks in finite-thickness plates

Stress-intensity factors for shallow and deep semi-elliptical surface cracks in plates subjected to tension are presented. To verify the accuracy of the three-dimensional finite-element models employed, convergence was studied by varying the number of degrees of freedom in the models from 1500 to 6900. The 6900 degrees of freedom used here were more than twice the number used in previously reported solutions

surf3d: A 3-D finite-element program for the analysis of surface and corner cracks in solids subjected to mode-1 loadings

A computer program, surf3d, that uses the 3D finite-element method to calculate the stress-intensity factors for surface, corner, and embedded cracks in finite-thickness plates with and without circular holes, was developed. The cracks are assumed to be either elliptic or part eliptic in shape. The computer program uses eight-noded hexahedral elements to model the solid. The program uses a skyline storage and solver. The stress-intensity factors are evaluated using the force method, the crack-opening displacement method, and the 3-D virtual crack closure methods. In the manual the input to and the output of the surf3d program are described. This manual also demonstrates the use of the program and describes the calculation of the stress-intensity factors. Several examples with sample data files are included with the manual. To facilitate modeling of the user's crack configuration and loading, a companion program (a preprocessor program) that generates the data for the surf3d called gensurf was also developed. The gensurf program is a three dimensional mesh generator program that requires minimal input and that builds a complete data file for surf3d. The program surf3d is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The stress intensity factors were calculated by a three dimensional finite element method. The finite element models employed singularity elements along the crack front and linear strain elements elsewhere. The models had about 6500 degrees of freedom. The stress intensity factors were evaluated from a nodal force method. An equation for the stress intensity factors was obtained form the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary integral equation method were in agreement (+ or - 2 percent) and those from another finite element were in fair agreement (+ or - 8 percent) with the present results

Stress-intensity factor equations for cracks in three-dimensional finite bodies

Empirical stress intensity factor equations are presented for embedded elliptical cracks, semi-elliptical surface cracks, quarter-elliptical corner cracks, semi-elliptical surface cracks at a hole, and quarter-elliptical corner cracks at a hole in finite plates. The plates were subjected to remote tensile loading. Equations give stress intensity factors as a function of parametric angle, crack depth, crack length, plate thickness, and where applicable, hole radius. The stress intensity factors used to develop the equations were obtained from three dimensional finite element analyses of these crack configurations

Prediction of fatigue crack-growth patterns and lives in three-dimensional cracked bodies

Fatigue crack growth patterns and lives for surface cracks, surface cracks at holes, and corner cracks at holes in three dimensional bodies were predicted using linear-elastic fracture mechanics concepts that were modified to account for crack-closure behavior. The predictions were made by using stress intensity factor equations for these crack configurations and the fatigue crack-growth (delta K against rate) relationship for the material of interest. The crack configurations were subjected to constant-amplitude fatigue loading under either remote tension or bending loads. The predicted crack growth patterns and crack growth lives for aluminum alloys agreed well with test data from the literature

Stress-Intensity Factors for Corner Cracks at the Edge of a Hole

Stress-intensity factors, calculated by a three-dimensional finite-element analysis, for shallow or deep quarter-elliptical corner cracks at the edge of a hole in a finite-thickness plate are presented. The plate was subjected to remote uniform tension, remote bending, or simulated pin loading in the hole. The crack depth-to-plate thickness ranged from 0.2 to 0.8, while the ratio of crack depth-to-plate crack length ranged from 0.2 to 2. The ratio of hole radius-to-plate thickness was held at 0.5. To verify the accuracy of the three-dimensional finite-element models empolyed, convergence studies were conducted (number of degrees of freedom ranged from 4400 to 9300). The stress-intensity factor variations along the crack front are presented and compared with other solutions from the literature