Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “embedded boundary 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.

124 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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↗
NASA NTRS1981

Trailing edge noise prediction from measured surface pressures

Trailing edge (TE) noise is investigated for the case of a two-dimensional airfoil embedded in a uniform low Mach number flow, and the usefulness of several TE noise theories is examined by applying them to the measured data. The TE noise spectra and directivity are quantitatively determined for the case of a high Reynolds number and a fully turbulent boundary layer. Parameters include angle of attack, flow velocity and TE bluntness. Evanescent wave theories by Chase (1975) and Chandiramani (1974) are compared to the present results and show good agreement. Agreement of the near field pressure scatter phenomenon analysis with measurements implies that the basic assumptions used in the analysis are correct, i.e., the turbulent boundary layer (TBL) flow passes the trailing edge into the wake region. No hydrodynamic wake shedding activity is confirmed for the two-sided TBL flow, and a method incorporating the principles of the coherent output power method is used to determine the sound field. The near field edge scatter model is found to not only establish optimum sizing of edge treatment for noise control, but also to separate and identify the scattered field from the incident hydrodynamic field.

Brooks, T. F.↗
NASA NTRS2013

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↗
NASA NTRS2017

Cart3D Simulations for the Second AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all test cases prescribed in the Second AIAA Sonic Boom Prediction Workshop. For each of the four nearfield test cases, we compute pressure signatures at specified distances and off-track angles, using an inviscid, embedded-boundary Cartesian-mesh flow solver with output-based mesh adaptation. The cases range in complexity from an axisymmetric body to a full low-boom aircraft configuration with a powered nacelle. For efficiency, boom carpets are decomposed into sets of independent meshes and computed in parallel. This also facilitates the use of more effective meshing strategies - each off-track angle is computed on a mesh with good azimuthal alignment, higher aspect ratio cells, and more tailored adaptation. The nearfield signatures generally exhibit good convergence with mesh refinement. We introduce a local error estimation procedure to highlight regions of the signatures most sensitive to mesh refinement. Results are also presented for the two propagation test cases, which investigate the effects of atmospheric profiles on ground noise. Propagation is handled with an augmented Burgers' equation method (NASA's sBOOM), and ground noise metrics are computed with LCASB.

Prediction↗
NASA NTRS2018

Design of Distributed Mail-Slot Propulsion System on a Hybrid Wingbody Aircraft

This paper deals with designing a thrust distribution strategy when a Turboelectric Distributed Propulsion (TeDP) system of 16 embedded propulsors is installed on an aerodynamically optimized hybrid wing-body configuration. This HWB previously designed to satisfy conditions of trim, longitudinally static stability and specific cargo space is employed as the baseline configuration for the current study of seeking an optimal propulsion/power system. According to the nature of the entrance flow condition for each distributed propulsion passage in hybrid wing-body aircraft, the ingested boundary layer thickness differs and results in different propulsive reaction. An optimal distribution of thrust and power output is determined by how the system utilizes the propulsive characteristics of each passage. The design space and the number of design variables are selected and described accordingly. An actuator disk model is employed to model thrust generation and shaft power from the propulsor. To carry out the optimization of the propulsion/power system on a computationally expensive CFD model, a Kriging method in conjunction with a Genetic Algorithm (GA) is applied. Throughout the design process, the propulsion performances of the sampled propulsion/power system are analyzed and compared to those of a clean flow engine. The performance metrics includes mass flow rate, fan pressure ratio besides the thrust and shaft power. Minimization of total shaft power from the distributed engine is performed at multiple thrust levels. The benefit of boundary layer ingestion propulsion system is quantified via comparison of thrust equivalent, shaft power and mass flow equivalent clean flow engines with CFD based system design.

Hybrid Wingbody↗
NASA NTRS1992

Crack-mouth displacements for semielliptical surface cracks subjected to remote tension and bending loads

The exact analytical solution for an embedded elliptical crack in an infinite body subjected to arbitrary loading was used in conjunction with the finite element alternating method to obtain crack-mouth-opening displacements (CMOD) for surface cracks in finite plates subjected to remote tension. Identical surface-crack configurations were also analyzed with the finite element method using 20-noded element for plates subjected to both remote tension and bending. The CMODs from these two methods generally agreed within a few percent of each other. Comparisons made with experimental results obtained from surface cracks in welded aluminum alloy specimens subjected to tension also showed good agreement. Empirical equations were developed for CMOD for a wide range of surface-crack shapes and sizes subjected to tension and bending loads. These equations were obtained by modifying the Green-Sneddon exact solution for an elliptical crack in an infinite body to account for finite boundary effects. These equations should be useful in monitoring surface-crack growth in tests and in developing complete crack-face-displacement equations for use in three-dimensional weight-function methods.

Raju, Ivatury S.↗
NASA NTRS2014

Using a Commercial Ethernet PHY Device in a Radiation Environment

This work involved placing a commercial Ethernet PHY on its own power boundary, with limited current supply, and providing detection methods to determine when the device is not operating and when it needs either a reset or power-cycle. The device must be radiation-tested and free of destructive latchup errors. The commercial Ethernet PHY's own power boundary must be supplied by a current-limited power regulator that must have an enable (for power cycling), and its maximum power output must not exceed the PHY's input requirements, thus preventing damage to the device. A regulator with configurable output limits and short-circuit protection (such as the RHFL4913, rad hard positive voltage regulator family) is ideal. This will prevent a catastrophic failure due to radiation (such as a short between the commercial device's power and ground) from taking down the board's main power. Logic provided on the board will detect errors in the PHY. An FPGA (field-programmable gate array) with embedded Ethernet MAC (Media Access Control) will work well. The error detection includes monitoring the PHY's interrupt line, and the status of the Ethernet's switched power. When the PHY is determined to be non-functional, the logic device resets the PHY, which will often clear radiation induced errors. If this doesn't work, the logic device power-cycles the FPGA by toggling the regulator's enable input. This should clear almost all radiation induced errors provided the device is not latched up.

Parks, Jeremy↗
NASA NTRSDate not supplied

Cartesian Mesh Simulations for the Third AIAA Sonic Boom Prediction Workshop

Simulation results are presented for all cases from the Third AIAA Sonic Boom Prediction Workshop. An inviscid, embedded-boundary Cartesian-mesh flow solver is used in conjunction with adjoint-based mesh adaptation to compute nearfield pressure signatures. Specialized techniques are applied to maximize accuracy and minimize cost on Cartesian meshes. The Richardson-based error estimate highlights regions of the signatures most sensitive to mesh refinement. Timing results and coarse, medium, and fine mesh sizes for nearfield cases demonstrate that the parallel decomposition approach is efficient in both computational time and wall-clock. Pressure signals are propagated to the ground using an augmented Burgers’ equation solver to predict boom carpets. Ground signatures and loudness metrics are presented for a standard atmosphere as well as more realistic atmospheric profiles, which affect overall noise levels and can significantly widen the boom carpet. Mesh convergence studies show that high sampling frequencies, around 500 kHz, are required for propagation, and the sampling frequency increases at large off-track angles with longer acoustic ray paths and propagation times. The numerical methods yield accurate results for predicting low sonic boom signatures while being among the least computationally expensive of the workshop.

ARMD↗
NASA NTRS1994

Evolution of dynamo-generated magnetic fields in accretion disks around compact and young stars

Geometrically thin, optically thick, turbulent accretion disks are believed to surround many stars. Some of them are the compact components of close binaries, while the others are throught to be T Tauri stars. These accretion disks must be magnetized objects because the accreted matter, whether it comes from the companion star (binaries) or from a collapsing molecular cloud core (single young stars), carries an embedded magnetic field. In addition, most accretion disks are hot and turbulent, thus meeting the condition for the MHD turbulent dynamo to maintain and amplify any seed field magnetic field. In fact, for a disk's magnetic field to persist long enough in comparison with the disk viscous time it must be contemporaneously regenerated because the characteristic diffusion time of a magnetic field is typically much shorter than a disk's viscous time. This is true for most thin accretion disks. Consequently, studying magentic fields in thin disks is usually synonymous with studying magnetic dynamos, a fact that is not commonly recognized in the literature. Progress in studying the structure of many accretion disks was achieved mainly because most disks can be regarded as two-dimensional flows in which vertical and radial structures are largely decoupled. By analogy, in a thin disk, one may expect that vertical and radial structures of the magnetic field are decoupled because the magnetic field diffuses more rapidly to the vertical boundary of the disk than along the radius. Thus, an asymptotic method, called an adiabatic approximation, can be applied to accretion disk dynamo. We can represent the solution to the dynamo equation in the form B = Q(r)b(r,z), where Q(r) describes the field distribution along the radius, while the field distribution across the disk is included in the vector function b, which parametrically depends on r and is normalized by the condition max (b(z)) = 1. The field distribution across the disk is established rapidly, while the radial distribution Q(r) evolves on a considerably longer timescale. It is this evolution that is the subject of this paper.

Stepinski, Tomasz F.↗
NASA NTRS2004

Design and Evaluation of Modifications to the NASA Langley Flow Impedance Tube

The need to minimize fan noise radiation from commercial aircraft engine nacelles continues to provide an impetus for developing new acoustic liner concepts. If the full value of such concepts is to be attained, an understanding of grazing flow effects is crucial. Because of this need for improved understanding of grazing flow effects, the NASA Langley Research Center Liner Physics Group has invested a large effort over the past decade into the development of a 2-D finite element method that characterizes wave propagation through a lined duct. The original test section in the Langley Grazing IncidenceTube was used to acquire data needed for implementation of this finite element method. This test section employed a stepper motor-driven axial-traversing bar, embedded in the wall opposite the test liner, to position a flush-mounted microphone at pre-selected locations. Complex acoustic pressure data acquired with this traversing microphone were used to educe the acoustic impedance of test liners using this 2-D finite element method and a local optimization technique. Results acquired in this facility have been extensively reported, and were compared with corresponding results from various U.S. aeroacoustics laboratories in the late 1990 s. Impedance data comparisons acquired from this multi-laboratory study suggested that it would be valuable to incorporate more realistic 3-D aeroacoustic effects into the impedance eduction methodology. This paper provides a description of modifications that have been implemented to facilitate studies of 3-D effects. The two key features of the modified test section are (1) the replacement of the traversing bar and its flush-mounted microphone with an array of 95 fixed-location microphones that are flush-mounted in all four walls of the duct, and (2) the inclusion of a suction device to modify the boundary layer upstream of the lined portion of the duct. The initial results achieved with the modified test section are provided in this report, and a comparison of these results with those achieved using the original test section is used to demonstrate that the data acquisition and analysis with the new test section can be confidently used for impedance eduction.

Jones, Michael G.↗
NASA NTRS2013

Diagnosing Aircraft Icing Potential from Satellite Cloud Retrievals

The threat for aircraft icing in clouds is a significant hazard that routinely impacts aviation operations. Accurate diagnoses and forecasts of aircraft icing conditions requires identifying the location and vertical distribution of clouds with super-cooled liquid water (SLW) droplets, as well as the characteristics of the droplet size distribution. Traditional forecasting methods rely on guidance from numerical models and conventional observations, neither of which currently resolve cloud properties adequately on the optimal scales needed for aviation. Satellite imagers provide measurements over large areas with high spatial resolution that can be interpreted to identify the locations and characteristics of clouds, including features associated with adverse weather and storms. This paper describes new techniques for interpreting cloud products derived from satellite data to infer the flight icing threat to aircraft. For unobscured low clouds, the icing threat is determined using empirical relationships developed from correlations between satellite imager retrievals of liquid water path and droplet size with icing conditions reported by pilots (PIREPS). For deep ice over water cloud systems, ice and liquid water content (IWC and LWC) profiles are derived by using the imager cloud properties to constrain climatological information on cloud vertical structure and water phase obtained apriori from radar and lidar observations, and from cloud model analyses. Retrievals of the SLW content embedded within overlapping clouds are mapped to the icing threat using guidance from an airfoil modeling study. Compared to PIREPS and ground-based icing remote sensing datasets, the satellite icing detection and intensity accuracies are approximately 90% and 70%, respectively, and found to be similar for both low level and deep ice over water cloud systems. The satellite-derived icing boundaries capture the reported altitudes over 90% of the time. Satellite analyses corresponding to the time and location of several recent aviation accidents and with icing PIREPS are also presented that reveal skill in identifying severe icing conditions. These results demonstrate the utility of satellite cloud retrievals for quantitatively diagnosing the potential for icing conditions on temporal and spatial scales that should be useful to the aviation community. Plans are being developed to deliver these new satellite products to the GOES-R Proving Ground in the near future so that they can be evaluated in operational applications.

Smith, William L., Jr.↗
NASA NTRS2017

Design of an Axisymmetric Afterbody Test Case for CFD Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flow fields must be supported by continuous improvements in fundamental, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, the underpinnings of a new test platform referred to herein as the NASA Axisymmetric Afterbody are presented. The devised body-of-revolution is a modular platform consisting of a forebody section and afterbody section, allowing for a range of flow behaviors to be studied on interchangeable afterbody geometries. A body-of-revolution offers advantages in shape definition and fabrication, in avoiding direct contact with wind tunnel sidewalls, and in tail-sting integration to facilitate access to higher Reynolds number tunnels. The current work is focused on validation of smooth-body turbulent flow separation, for which a six-parameter body has been developed. A priori RANS computations are reported for a risk-reduction test configuration in order to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. RANS studies of the effects of forebody nose (with/without) and wind tunnel boundary (slip/no-slip) on the selected afterbody are presented. Representative modeling issues that can be explored with this configuration are the effect of higher Reynolds number on separation behavior, flow physics of the progression from attached to increasingly-separated afterbody flows, and the effect of embedded longitudinal vortices on turbulence structure.

Disotell, Kevin J.↗
NASA NTRSDate not supplied

Progress on the Development of a Step Height Sensor for the BOLT Flight Vehicle

This presentation will summarize progress on the development of a step height sensor for the BOLT flight vehicle. Two methods have been proposed to measure the step size in this effort.The first method involves the application of miniaturized camera embedded in the flight vehicle that views the step through a rearward-facing port (shown in Fig. 1a inset). The second method involves the application of a linear fiber optic array bundle to measure the intensity of light reflected back from the aft body step (concept shown in Fig. 1b). Preliminary testing with the miniaturized camera has been performed in the lab and shows that displacements of ≤10 μm can be measured using this method. Final testing with both methods will be performed on the test stand shown in Fig. 1a. This test stand will include a fine adjustment vertical translation stage that will be used to adjust the step height in steps of 1-2 μm. A high-magnification imaging system will also be used to verify the step height. The miniaturized camera and linear fiber optic array bundle will be mounted in one of several port holes with varying degrees of inclination relative to the top surface in the opposing aft-section slab. BOLT flight vehicle details and information on the steps can be found in Ref. [1].

BOLT↗
NASA NTRS2024

Accurate Assessment of Land-Atmosphere Coupling in Climate Models Requires High Frequency Data Output

Land-atmosphere (L-A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centred climate anomalies. Local L-A coupling (LoCo) metrics capture relevant L-A processes, highlighting the impact of soil and vegetation states on surface flux partitioning, and the impact of surface fluxes on boundary layer (BL) growth, development, and entrainment of air above the BL. A primary goal of the Climate Process Team on Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L-A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land-atmosphere interactions span time scales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability of behavioural regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyse, and better understand L-A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behaviour. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant time scales, (ii) enable demonstration of inter-model spread in L-A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

Kirsten L. Findell↗
NASA NTRS2018

Aerodynamic Design and Optimization of Fan Stage for Boundary Layer Ingestion Propulsion System

The present paper addresses the process of preliminary design of a low-pressure fan and outlet guide vane (OGV) of a boundary layer ingestion (BLI) propulsion system. The tail-cone thruster systems of NASA's STARC_ABL (Single-aisle Turboelectric Aircraft with an Aft Boundary-Layer propulsor) adopts an axi-symmetric BLI type inlet as opposed to other embedded engine systems. Thus, the focus of the present work is placed on maximizing the efficiency of the fan and OGV stages under a significant radial distortion. A parameterization with B-spline function for camber line angles, metal chord, thickness distribution and stacking axis of blades is presented. The flowpath lines are also parameterized by B-spline function and aggregated in the design system of blades. The design optimization with evolutionary algorithm is performed with constraints of fan pressure ratio, OGV exit swirl angle and nozzle exit properties. The inlet conditions for the turbo-machinery CFD (Computational Fluid Dynamics) domain and the design goal of the fan stage are driven by a propulsion airframe integration (PAI) model that uses a 3-D unstructured RANS (Reynolds Average Navier Stokes) solver and actuator disk model. The expected power saving of the BLI propulsor is quantified via PAI analysis and the resulting preliminary design of the fan stages is compared with a clean-inlet flow propulsor.

Turbo-machinery↗
NASA NTRS2017

3-D-Observation of Matrix of MIL 090657 Meteorite by Absorption-Phase Tomography

MIL 090657 meteorite (CR2.7) is one of the least altered primitive carbonaceous c hondrites [1]. This meteorite has amorphous silicates like GEMS (glass with embedded metal and sulfide), which are characteristically contained in cometary dust, in matrix [2,3] as with the Paris meteorite [4]. Three lithologies have been recognized; lithology-1 (L 1) dominated by submicron anhydrous silicates, lithology-2 (L2) by GEMS-like amorphous silicates and lithology-3 (L3) by phyllosilicates [2]. Organic materials are abundant in L 1 and L2 [2,3]. L 1 and L2 were further divided into sub-lithology respectively based on their textures and compositions [5]. These studies were performed by 2D SEM and TEM observations of sample surfaces and thin sections that are unable to reveal what constitute each lithology and how these lithologies are distributed and related to each other. This information will provide important insights into alteration and aggregation processes on asteroids and in the early solar nebula. In this study, MIL 090657 matrix was examined in 3D using two types of X-ray tomography; DET (dual-energy tomography) [6] and SIXM (scanning-imaging X-ray microscopy) [7]. Mineral phases can be discriminated based on absorption contrasts at two different X-ray energies in DET. In SIXM, materials composed of light elements such as water or organic materials can be identified based on phase and absorption contrasts. By combining these methods, we can discriminate not only organic materials from voids but also hydrous alteration products, such as hydrated silicates and carbonates, from anhydrous minerals [8]. In this study, we first observed cross sections of MIL 090657 matrix fragments C1 00 mm) in detail using FE-SEM/ EDS. Based on the results, three house-shaped samples (3 0 -50 mm) were extracted from L 1, L2 and their boundary (H1, H3 and H5, respectively) using FIB. 3D imaging of these samples were conducted at BL47XU of SPring-8, a synchrotron radiation facility, with -30-40 nm/ voxel and - 70-80 nm/ voxel at 7keV and 8keV in DET and -100 nm/ voxel at 8keV in SIXM. We found new lithologies that we named L4, L5 and L6 in H1 and H3 in addition to L 1 and L2. L4. L5 and L6 are mainly composed of probably phyllosilicates with different Fe contents. Sulfide and framboidal magnetite were recognized in L4. L5 includes magnetite and carbonate and L6 includes anhydrous silicates having cracks inside. L 1, L2, L4 and L5 are porous while few voids were observed in L6. L4 adjoins to L 1 with boundary, which is not very distinct. L2, L5 and L6 adjoin to each other, and the boundaries of L6 with L2 and L5 are clear. In H5, coarse mineral grains C5-1 0 mm) such as Fe-metal and enstatite are present in L 1 and L2. L 1-L2 boundary is not sharp in 3 D. In conclusion, we found a variety of lithologies by 3D observation for the first time, suggesting that t he MIL 0 90657 meteorite experienced complex alteration and aggregation histories. As L2 is dominated by amorphous silicates, which are extremely susceptible to aqueous alteration, t his is p resumed to be t he most primitive lithology. The contact between L2 and phyllosilicate-bearing lithologies (L5 and L6) with clear boundaries indicates that they were aggregated after aqueous alteration of L5 and L6. The indistinct boundary between L 1 and L2 is suggesting that these two lithologies might originally be the same aggregate composed of amorphous silicates and coarse mineral grains. L 1 might have experienced weak aqueous alteration followed by mild thermal alteration [2], while L2 did not undergo aqueous alteration.

Miyama, Sugimoto↗