Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dimensional analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

Nonlinear evolution of interacting oblique waves on two-dimensional shear layers

The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Goldstein, M. E.↗

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗

Rarefied gas flow through two-dimensional nozzles

A kinetic theory analysis is made of the flow of a rarefied gas from one reservoir to another through two-dimensional nozzles with arbitrary curvature. The Boltzmann equation simplified by a model collision integral is solved by means of finite-difference approximations with the discrete ordinate method. The physical space is transformed by a general grid generation technique and the velocity space is transformed to a polar coordinate system. A numerical code is developed which can be applied to any two-dimensional passage of complicated geometry for the flow regimes from free-molecular to slip. Numerical values of flow quantities can be calculated for the entire physical space including both inside the nozzle and in the outside plume. Predictions are made for the case of parallel slots and compared with existing literature data. Also, results for the cases of convergent or divergent slots and two-dimensional nozzles with arbitrary curvature at arbitrary knudsen number are presented.

De Witt, Kenneth J.↗

Efficient computation of the stability of three-dimensional compressible boundary layers

Methods for the computer analysis of the stability of three-dimensional compressible boundary layers are discussed and the user-oriented Compressible Stability Analysis (COSAL) computer code is described. The COSAL code uses a matrix finite-difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, M. R.↗

Stability of Thin Liquid Sheet Flows

A two-dimensional, linear stability analysis of a thin nonplanar liquid sheet flow in vacuum is carried out. A sheet flow created by a narrow slit of W and tau attains a nonplanar cross section as a consequence of cylinders forming on the sheet edge under the influence of surface tension forces. The region where these edge cylinders join the sheet is one of high curvature, and this is found to be the location where instability is most likely to occur. The sheet flow is found to be unstable, but with low growth rates for symmetric wave disturbances and high growth rates for antisymmetric disturbances. By combining the symmetric and antisymmetric disturbance modes, a wide range of stability characteristics is obtained. The product of unstable growth rate and flow time is proportional to the width-to-thickness ratio of the sift generating the sheet Three-dimensional effects can alter these results, particularly when the sheet length-to-width ratio is not much greater than unity.

McConley, Marc W.↗

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗

Hierarchical Strategy for Rapid Analysis Environment

A new philosophy is developed wherein the hierarchical definition of data is made use of in creating a better environment to conduct analyses of practical problems. This system can be adapted to conduct virtually any type of analysis, since this philosophy is not bound to any specific kind of analysis. It provides a framework to manage different models and its results and more importantly, the interaction between the different models. Thus, it is ideal for many types of finite element analyses like globalAoca1 analysis and those that involve multiple scales and fields. The system developed during the course of this work is just a demonstrator of the basic concepts. A complete implementation of this strategy could potentially make a major impact on the way analyses are conducted. It could considerably reduce the time frame required to conduct the analysis of real-life problems by efficient management of the data involved and reducing the human effort involved. It also helps in better decision making because of more ways to interpret the results. The strategy has been currently implemented for structural analysis, but with more work it could be extended to other fields of science when the finite element method is used to solve the differential equations numerically. This report details the work that has been done during the course of this project and its achievements and results. The following section discusses the meaning of the word hierarchical and the different references to the term in the literature. It talks about the development of the finite element method, its different versions and how hierarchy has been used to improve the methodology. The next section describes the hierarchical philosophy in detail and explains the different concepts and terms associated with it. It goes on to describe the implementation and the features of the demonstrator. A couple of problems are analyzed using the demonstrator program to show the working of the system. The two problems considered are two dimensional plane stress analysis problems. The results are compared with those obtained using conventional analysis. The different challenges faced during the development of this system are discussed. Finally, we conclude with suggestions for future work to add more features and extend it to a wider range of problems.

Whitcomb, John↗

Composite transport wing technology development

The design, fabrication, testing, and analysis of stiffened wing cover panels to assess damage tolerance criteria are discussed. The damage tolerance improvements were demonstrated in a test program using full-sized cover panel subcomponents. The panels utilized a hard skin concept with identical laminates of 44-percent 0-degree, 44-percent plus or minus 45-degree, and 12-percent 90-degree plies in the skins and stiffeners. The panel skins were impacted at midbay between the stiffeners, directly over the stiffener, and over the stiffener flange edge. The stiffener blades were impacted laterally. Impact energy levels of 100 ft-lb and 200 ft-lb were used. NASTRAN finite-element analyses were performed to simulate the nonvisible damage that was detected in the panels by nondestructive inspection. A closed-form solution for generalized loading was developed to evaluate the peel stresses in the bonded structure. Two-dimensional delamination growth analysis was developed using the principle of minimum potential energy in terms of closed-form solution for critical strain. An analysis was conducted to determine the residual compressive stress in the panels after impact damage, and the analytical predictions were verified by compression testing of the damaged panels.

Madan, Ram C.↗

The two-dimensional diagnostic model for tropospheric OH - An uncertainty analysis

Tropospheric OH is believed to play a major role in atmospheric photochemistry, and, via its chemical interactions OH ultimately may affect the climate, stratospheric ozone levels, and the pH of rain. A valuable tool for predicting the global OH abundance as a function of latitude and altitude is the two-dimensional diagnostic model for tropospheric OH. This model can also be applied to budget studies of species which interact with OH, such as CH4, CO, and O3. In addition to laboratory-obtained kinetic data, important variables include the rate constants for several reactions, the rate of heterogeneous removal of soluble species, and the global abundance and spatial variability in H2O, O3, CO, and NOx. The uncertainties in these parameters imply over a factor of 3 variability in the calculated global OH concentration as well as significant uncertainties in the latitudinal and vertical OH distribution.

Chameides, W. L.↗

Two-Dimensional Array Beam Scanning Via Externally and Mutually Injection Locked Coupled Oscillators

Some years ago, Stephan proposed an approach to one dimensional (linear) phased array beam steering which requires only a single phase shifter. This involves the use of a linear array of voltage-controlled electronic oscillators coupled to nearest neighbors. The oscillators are mutually injection locked by controlling their coupling and tuning appropriately. Stephan's approach consists of deriving two signals from a master oscillator, one signal phase shifted with respect to the other by means of a single phase shifter. These two signals are injected into the end oscillators of the array. The result is a linear phase progression across the oscillator array. Thus, if radiating elements are connected to each oscillator and spaced uniformly along a line, they will radiate a beam at an angle to that line determined by the phase gradient which is, in turn, determined by the phase difference between the injection signals.The beam direction is therefore controlled by adjusting this phase difference. Recently, Pogorzelski and York presented a formulation which facilitates theoretical analysis of the above beam steering technique. This was subsequently applied by Pogorzelski in analysis of two dimensional beam steering using perimeter detuning of a coupled oscillator array. The formulation is based on a continuum model in which the oscillator phases are represented by a continuous function satisfying a partial differential equation of diffusion type. This equation can be solved via the Laplace transform and the resulting solution exhibits the dynamic behavior of the array as the beam is steered. Stephan's beam steering technique can be similarly generalized to two-dimensional arrays in which the beam control signals are applied to the oscillators on the perimeter of the array. In this paper the continuum model for this two-dimensional case is developed and the dynamic solution for the corresponding aperture phase function is obtained. The corresponding behavior of the resulting far-zone radiation pattern is displayed as well.

Pogorzelski, Ronald J.↗

Differentiating Double-Layer, Psuedocapacitance, and Battery-like Mechanisms by Analyzing Impedance Measurements in Three Dimensions

Electrochemical energy storage arises from processes that are broadly categorized as capacitive, pseudocapacitive, or battery-like. Advanced charge-storing materials that are designed to deliver high capacity at a high rate often exhibit a multiplicity of such mechanisms, which complicates the understanding of their charge-storage behavior. In this paper, we apply a “3D Bode analysis” technique to identify key descriptors for fast Li-ion storage processes, where AC impedance data, such as the real capacitance ($C'$) or phase angle ($Φ$), are represented versus the frequency ($f$) and a third independent variable, the applied DC cell voltage. For double-layer processes, a near-constant $C'$ or $Φ$ is supported across the entire voltage range, and the decrease in these values shows a near-linear decrease at higher $f$. For pseudocapacitance, an increase in $C'$ is delivered, accompanied by high $C'$ retention at higher $f$ compared to double-layer processes. Interestingly, the lower $Φ$ values, where $C'$ is highest, suggest that this is a key descriptor for pseudocapacitance, where high-rate charge storage is still facilitated within a kinetically limited regime. For battery-like processes, a high $C'$ is only observed at the voltage at which the material stores charge, while outside that voltage, $C'$ is negligible. The three-dimensional (3D) Bode analysis allows charge-storage dynamics to be mapped out in great detail with more delineation between mechanisms compared to the more frequently deployed kinetic analyses derived from cyclic voltammetry.

25 ENERGY STORAGE↗

Strain-energy-release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtural crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finite-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

Strain energy release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtual crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

Stability analysis for laminar flow control, part 1

The basic equations for the stability analysis of flow over three dimensional swept wings are developed and numerical methods for their solution are surveyed. The equations for nonlinear stability analysis of three dimensional disturbances in compressible, three dimensional, nonparallel flows are given. Efficient and accurate numerical methods for the solution of the equations of stability theory were surveyed and analyzed.

Benney, D. J.↗

Three-Dimensional Unsteady Simulation of a Modern High Pressure Turbine Stage Using Phase Lag Periodicity: Analysis of Flow and Heat Transfer

Unsteady three-dimensional RANS simulations have been performed on a highly loaded transonic turbine stage and results are compared to steady calculations as well as experiment. A low Reynolds number k- turbulence model is employed to provide closure for the RANS system. A phase-lag boundary condition is used in the periodic direction. This allows the unsteady simulation to be performed by using only one blade from each of the two rows. The objective of this paper is to study the effect of unsteadiness on rotor heat transfer and to glean any insight into unsteady flow physics. The role of the stator wake passing on the pressure distribution at the leading edge is also studied. The simulated heat transfer and pressure results agreed favorably with experiment. The time-averaged heat transfer predicted by the unsteady simulation is higher than the heat transfer predicted by the steady simulation everywhere except at the leading edge. The shock structure formed due to stator-rotor interaction was analyzed. Heat transfer and pressure at the hub and casing were also studied. Thermal segregation was observed that leads to the heat transfer patterns predicted by steady and unsteady simulations to be different.

Shyam, Vikram↗

Finding Hidden Patterns in High Resolution Wind Flow Model Simulations

Wind flow data is critical in terms of investment decisions and policy making. High resolution data from wind flow model simulations serve as a supplement to the limited resource of original wind flow data collection. Given the large size of data, finding hidden patterns in wind flow model simulations are critical for reducing the dimensionality of the analysis. In this work, we first perform dimension reduction with two autoencoder models: the CNN-based autoencoder (CNN-AE) [1], and hierarchical autoencoder (HIER-AE) [2], and compare their performance with the Principal Component Analysis (PCA). We then investigate the super-resolution of the wind flow data. By training a Generative Adversarial Network (GAN) with 300 epochs, we obtained a trained model with 2× resolution enhancement. We compare the results of GAN with Convolutional Neural Network (CNN), and GAN results show finer structure as expected in the data field images. Also, the kinetic energy spectra comparisons show that GAN outperforms CNN in terms of reproducing the physical properties for high wavenumbers and is critical for analysis where high-wavenumber kinetics play an important role.

97 MATHEMATICS AND COMPUTING↗

Effect of heat treatment on microstructural evolution and hardness homogeneity in laser powder bed fusion of alloy 718

The thermal history developed in laser powder bed fusion (LPBF) processes has been shown to be complex resulting in equally complex microstructures and mechanical properties. Microscopic observations and Vickers micro-hardness mapping measurements were carried out on diff ;erent section planes of LPBF alloy 718 cuboids. Three-dimensional finite element analysis was used to simulate thermal history and to predict the residual stress distribution in the as-built material. Computational thermodynamics was used to predict the micro-segregation and nucleation driving force of various phases in the bulk and in segregated regions. Varied heat-treatments such as simulated hot isostatic pressing, and double aging were applied. Their influence on the microstructure, microsegregation, precipitate formation, and micro-hardness variations of LPBF alloy 718 were investigated. Hardness map results showed heterogeneous micro-hardness on the xy- and xz-planes of the as-built parts where the bottom plane and center regions had larger hardness of ~315 HV 0.5 while the top plane and contours showed hardness of ~300 HV 0.5 . It was found that the aging treatment increased the overall hardness of the as-built condition from ~310 HV 0.5 to 470 HV 0.5 but also increased the hardness gradient throughout the coupon. After simulated hot isostatic pressing process (i.e., without applied pressure) at 1020 °C for 4 h followed by water quench (HIPWQ), the hardness gradient and hardness was minimized (~210 HV 0.5 ) as the microstructure transitioned from heterogeneous columnar grains in the as-built condition to more uniform recrystallized grains. A double aging treatment was applied to enhance hardness from ~210 HV 0.5 to ~440 HV 0.5 . Finally, HIPWQ followed by double aging produced a homogeneous microstructure and more uniform hardness map with enhanced mechanical properties in LPBF alloy 718 coupons.

36 MATERIALS SCIENCE↗

3D pore structure reconstruction of fluff region in top of EBR-II irradiated fuel by X-ray micro-computed tomography

Metallic fuels comprised of U-Zr binary and U-Pu-Zr ternary alloys are a promising fuel candidate for sodium-cooled fast reactors and possess many advantages over ceramic fuels such as higher uranium density and improved thermal conductivity. Fluff is a highly porous structure forming at the top of metallic fuel slugs during irradiation and was extensively observed in past sodium-cooled fast reactor (SFR) experiments such as Idaho National Laboratory’s Experimental Breeder Reactor II (EBR-II). Previous investigations have shown that fluff is rich in fissile atoms and hence could have neutronics implications, but its formation mechanisms are not understood. In this work, two cube-shaped lift outs were taken from the fluff and bulk fuel regions of a high burnup U-19Pu-10Zr composition EBR-II fuel slug and characterized by micro computed tomography (Micro-CT) and energy dispersive x-ray spectroscopy (EDS). This was used to characterize differences in microscale porosity and composition between fluff and bulk fuel, and to support or oppose postulated fluff formation mechanisms. Surface EDS showed minimal differences between the elemental compositions of the fluff and fuel cubes. However, three-dimensional (3D) pore analysis from the tomography scans showed significant differences between the fluff and fuel cube pore size distributions at the microscale level. The high porosity and irregular pore distribution seen in macroscale fluff was also present within the micro-scale fluff cube, having a much higher porosity than the fuel cube which was concentrated in a few large pore networks. In conclusion, these results indicate support for a creep driven mechanism of fluff formation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗