Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary element method”

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 865 records · Page 48

Heat transfer in the turbulent boundary layer with a step change in surface roughness

Measurements of Stanton numbers, velocity profiles, temperature profiles, and turbulence intensity profiles are reported for turbulent flat plate boundary layer flows with a step change in surface roughness. The first 0.9 m length of the test surface is roughened with 1.27 mm diameter hemispheres spaced 2 base diameters apart in a staggered array. The remaining 1.5 m length is smooth. The experiments show that the step change from a rough to a smooth surface has a dramatic effect on the convective heat transfer. In many cases, the Stanton number drops below the smooth-wall correlation immediately downstream of the change in roughness. The Stanton number measurements are compared with predictions using the discrete element method with excellent results.

Taylor, Robert P.↗

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗

Two-dimensional finite-element temperature variance analysis

The finite element method is extended to thermal analysis by forming a variance analysis of temperature results so that the sensitivity of predicted temperatures to uncertainties in input variables is determined. The temperature fields within a finite number of elements are described in terms of the temperatures of vertices and the variational principle is used to minimize the integral equation describing thermal potential energy. A computer calculation yields the desired solution matrix of predicted temperatures and provides information about initial thermal parameters and their associated errors. Sample calculations show that all predicted temperatures are most effected by temperature values along fixed boundaries; more accurate specifications of these temperatures reduce errors in thermal calculations.

Heuser, J. S.↗

Reynolds number influences in aeronautics

Reynolds number, a measure of the ratio of inertia to viscous forces, is a fundamental similarity parameter for fluid flows and therefore, would be expected to have a major influence in aerodynamics and aeronautics. Reynolds number influences are generally large, but monatomic, for attached laminar (continuum) flow; however, laminar flows are easily separated, inducing even stronger, non-monatomic, Reynolds number sensitivities. Probably the strongest Reynolds number influences occur in connection with transitional flow behavior. Transition can take place over a tremendous Reynolds number range, from the order of 20 x 10(exp 3) for 2-D free shear layers up to the order of 100 x 10(exp 6) for hypersonic boundary layers. This variability in transition behavior is especially important for complex configurations where various vehicle and flow field elements can undergo transition at various Reynolds numbers, causing often surprising changes in aerodynamics characteristics over wide ranges in Reynolds number. This is further compounded by the vast parameterization associated with transition, in that any parameter which influences mean viscous flow development (e.g., pressure gradient, flow curvature, wall temperature, Mach number, sweep, roughness, flow chemistry, shock interactions, etc.), and incident disturbance fields (acoustics, vorticity, particulates, temperature spottiness, even electro static discharges) can alter transition locations to first order. The usual method of dealing with the transition problem is to trip the flow in the generally lower Reynolds number wind tunnel to simulate the flight turbulent behavior. However, this is not wholly satisfactory as it results in incorrectly scaled viscous region thicknesses and cannot be utilized at all for applications such as turbine blades and helicopter rotors, nacelles, leading edge and nose regions, and High Altitude Long Endurance and hypersonic airbreathers where the transitional flow is an innately critical portion of the problem.

Bushnell, Dennis M.↗

Micromechanical models for textile composites

Numerical and analytical micromechanical models are presented to predict the thermoelastic behavior of a textile composite. In the numerical model, the unit-cell is discretized with finite elements, and periodic boundary conditions are imposed between opposite faces of the unit-cell. For a thin textile composite, stress gradients effects through the thickness are demonstrated. The consequent difference in the stiffness and strength behavior of thick and thin composites are discussed. The numerical model is implemented to predict 3-D thermo-elastic constants for a thick textile composite, and the plate thermo-mechanical properties for a thin textile composite. The numerical model is extended to compute the thermal residual microstresses due to processing to predict the composite failure envelopes. An analytical model - Selective Averaging Method (SAM) - is proposed, which is based on a judicious combination of stiffness and compliance averaging to estimate the 3-D elastic constants. Both the models are tested and verified for several examples by comparing the stiffness properties with elasticity solutions and available results.

Sankar, Bhavani V.↗

Techniques Optimized for Reducing Instabilities in Advanced Nickel-Base Superalloys for Turbine Blades

The High-Speed Research (HSR) Airfoil Alloy program developed fourth-generation single-crystal superalloys with up to an 85 F increase in creep rupture capability over current production airfoil alloys. Recent results have been generated at the NASA Glenn Research Center on these fourth-generation alloys, but in coated form, for subsonic turbine blade applications under NASA's Ultra-Efficient Engine Technology (UEET) Program. One goal for UEET is to optimize the airfoil alloy/thermal barrier coating system for 3100 F turbine inlet temperatures. The state-of-the art turbine blade airfoil system consists of a superalloy single crystal that provides the basic mechanical performance of the airfoil. A thermal barrier coating is used to reduce the temperature of the base superalloy, and a bondcoat is deposited between the base material and the thermal barrier coating. The bondcoat improves the oxidation and corrosion resistance of the base superalloy and improves the spallation resistance of the thermal barrier coating. A commercial platinum aluminide bondcoat was applied to the HSR-developed alloys, and a diffusion zone developed as a result of interaction between the bondcoat and the superalloy. Optimized strength is obtained for superalloys when the refractory element content is high and the limits of microstructural stability are approached or exceeded slightly. For fourthgeneration alloys, instability leads to the formation of topologically close packed (TCP) phases, which form internally in the superalloy, and a secondary reaction zone (SRZ), which forms under the diffusion zone. There was a concern that excessive quantities of either TCP or SRZ might decrease the mechanical properties of the superalloy, with SRZ thought to be particularly detrimental and its formation unpredictable. Thus, an SRZreduction effort was initiated in the NASA UEET Program so that methods developed during the HSR project could be optimized further to reduce or eliminate the SRZ. An SRZ is a three-phase constituent composed of TCP and stringers of gamma phase in a matrix of gamma prime. An incoherent grain boundary separates the SRZ from the gammagamma prime microstructure of the superalloy. The SRZ is believed to form as a result of local chemistry changes in the superalloy due to the application of the diffusion aluminide bondcoat. Locally high surface stresses also appear to promote the formation of the SRZ. Thus, techniques that change the local alloy chemistry or reduce surface stresses have been examined for their effectiveness in reducing SRZ. These SRZ-reduction steps are performed on the test specimen or the turbine blade before the bondcoat is applied. Stressrelief heat treatments developed at NASA Glenn have been demonstrated to reduce significantly the amount of SRZ that develops during subsequent high-temperature exposures. Stress-relief heat treatments reduce surface stresses by recrystallizing a thin surface layer of the superalloy. However, in alloys with very high propensities to form SRZ, stress relief heat treatments alone do not eliminate SRZ entirely. Thus, techniques that modify the local chemistry under the bondcoat have been emphasized and optimized successfully at Glenn. One such technique is carburization, which changes the local chemistry by forming submicron carbides near the surface of the superalloy. Detailed characterizations have demonstrated that the depth and uniform distribution of these carbides are enhanced when a stress relief treatment and an appropriate surface preparation are employed in advance of the carburization treatment. Even in alloys that have the propensity to develop a continuous SRZ layer beneath the diffusion zone, the SRZ has been completely eliminated or reduced to low, manageable levels when this combination of techniques is utilized. Now that the techniques to mitigate SRZ have been established at Glenn, TCP phase formation is being emphasized in ongoing work under the UEET Program. The limitsf stability of the fourth-generation alloys with respect to TCP phase formation are currently being defined along with high-temperature creep rupture properties. In addition, a regression model is being developed at Glenn for the prediction of the presence of TCP phase in the microstructure and SRZ under the diffusion zone. The model is based on a design-of-experiments methodology with emphasis on the potential synergistic effects of alloying elements.

MacKay, Rebecca A.↗

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow; (2) the energy equation for the structural-thermal response; and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two application studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗

Evaluation of an adaptive unstructured remeshing technique for integrated fluid-thermal-structural analysis

An adaptive unstructured remeshing technique is evaluated for integrated fluid-thermal-structural analysis. The technique is combined with the finite element method to solve: (1) the Navier-Stokes equations for high-speed compressible flow, (2) the energy equation for the structural-thermal response, and (3) the quasi-static equilibrium equations for the structural response. The remeshing technique and the analysis solution procedure are described. The effectiveness of the approach is evaluated with two applications studies. The flow analysis of Mach 8 shock-shock interference on a three-inch-diameter cylinder is used as the first application study to demonstrate the capability of the remeshing technique and to examine proper remeshing indicators for the inviscid and boundary layer regions. The applicability of the approach for the thermal and structural analyses of the structure is evaluated in the second application study of a 0.25-inch-diameter convectively cooled leading edge subjected to intense aerodynamic heating. Issues associated with remeshing indicators for thermal stress problems are identified.

Dechaumphai, Pramote↗

Theory and Experiment of Multielement Airfoils: A Comparison

A detailed comparison of computed and measured pressure distributions, velocity profiles, transition onset, and Reynolds shear stresses for multi-element airfoils is presented. It is shown that the transitional k-zeta model, which is implemented into CFL3D, does a good job of predicting pressure distributions, transition onset, and velocity profiles with the exception of velocities in the slat wake region. Considering the fact that the hot wire used was not fine enough to resolve Reynolds stresses in the boundary layer, comparisons of turbulence stresses varied from good to fair. It is suggested that the effects of unsteadiness be thoroughly evaluated before more complicated transition/turbulence models are used. Further, it is concluded that the present work presents a viable and economical method for calculating laminar/transitional/turbuient flows over complex shapes without user interface.

Czerwiec, Ryan↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Effects of local vibrations on the dynamics of space truss structures

The paper discusses the influence of local member vibrations on the dynamics of repetitive space truss structures. Several focus problems wherein local member vibration modes are in the frequency range of the global truss modes are discussed. Special attention is given to defining methods that can be used to identify the global modes of a truss structure amidst many local modes. Significant interactions between the motions of local member vibrations and the global behavior are shown to occur in truss structures when: (1) the natural frequencies of the individual members for clamped-clamped boundary conditions are in the vicinity of the global truss frequency; and (2) the total mass of the individual members represents a large portion of the mass of the whole structure. The analysis is carried out with a structural analysis code which uses exact member theory. The modeling detail required using conventional finite element codes to adequately represent such a class of problems is examined. The paper concludes with some practical considerations for the design and dynamic testing of structures which might exhibit such behavior.

Warnaar, Dirk B.↗

Preliminary large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids

Large-eddy simulation (LES) has matured to the point where application to complex flows is desirable. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential. In the past two years, important building blocks have been put into place making possible a careful study of LES on unstructured grids. The first building block was an efficient mesh generator which allowed the placement of points according to smooth variation of physical length scales. This variation of length scales is in all three directions independently, which allows a large reduction in points when compared to structured-grid methods, which can only vary length scales in one direction at a time. The second building block was the development of a dynamic model appropriate for unstructured grids. The principle obstacle was the development of an unstructured-grid filtering operator. In the past year, some of the new filters developed by Jansen have been implemented into a highly parallelized finite element code based on the Galerkin/least-squares finite element method. We have chosen the NACA 4412 airfoil at maximum lift as the first simulation for a variety of reasons. First, it is a problem of significant interest since it would be the first LES of an aircraft component. Second, this flow has been the subject of three experimental studies. The third reason for considering this flow is the variety of flow features which provide an important test of the dynamic model. Only the dynamic model can be expected to perform satisfactorily in this variety of situations: from the laminar regions where it must not modify the flow at all to the turbulent boundary layers and wake where it must represent a wide variety of subgrid-scale structures. The flow configuration we have chosen is that of Wadcock (1987) at Reynolds number based on chord Re(sub c) = u(sub infinity)c/v = 1.64 x 10(exp 6), Mach number M = 0.2, and 12 deg angle of attack.

Jansen, Kenneth↗

Graphical User Interface Development for Representing Air Flow Patterns

In the Turbine Branch, scientists carry out experimental and computational work to advance the efficiency and diminish the noise production of jet engine turbines. One way to do this is by decreasing the heat that the turbine blades receive. Most of the experimental work is carried out by taking a single turbine blade and analyzing the air flow patterns around it, because this data indicates the sections of the turbine blade that are getting too hot. Since the cost of doing turbine blade air flow experiments is very high, researchers try to do computational work that fits the experimental data. The goal of computational fluid dynamics is for scientists to find a numerical way to predict the complex flow patterns around different turbine blades without physically having to perform tests or costly experiments. When visualizing flow patterns, scientists need a way to represent the flow conditions around a turbine blade. A researcher will assign specific zones that surround the turbine blade. In a two-dimensional view, the zones are usually quadrilaterals. The next step is to assign boundary conditions which define how the flow enters or exits one side of a zone. way of setting up computational zones and grids, visualizing flow patterns, and storing all the flow conditions in a file on the computer for future computation. Such a program is necessary because the only method for creating flow pattern graphs is by hand, which is tedious and time-consuming. By using a computer program to create the zones and grids, the graph would be faster to make and easier to edit. Basically, the user would run a program that is an editable graph. The user could click and drag with the mouse to form various zones and grids, then edit the locations of these grids, add flow and boundary conditions, and finally save the graph for future use and analysis. My goal this summer is to create a graphical user interface (GUI) that incorporates all of these elements. I am writing the program in Java, a language that is portable among platforms, because it can run on different operating systems such as Windows and Unix without having to be rewritten. I had no prior experience of programming in Java at the start of my internship; I am continuously learning as I create the program. I have written the part of the program that enables a user to draw several zones, edit them, and store their locations. The next phase of my project is to allow the user to click on the side of a zone and create a boundary condition for it. A previous intern wrote a program that allows the user to input boundary conditions. I can integrate the two programs to create a larger, more usable program. After that, I will develop a way for the user to save the graph for future reference. Another eventual goal is to make the GUI capable of creating three-dimensional zones as well. Researchers such as my mentor, Dr. David Ashpis, need a quick, user-friendly

Chaudhary, Nilika↗

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

Usage of ChatGPT for Engineering Design and Analysis Tool Development

ChatGPT, a generative AI large language model, has recently captured significant attention in both the computer science community and the broader public domain. It has demonstrated a wide range of capabilities, from answering simple questions to writing fully functional computer code. This study spotlights both the capabilities and limitations of ChatGPT when addressing engineering problems. The model's capacity to generate practical engineering tools is highlighted through an example of a prompt that leads to an interactive plotting tool, enabling the examination of the fluid boundary layer around a fan blade. Subsequently, the paper also uncovers potential pitfalls in ChatGPT’s application, shown through an unsuccessful attempt to use ChatGPT to automate a process in Ansys Workbench through scripting. The research further investigates ChatGPT's proficiency in addressing inquiries and providing explanations about the functionalities of OpenMDAO, an open-source, multidisciplinary design, analysis, and optimization tool developed at NASA Glenn Research Center. Finally, an optimization methodology, developed with ChatGPT’s help, is applied to the structural optimization of a fan blade. The developed optimization method utilizes T-Blade3 for geometry generation, Ansys Mechanical for meshing and finite element analysis, and sci-kit learn’s MLPRegressor method to generate a trained neural network model of the design space. OpenMDAO is then used to find the optimal point within the design space. The outcome is a significant reduction in stress in the optimized model—less than one-fifth of the stress value in the baseline model.

Design↗

Usage of ChatGPT for Engineering Design and Analysis Tool Development

ChatGPT, a generative AI large language model, has recently captured significant attention in both the computer science community and the broader public domain. It has demonstrated a wide range of capabilities, from answering simple questions to writing fully functional computer code. This study spotlights both the capabilities and limitations of ChatGPT when addressing engineering problems. The model's capacity to generate practical engineering tools is highlighted through an example of a prompt that leads to an interactive plotting tool, enabling the examination of the fluid boundary layer around a fan blade. Subsequently, the paper also uncovers potential pitfalls in ChatGPT’s application, shown through an unsuccessful attempt to use ChatGPT to automate a process in Ansys Workbench through scripting. The research further investigates ChatGPT's proficiency in addressing inquiries and providing explanations about the functionalities of OpenMDAO, an open-source, multidisciplinary design, analysis, and optimization tool developed at NASA Glenn Research Center. Finally, an optimization methodology, developed with ChatGPT’s help, is applied to the structural optimization of a fan blade. The developed optimization method utilizes T-Blade3 for geometry generation, Ansys Mechanical for meshing and finite element analysis, and sci-kit learn’s MLPRegressor method to generate a trained neural network model of the design space. OpenMDAO is then used to find the optimal point within the design space. The outcome is a significant reduction in stress in the optimized model—less than one-fifth of the stress value in the baseline model.

Design↗

System-state and operating condition sensitive control method and apparatus for electric power delivery systems

This invention provides a method and apparatus for determining a precise switching sequence for the power switching elements of electric power delivery systems of the on-off switching type and which enables extremely fast transient response, precise regulation and highly stable operation. The control utilizes the values of the power delivery system power handling network components, a desired output characteristic, a system timing parameter, and the externally imposed operating conditions to determine where steady state operations should be in order to yield desired output characteristics for the given system specifications. The actual state of the power delivery system is continuously monitored and compared to a state-space boundary which is derived from the desired equilibrium condition, and from the information obtained from this comparison, the system is moved to the desired equilibrium condition in one cycle of switching control. Since the controller continuously monitors the power delivery system's externally imposed operating conditions, a change in the conditions is immediately sensed and a new equilibrium condition is determined and achieved, again in a single cycle of switching control.

Burns, III, William Wesley↗