Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary”

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 217 records · Page 12

An Experimental and Numerical Study of Roughness-Induced Instabilities in a Mach 3.5 Boundary Layer

Progress on a joint experimental and numerical study of laminar-to-turbulent transition induced by an isolated roughness element in a high-speed laminar boundary layer is reported in this paper. The numerical analysis suggests that transition is driven by the instability of high- and low-speed streaks embedded in the wake of the isolated roughness element. In addition, spatial stability analysis revealed that the wake flow supports multiple modes (even and odd) of convective instabilities that experience strong enough growth to cause transition. The experimental measurements, which included hot-wire and pitot-probe surveys, confirmed the existence of embedded high- and low-speed streaks in the roughness wake. Furthermore, the measurements indicate the presence of both even and odd modes of instability, although their relative magnitude depends on the specifics of the roughness geometry and flow conditions (e.g., the value of Re(sub kk) or k/delta. For the two test cases considered in the measurements (Re(sub kk) values of 462 and 319), the even mode and the odd mode were respectively dominant and appear to play a primary role in the transition process.

Kegerise, Michael A.↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗

Embedded shear layer computations for increased drag reduction

One of the most promising methods of minimizing drag is the reduction of skin friction by injection of low momentum fluid into the near-wall region of turbulent boundary layer flows. This method could be made more effective by limiting the spread rate of the resulting mixing region. In order to achieve a better understanding of how this goal might be achieved, numerical investigations of the relevant fluid dynamic processes governing these regions have been conducted. A compact finite-difference algorithm has been applied to the complete form of the governing conservation equations for a two-dimensional laminar mixing layer. The ability of this computational approach to model successfully the formation and interaction of the large scale vortical structures which dominate such flow fields is verified in the present study. Parameters which affect the spread rate of the mixing region are also identified. In addition, the relative importance of viscous and momentum transport effects in the vortex interactions is determined.

Gatski, T. B.↗

On sound radiation from the trailing edge of an isolated airfoil in a uniform flow

An experimental investigation has been conducted to determine the physical process of the noise production from the trailing edge of an isolated two-dimensional airfoil embedded in a low turbulence uniform mean flow. The Reynolds number of the airfoil based on chord was greater than 1,000,000 and that the boundary layer was fully turbulent at the trailing edge. Smoke injection technique and spark shadowgraphy were used to study the structural features of the boundary layer. Two point joint statistical analyses were made on the surface pressure, the radiated sound and the relationship between the two quantities. Operating conditions included two free stream velocities and two angles of attack for both a naturally developed and an artificially tripped boundary layer. Flow visualizations revealed the existence of large scale coherent motions in the outer region of the boundary layer. The observed coherent motions had scales of the order of boundary layer thickness and a convection velocity near the free stream velocity. The production process of the airfoil trailing edge noise was determined to be the convection of the large scale coherent eddies over the trailing edge of the airfoil. The noise field so generated was found to be rather coherent and dipole-like.

Yu, J. C.↗

Igneous petrology of the new ureilites Nova 001 and Nullarbor 010

The Nova 001 (= Nuevo Mercurio (b)) and Nullarbor 010 meteorites are ureilites, both of which contain euhedral graphite crystals. The bulk of the meteorites are olivine (Fo79) and pyroxenes (Wo9En73Fs18, Wo3En77Fs20), with a few percent graphite and minor amounts of troilite, Ni-Fe metal, and possibly diamond. The rims of olivine grains are reduced (to Fo91) and contain abundant blebs of Fe metal. Silicate mineral grains are equant, anhedral, up to 2 mm across, and lack obvious preferred orientations. Euhedral graphite crystals (to 1 mm x 0.3 mm) are present at silicate grain boundaries, along boundaries and protruding into the silicates, and entirely within silicate mineral grains. Graphite euhedra are also present as radiating clusters and groups of parallel plates grains embedded in olivine; no other ureilite has comparable graphite textures. Minute lumps within graphite grains are possible diamond, inferred to be a result of shock. Other shock effects are limited to undulatory extinction and fracturing. Both ureilites have been weathered significantly. Considering their similar mineralogies, identical mineral compositions, and identical unusual textures, Nova 001 and Nullarbor 010 are probably paired. Based on olivine compositions, Nova 001 and Nullarbor 010 are in Group 1 (FeO-rich) of Berkley et al. (1980). Silicate mineral compositions are consistent with those of others known ureilites. The presence of euhedral graphite crystals within the silicate minerals is consistent with an igneous origin, and suggests that large proportions of silicate magma were present locally and crystallized in situ.

Triman, Allan H.↗

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

Geotail Measurements Compared with the Motions of High-Latitude Auroral Boundaries during Two Substorms

Geotail plasma and field measurements at -95 R(sub E) are compared with extensive ground-based, near-Earth, and geosynchronous measurements to study relationships between auroral activity and magnetotail dynamics during the expansion phases of two substorms. The studied intervals are representative of intermittent, moderate activity. The behavior of the aurora and the observed effects at Geotail for both events are harmonized by the concept of the activation of near-Earth X lines (NEXL) after substorm onsets, with subsequent discharges of one or more plasmoids down the magnetotail. The plasmoids must be viewed as three-dimensional structures which are spatially limited in the dawn-dusk direction. Also, reconnection at the NEXL must proceed at variable rates on closed magnetic field lines for significant times before beginning to reconnect lobe flux. This implies that the plasma sheet in the near-Earth magnetotail is relatively thick in comparison with an embedded current sheet and that both the NEXL and distant X line can be active simultaneously. Until reconnection at the NEXL engages lobe flux, the distant X line maintains control of the poleward auroral boundary. If the NEXL remains active after reaching the lobe, the auroral boundary can move poleward explosively. The dynamics of high-latitude aurora in the midnight region thus provides a means for monitoring these processes and indicating when significant lobe flux reconnects at the NEXL.

Maynard, N. C.↗

A numerical model for coupling between atomization and spray dynamics in liquid rocket thrust chambers

This paper describes a novel method of coupling the atomization and spray combustion processes encountered in coaxial injection elements of liquid rocket engine thrust chambers. This method is based on the Jet-Embedding technique in which the liquid jet core equations and the gas phase equations are solved separately. The liquid and gas phase solutions, however, are coupled through the boundary conditions at the interface between the phases. The computational grid for the gas phase calculations are adapted to the shape of the liquid jet core. The axial variation of droplet sizes are calculated using a stability analysis appropriate for the atomization regime of liquid jet breakup. The predictions of this method have been validated with experimental data on low speed water jets. Using this method, calculations are performed for the SSME fuel preburner single injector flow field. The results obtained are in good agreement with the predictions of the volume-of-fluid method.

Giridharan, M. G.↗

Large-Scale Boundary Layer Ingesting Propulsor Research

NASA’s Advanced Air Transport Technology (AATT) project is investigating boundary layer ingesting (BLI) propulsors for advanced subsonic commercial vehicle concepts to enable the reduction of fuel burn. A multidisciplinary team of researchers from NASA, United Technologies Research Center (UTRC), Virginia Polytechnic University, and the Air Force Arnold Engineering Development Complex developed and tested an embedded BLI inlet and distortion-tolerant fan (BLI2DTF) system in the NASA Glenn Research Center (GRC) 8- foot by 6-foot (8x6) transonic wind tunnel. The test demonstrated the component performance goals necessary for an overall fuel burn reduction of 3 to 5 percent on a large hybrid wing body (HWB) aircraft. Special test equipment, including a raised floor with flow effectors and a bleed system, was developed for use in the 8x6 to produce the appropriate incoming boundary layer representative of an HWB application. Detailed measurements were made to determine the inlet total pressure loss and distortion, fan stage efficiency, and aeromechanic performance including blade vibration stress and displacement response. Results from this test were used as input to a vehicle-level system study performed by the AATT project to assess the impact of BLI on an alternative advanced concept aircraft referred to as the NASA D8 (ND8), which is somewhat similar to the HWB in its integration of the propulsor. This paper will provide an overview of the project timeline, special test equipment needed in the wind tunnel to develop the appropriate incoming boundary layer, and the difficulties in designing a propulsor for the test. The paper will conclude with some representative aerodynamic and aeromechanic data from the test itself and conclude with how this data was used in the ND8 system study.

Aerodynamic performance↗

Large-Scale Boundary Layer Ingesting Propulsor Research

_NASA's Advanced Air Transport Technology (AATT) project is investigating boundary layer ingesting (BLI) propulsors for advanced subsonic commercial vehicle concepts to enable the reduction of fuel burn. A multidisciplinary team of researchers from NASA, United Technologies Research Center (UTRC), Virginia Polytechnic University, and the Air Force Arnold Engineering Development Complex developed and tested an embedded BLI inlet and distortion-tolerant fan (BLI2DTF) system in the NASA Glenn Research Center (GRC) 8-foot by 6-foot (8x6) transonic wind tunnel. The test demonstrated the component performance goals necessary for an overall fuel burn reduction of 3 to 5 percent on a large hybrid wing body (HWB) aircraft. Special test equipment, including a raised floor with flow effectors and a bleed system, was developed for use in the 8x6 to produce the appropriate incoming boundary layer representative of an HWB application. Detailed measurements were made to determine the inlet total pressure loss and distortion, fan stage efficiency, and aeromechanic performance including blade vibration stress and displacement response. Results from this test were used as input to a vehicle-level system study performed by the AATT project to assess the impact of BLI on an alternative advanced concept aircraft referred to as the NASA D8 (ND8), which is somewhat similar to the HWB in its integration of the propulsor. This paper will provide an overview of the project timeline, special test equipment needed in the wind tunnel to develop the appropriate incoming boundary layer, and the difficulties in designing a propulsor for the test. The paper will conclude with some representative aerodynamic and aeromechanic data from the test itself and conclude with how this data was used in the ND8 system study.

Celestina, Mark L.↗

TPSAS-NF1676L-17800-DND

Currently, there are two national challenge problems that guide much of the research in durability and damage tolerance at NASA Langley. The first, Airframe Digital Twin, is a concept that combines as-built vehicle components, as-experienced loads and environments, and other vehicle-specific characteristics to enable ultrahigh fidelity modeling of aircraft and spacecraft throughout their service lives. The second, Materials Genome Initiative, is an analog to the Human Genome Project, and is intended to improve the rate at which materials scientists can discover, understand fundamental physics, and improve material systems. This presentation will highlight several research projects ongoing at NASA Langley that are in support of the above challenge problems. Two of those topics will be the subject of detailed discussion. First, investigations of microstructurally-small fatigue cracking (MSFC) in Al-2Cu and Al-4Cu, fabricated in-house, will be presented. Single- and oligo-crystals of Al-Cu specimens were loaded in uniaxial fatigue, while high-resolution in-situ measurements of deformation were made using image correlation (IC) in a scanning-electron microscope (SEM). The Al-Cu specimens were then replicated as crystal plasticity finite element models (CPFEM), where evolution of slip localization near grain boundaries was computed. Comparison among experiment and CPFEM is made. In addition, XRay diffraction measurements of the as-fabricated specimens were carried out, where direct measurements of the embedded copper precipitates were made, and their influence on growing MSFCs were directly observed. The second main topic will illustrate ongoing work in the area of so-called damage-sensing particles. In this work, shape-memory alloys are embedded in an aluminum alloy matrix. Upon the propagation of a fatigue crack, these particles undergo a strain-induced phase transformation which is detected using an acoustic sensor, providing real-time information on propagating cracks. Experiments and simulations regarding the development of this system will also be detailed.

Jacob Hochhalter↗

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

Distributed Turboelectric Propulsion for Hybrid Wing Body Aircraft

Meeting future goals for aircraft and air traffic system performance will require new airframes with more highly integrated propulsion. Previous studies have evaluated hybrid wing body (HWB) configurations with various numbers of engines and with increasing degrees of propulsion-airframe integration. A recently published configuration with 12 small engines partially embedded in a HWB aircraft, reviewed herein, serves as the airframe baseline for the new concept aircraft that is the subject of this paper. To achieve high cruise efficiency, a high lift-to-drag ratio HWB was adopted as the baseline airframe along with boundary layer ingestion inlets and distributed thrust nozzles to fill in the wakes generated by the vehicle. The distributed powered-lift propulsion concept for the baseline vehicle used a simple, high-lift-capable internally blown flap or jet flap system with a number of small high bypass ratio turbofan engines in the airframe. In that concept, the engine flow path from the inlet to the nozzle is direct and does not involve complicated internal ducts through the airframe to redistribute the engine flow. In addition, partially embedded engines, distributed along the upper surface of the HWB airframe, provide noise reduction through airframe shielding and promote jet flow mixing with the ambient airflow. To improve performance and to reduce noise and environmental impact even further, a drastic change in the propulsion system is proposed in this paper. The new concept adopts the previous baseline cruise-efficient short take-off and landing (CESTOL) airframe but employs a number of superconducting motors to drive the distributed fans rather than using many small conventional engines. The power to drive these electric fans is generated by two remotely located gas-turbine-driven superconducting generators. This arrangement allows many small partially embedded fans while retaining the superior efficiency of large core engines, which are physically separated but connected through electric power lines to the fans. This paper presents a brief description of the earlier CESTOL vehicle concept and the newly proposed electrically driven fan concept vehicle, using the previous CESTOL vehicle as a baseline.

Kim, Hyun Dae↗

Advanced development of the boundary element method for steady-state heat conduction

Considerable progress has been made in recent years toward advancing the state-of-the-art in solid mechanics boundary element technology. In the present work, much of this new technology is applied in the development of a general-purpose boundary element method (BEM) for steady-state heat conduction. In particular, the BEM implementation involves the use of higher-order conforming elements, self-adaptive integration and multi-region capability. Two- and three-dimensional, as well as axisymmetric analysis, are incorporated within a unified framework. In addition, techniques are introduced for the calculation of boundary flux, and for the inclusion of thermal resistance across interfaces. As a final extension, an efficient formulation is developed for the analysis of solid three-dimensional bodies with embedded holes. For this last class of problems, the new BEM formulation is particularly attractive, since use of the alternatives (i.e. finite element or finite difference methods) is not practical. A number of detailed examples illustrate the suitability and robustness of the present approach for steady-state heat conduction.

Dargush, G. F.↗

Direct Reconstruction of Ablative Thermal Protection System Aeroheating Using A Green's Function Approach

A Green’s function inverse heat transfer (IHT) approach is used to reconstruct the surface heating conditions on ablative thermal protection system (TPS) materials from embedded heat flux sensor and temperature probe measurements. The approach models the temperature time-history at the measurement location as a discrete linear system, allowing for the heat flux boundary condition to be recovered directly without the need for time-marching schemes. The effects of material decomposition and pyrolysis gas transport are modeled using an energy source/sink analogue. The performance of the reconstruction approach is analyzed on a 1D test case representative of an atmospheric entry heating scenario. The approach can recover the TPS surface heat flux to within 4% of the input heating condition with a computation time of 2-3 seconds (>3 orders of magnitude faster than current time-marching IHT methods). As a byproduct of the surface heating reconstruction, the algorithm also captures the surface pyrolysis gas mass flux and solid decomposition at multiple through-thickness locations within the TPS.

Kenneth McAfee↗

Mesoscale spiral vortex embedded within a Lake Michigan snow squall band - High resolution satellite observations and numerical model simulations

It is known that Great Lakes snow squall convection occurs in a variety of different modes depending on various factors such as air-water temperature contrast, boundary-layer wind shear, and geostrophic wind direction. An exceptional and often neglected source of data for mesoscale cloud studies is the ultrahigh resolution multispectral data produced by Landsat satellites. On October 19, 1972, a clearly defined spiral vortex was noted in a Landsat-1 image near the southern end of Lake Michigan during an exceptionally early cold air outbreak over a still very warm lake. In a numerical simulation using a three-dimensional Eulerian hydrostatic primitive equation mesoscale model with an initially uniform wind field, a definite analog to the observed vortex was generated. This suggests that intense surface heating can be a principal cause in the development of a low-level mesoscale vortex.

Lyons, Walter A.↗

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗