Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variable Geometry”

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,459 records · Page 81

An Efficient, Multi-Layered Crown Delineation Algorithm for Mapping Individual Tree Structure Across Multiple Ecosystems

Deriving individual tree information from discrete return, small footprint LiDAR data may improve forest above ground biomass estimates, and provide tree-level information that is important in many ecological studies. Several crown delineation algorithms have been developed to extract individual tree information from LiDAR point clouds or rasterized canopy height models (CHM), but many of these algorithms have difficulty discriminating between overlapping crowns, and also may fail to detect understory trees. Our approach uses a watershed based delineation of a CHM, which is subsequently refined using the LiDAR point cloud. Individual tree detection was validated with stem mapped field data from the Smithsonian Environmental Research Center (SERC), Maryland, and on a plot and stand level through comparisons of stem density and basal area to delineated metrics at both SERC and a study area in the Sierra Nevada, California. For individual tree detection, the algorithm correctly identified 70% of dominant trees, 58% of co-dominant trees, 35% of intermediate trees and 21% of suppressed trees at SERC. The algorithm had difficulty distinguishing between crowns of small, dense understory trees of approximately the same height. Delineated crown volume alone explained 53% and 84% of the variability in basal area at the SERC and Sierra Nevada sites, respectively. The algorithm produced crown area distributions comparable to diameter at breast height (DBH) size class distributions observed in the field in both study sites. The algorithm detected understory crowns better in the conifer-dominated Sierra Nevada site than in the closed-canopy deciduous site in Maryland. The ability for the algorithm to reproduce both accurate tree size distributions and individual crown geometries in two dissimilar and complex forests suggests great promise for applicability to a wide range of forest systems.

LiDAR↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine Learning Self-Diffusion Prediction for Lennard-Jones Fluids in Pores

Predicting the diffusion coefficient of fluids under nanoconfinement is important for many applications including the extraction of shale gas from kerogen and product turnover in porous catalysts. Due to the large number of important variables, including pore shape and size, fluid temperature and density, and the fluid–wall interaction strength, simulating diffusion coefficients using molecular dynamics (MD) in a systematic study could prove to be prohibitively expensive. Here, we use machine learning models trained on a subset of MD data to predict the self-diffusion coefficients of Lennard-Jones fluids in pores. Our MD data set contains 2280 simulations of ideal slit pore, cylindrical pore, and hexagonal pore geometries. We use the forward feature selection method to determine the most useful features (i.e., descriptors) for developing an artificial neutral network (ANN) model with an emphasis on easily acquired features. Our model shows good predictive ability with a coefficient of determination (i.e., R 2 ) of ~0.99 and a mean squared error of ~2.9 × 10 –5 . Finally, we propose an alteration to our feature set that will allow the ANN model to be applied to nonideal pore geometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improving River Routing Using a Differentiable Muskingum‐Cunge Model and Physics‐Informed Machine Learning

Recently, rainfall-runoff simulations in small headwater basins have been improved by methodological advances such as deep neural networks (NNs) and hybrid physics-NN models—particularly, a genre called differentiable modeling that intermingles NNs with physics to learn relationships between variables. However, hydrologic routing simulations, necessary for simulating floods in stem rivers downstream of large heterogeneous basins, had not yet benefited from these advances and it was unclear if the routing process could be improved via coupled NNs. We present a novel differentiable routing method (δMC-Juniata-hydroDL2) that mimics the classical Muskingum-Cunge routing model over a river network but embeds an NN to infer parameterizations for Manning's roughness (n) and channel geometries from raw reach-scale attributes like catchment areas and sinuosity. The NN was trained solely on downstream hydrographs. Synthetic experiments show that while the channel geometry parameter was unidentifiable, n can be identified with moderate precision. With real-world data, the trained differentiable routing model produced more accurate long-term routing results for both the training gage and untrained inner gages for larger subbasins (>2,000 km2) than either a machine learning model assuming homogeneity, or simply using the sum of runoff from subbasins. The n parameterization trained on short periods gave high performance in other periods, despite significant errors in runoff inputs. The learned n pattern was consistent with literature expectations, demonstrating the framework's potential for knowledge discovery, but the absolute values can vary depending on training periods. The trained n parameterization can be coupled with traditional models to improve national-scale hydrologic flood simulations.

54 ENVIRONMENTAL SCIENCES↗

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Navier-Stokes analysis of turbine blade heat transfer

Comparisons with experimental heat transfer and surface pressures were made for seven turbine vane and blade geometries using a quasi-three-dimensional thin-layer Navier-Stokes analysis. Comparisons are made for cases with both separated and unseparated flow over a range of Reynolds numbers and freestream turbulence intensities. The analysis used a modified Baldwin-Lomax turbulent eddy viscosity mode. Modifications were made to account for the effects of: (1) freestream turbulence on both transition and leading edge heat transfer; (2) strong favorable pressure gradients on relaminarization; and (3) variable turbulent Prandtl number heat transfer. In addition, the effect of heat transfer on the near wall model of Deissler is compared with the Van Driest model.

Boyle, R. J.↗

Navier-Stokes analysis of turbine blade heat transfer

Comparisons with experimental heat transfer and surface pressures were made for seven turbine vane and blade geometries using a quasi-three-dimensional thin-layer Navier-Stokes analysis. Comparisons are made for cases with both separated and unseparated flow over a range of Reynolds numbers and freestream turbulence intensities. The analysis used a modified Baldwin-Lomax turbulent eddy viscosity mode. Modifications were made to account for the effects of: (1) freestream turbulence on both transition and leading edge heat transfer; (2) strong favorable pressure gradients on relaminarization; and (3) variable turbulent Prandtl number heat transfer. In addition, the effect of heat transfer on the near wall model of Deissler is compared with the Van Driest model.

Boyle, R. J.↗

Multiangle Implementation of Atmospheric Correction (MAIAC):: Radiative Transfer Basis and Look-up Tables - 1

This paper describes a radiative transfer basis of the algorithm MAIAC which performs simultaneous retrievals of atmospheric aerosol and bidirectional surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS). The retrievals are based on an accurate semianalytical solution for the top-of-atmosphere reflectance expressed as an explicit function of three parameters of the Ross-Thick Li-Sparse model of surface bidirectional reflectance. This solution depends on certain functions of atmospheric properties and geometry which are precomputed in the look-up table (LUT). This paper further considers correction of the LUT functions for variations of surface pressure/height and of atmospheric water vapor, which is a common task in the operational remote sensing. It introduces a new analytical method for the water vapor correction of the multiple ]scattering path radiance. It also summarizes the few basic principles that provide a high efficiency and accuracy of the LUT ]based radiative transfer for the aerosol/surface retrievals and optimize the size of LUT. For example, the single-scattering path radiance is calculated analytically for a given surface pressure and atmospheric water vapor. The same is true for the direct surface-reflected radiance, which along with the single-scattering path radiance largely defines the angular dependence of measurements. For these calculations, the aerosol phase functions and kernels of the surface bidirectional reflectance model are precalculated at a high angular resolution. The other radiative transfer functions depend rather smoothly on angles because of multiple scattering and can be calculated at coarser angular resolution to reduce the LUT size. At the same time, this resolution should be high enough to use the nearest neighbor geometry angles to avoid costly three ]dimensional interpolation. The pressure correction is implemented via linear interpolation between two LUTs computed for the standard and reduced pressure levels. A linear mixture and a modified linear mixture methods are used to represent different aerosol types in the aerosol/surface retrievals from several base models of the fine and coarse aerosol fractions. In summary, the developed LUT algorithm allows fast high-accuracy simulations of the outgoing radiance with full variability of the atmospheric and surface bidirectional reflectance properties for the aerosol/surface remote sensing.

Lyapustin, Alexei↗

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗

Numerical analysis of experimental studies of methane hydrate dissociation induced by depressurization in a sandy porous medium

Methane Hydrates (MHs) are a promising energy source abundantly available in nature. Understanding the complex processes of MH formation and dissociation is critical for the development of safe and efficient technologies for energy recovery. Many laboratory and numerical studies have investigated these processes using synthesized MH-bearing sediments. A near-universal issue encountered in these studies is the spatial heterogeneous hydrate distribution in the testing apparatus. In the absence of direct observations (e.g. using X-ray computed tomography) coupled with real time production data, the common assumption made in almost all numerical studies is a homogeneous distribution of the various phases. In an earlier study (Yin et al., 2018) that involved the numerical description of a set of experiments on MH-formation in sandy medium using the excess water method, we showed that spatially heterogeneous phase distribution is inevitable and significant. In the present study, we use as a starting point the results and observations at the end of the MH formation and seek to numerically reproduce the laboratory experiments of depressurization-induced dissociation of the spatially-heterogeneous MH distribution. This numerical study faithfully reproduces the geometry of the laboratory apparatus, the initial and boundary conditions of the system, and the parameters of the dissociation stimulus, capturing accurately all stages of the experimental process. Using inverse modelling (history-matching) that minimized deviations between the experimental observations and numerical predictions, we determined the values of all the important flow, thermal, and kinetic parameters that control the system behaviour, which yielded simulation results that were in excellent agreement with the measurements of key monitored variables, i.e. pressure, temperature, cumulative production of gas and water over time. We determined that at the onset of depressurization (when the pressure drop – the driving force of dissociation – is at its maximum), the rate of MH dissociation approaches that of an equilibrium reaction and is limited by the heat transfer from the system surroundings. Furthermore, as the effect of depressurization declines over time, the dissociation reaction becomes kinetically limited despite significant heat inflows from the boundaries, which lead to localized temperature increases in the reactor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bridging Scales in Bioenergy and Catalysis: A Review of Mesoscale Modeling Applications, Methods, and Future Directions

Between the molecular and reactor scales, which are familiar to the chemical engineering community, lies an intermediate regime, here termed the “mesoscale,” where transport phenomena and reaction kinetics compete on similar time scales. Bioenergy and catalytic processes offer particularly important examples of mesoscale phenomena owing to their multiphase nature and the complex, highly variable porosity characteristic of biomass and many structured catalysts. In this review, we overview applications and methods central to mesoscale modeling as they apply to reaction engineering of biomass conversion and catalytic processing. A brief historical perspective is offered to put recent advances in context. Applications of mesoscale modeling are described, and several specific examples from biomass pyrolysis and catalytic upgrading of bioderived intermediates are highlighted. Methods including reduced order modeling, finite element and finite volume approaches, geometry construction and import, and visualization of simulation results are described; in each category, recent advances, current limitations, and areas for future development are presented. Owing to improved access to high-performance computational resources, advances in algorithm development, and sustained interest in reaction engineering to sustainably meet societal needs, we conclude that a significant upsurge in mesoscale modeling capabilities is on the horizon that will accelerate design, deployment, and optimization of new bioenergy and catalytic technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deep Learning Coordinate-Free Quantum Chemistry

Computing quantum chemical properties of small molecules and polymers can provide insights valuable to physicists, chemists and biologists when designing new materials, catalysts, biological probes and drugs. Deep learning can compute quantum chemical properties accurately in a fraction of the time required by commonly used methods such as density functional theory (DFT). However, many of these deep learning architectures require energy minimized molecular geometries as input, which is also computationally expensive, and decreasing the reproducibility and throughput of these methods. In this study, we demonstrate that accurate quantum chemical computations can be performed without optimized geometries by operating in the coordinate-free domain using deep learning on graph encodings. Furthermore, we also find that the choice of graph-encoding architecture substantially affects the performance of these methods. The Wave architecture outperforms graph convolution architectures, particularly on complex molecules. Furthermore, the structures of these graph encoding architectures provide an opportunity to probe an important, outstanding question in quantum mechanics: What types of quantum chemical properties can be represented by local-variable models? We find that Wave, a local-variable model, is more accurately calculates quantum chemical properties. Graph convolutional architectures require global variables, and are not as effective as as Wave. We anticipate that coordinate-free, deep-learning models of quantum chemistry will become valuable tools in chemistry and biology, enabling researchers to rapidly screen chemical databases or identify new molecules using automated, de-novo design algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Acoustic and aerodynamic performance investigation of inverted velocity profile coannular plug nozzles

The results of model scale parametric static and wind tunnel aerodynamic performance tests on unsuppressed coannular plug nozzle configurations with inverted velocity profile are discussed. The nozzle configurations are high-radius-ratio coannular plug nozzles applicable to dual-stream exhaust systems typical of a variable cycle engine for Advanced Supersonic Transport application. In all, seven acoustic models and eight aerodynamic performance models were tested. The nozzle geometric variables included outer stream radius ratio, inner stream to outer stream ratio, and inner stream plug shape. When compared to a conical nozzle at the same specific thrust, the results of the static acoustic tests with the coannular nozzles showed noise reductions of up to 7 PNdB. Extensive data analysis showed that the overall acoustic results can be well correlated using the mixed stream velocity and the mixed stream density. Results also showed that suppression levels are geometry and flow regulation dependent with the outer stream radius ratio, inner stream-to-outer stream velocity ratio and inner stream velocity ratio and inner stream plug shape, as the primary suppression parameters. In addition, high-radius ratio coannular plug nozzles were found to yield shock associated noise level reductions relative to a conical nozzle. The wind tunnel aerodynamic tests showed that static and simulated flight thrust coefficient at typical takeoff conditions are quite good - up to 0.98 at static conditions and 0.974 at a takeoff Mach number of 0.36. At low inner stream flow conditions significant thrust loss was observed. Using an inner stream conical plug resulted in 1% to 2% higher performance levels than nozzle geometries using a bent inner plug.

Knott, P. R.↗

Geological mapping from spaceborne imaging radars Kentucky-Virginia, USA

Radar images (at wavelength 23.5 cm) of a 50-km-wide swath across Kentucky and Virginia obtained with the Shuttle Imaging Radar experiment (SIR-A) in 1981 and with the Seasat SAR in 1978 are compared. Image tone and texture, lineament mapping, drainage mapping, and the effects of illumination geometry and incidence angle are considered, and sample Landsat images are evaluated. The dominant backscatter effect in the SIR-A images is found to facilitate the mapping of steeply sloping terranes and lineaments shorter than the Seasat length resolution limit of about 15 km. It is determined that optimum enhancement of topographic features is obtained when the radar look angle exceeds the surface slope angle by a discrete amount, avoiding layover or relief displacement. A variable-look-angle radar is needed to maintain low incidence angles in regions with widely varying slope angles, as illustrated by the Landsat MSS images.

Ford, J. P.↗

Spectral behavior of Hercules X-1 - Its long-term variability and pulse phase spectroscopy

Results are presented of the analysis of the broad-energy-band Her X-1 spectra obtained by HEAO 1. Data from 13 pointed observations covered three 35-day cycles in February, August, and September of 1978. The pulse phase-averaged spectrum did not show any apparent variation during the main-on state, except for the overall intensity. The pulse phase-resolved spectrum, on the other hand, showed a moderately strong phase dependence in each of the main-on state pointings. This may imply that the orientation of the emission region is not changing on a monthly time frame, but the geometry of that region may evolve with time.

Soong, Y.↗

Boundary Layer Transition Protuberance Tests at NASA JSC Arc-Jet Facility

A series of arc-jet tests in support of the Shuttle Orbiter Boundary Layer Transition flight experiment was conducted in the Channel Nozzle of the NASA Johnson Space Center Atmospheric Reentry Materials and Structures Facility. The boundary layer trip was a protrusion of a certain height and geometry fabricated as part of a 6"x6" tile insert, a special test article made of the Boeing Rigid Insulation tile material and coated with the Reaction Cured Glass used for the bottom fuselage tiles of the Space Shuttle Orbiter. A total of five such tile inserts were manufactured: four with the 0.25-in. trip height, and one with the 0.35-in. trip height. The tile inserts were interchangeably installed in the center of the 24"x24" variable configuration tile array mounted in the 24"x24" test section of the channel nozzle. The objectives of the test series were to demonstrate that the boundary layer trip can safely withstand the Space Shuttle Orbiter flight-like re-entry environments and provide temperature data on the protrusion surface, surfaces of the nearby tiles upstream and downstream of the trip, as well as the bond line between the tiles and the structure. The targeted test environments were defined for the tip of the protrusion, away from the nominal surface of the tile array. The arc jet test conditions were approximated in order to produce the levels of the free stream total enthalpy at the protrusion height similar to those expected in flight. The test articles were instrumented with surface, sidewall and bond line thermocouples. Additionally, Tempilaq temperature-indicating paint was applied to the nominal tiles of the tile array in locations not interfering with the protrusion trip. Five different grades of paint were used that disintegrate at different temperatures between 1500 and 2000 deg F. The intent of using the paint was to gauge the RCG-coated tile surface temperature, as well as determine its usefulness for a flight experiment. This paper provides an overview of the channel nozzle arc jet, test articles and test conditions, as well as the results of the arc-jet tests including the measured temperature response of the test articles, their pre- and post-test surface scans, condition of the thermal paint, and continents on the protrusion tip heating achieved in tests compared to the computational fluid dynamics predictions.

Larin, M. E.↗

Air traffic management evaluation tool

Method and system for evaluating and implementing air traffic management tools and approaches for managing and avoiding an air traffic incident before the incident occurs. The invention provides flight plan routing and direct routing or wind optimal routing, using great circle navigation and spherical Earth geometry. The invention provides for aircraft dynamics effects, such as wind effects at each altitude, altitude changes, airspeed changes and aircraft turns to provide predictions of aircraft trajectory (and, optionally, aircraft fuel use). A second system provides several aviation applications using the first system. These applications include conflict detection and resolution, miles-in trail or minutes-in-trail aircraft separation, flight arrival management, flight re-routing, weather prediction and analysis and interpolation of weather variables based upon sparse measurements.

Sridhar, Banavar↗

Automated Generation of Finite-Element Meshes for Aircraft Conceptual Design

This paper presents a novel approach for automated generation of fully connected finite-element meshes for all internal structural components and skins of a given wing-body geometry model, controlled by a few conceptual-level structural layout parameters. Internal structural components include spars, ribs, frames, and bulkheads. Structural layout parameters include spar/rib locations in wing chordwise/spanwise direction and frame/bulkhead locations in longitudinal direction. A simple shell thickness optimization problem with two load conditions is used to verify versatility and robustness of the automated meshing process. The automation process is implemented in ModelCenter starting from an OpenVSP geometry and ending with a NASTRAN 200 solution. One subsonic configuration and one supersonic configuration are used for numerical verification. Two different structural layouts are constructed for each configuration and five finite-element meshes of different sizes are generated for each layout. The paper includes various comparisons of solutions of 20 thickness optimization problems, as well as discussions on how the optimal solutions are affected by the stress constraint bound and the initial guess of design variables.

Li, Wu↗