Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent discretization”

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 37 records · Page 2

High-Fidelity Multidisciplinary Design Optimization Methodology with Application to Rotor Blades

A multidisciplinary design optimization procedure has been developed and applied to rotorcraft simulations involving tightly-coupled, high-fidelity computational fluid dynamics and comprehensive analysis. A discretely-consistent, adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows on unstructured, dynamic, overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system for high-fidelity rotorcraft analysis is verified; simulation results exhibit good agreement with established solutions. A constrained gradient-based design optimization for a HART-II rotorcraft configuration is demonstrated. The computational cost for individual components of the multidisciplinary sensitivity analysis is assessed and improved.

Wang, Li↗

High-fidelity Multidisciplinary Sensitivity Analysis and Design Optimization for Rotorcraft Applications

A multidisciplinary sensitivity analysis of rotorcraft simulations involving tightly coupled high-fidelity computational fluid dynamics and comprehensive analysis solvers is presented and evaluated. A sensitivity-enabled fluid dynamics solver and a nonlinear flexible multibody dynamics solver are coupled to predict aerodynamic loads and structural responses of helicopter rotor blades. A discretely consistent adjoint-based sensitivity analysis available in the fluid dynamics solver provides sensitivities arising from unsteady turbulent flows and unstructured dynamic overset meshes, while a complex-variable approach is used to compute structural sensitivities with respect to aerodynamic loads. The multidisciplinary sensitivity analysis is conducted through integrating the sensitivity components from each discipline of the coupled system. Accuracy of the coupled system is validated by conducting simulations for a benchmark rotorcraft model and comparing solutions with established analyses and experimental data. Sensitivities of lift computed by the multidisciplinary sensitivity analysis are verified by comparison with the sensitivities obtained by complex-variable simulations. Finally the multidisciplinary sensitivity analysis is applied to a constrained gradient-based design optimization for a HART-II rotorcraft configuration.

Wang, Li↗

High-Fidelity Multidisciplinary Sensitivity Analysis Framework for Multipoint Rotorcraft Optimization

A multidisciplinary, gradient-based sensitivity-analysis methodology is evaluated for optimization of rotorcraft configurations. The tightly coupled discipline models include physics-based fluid dynamics and rotorcraft comprehensive analysis. A discretely-consistent adjoint method accounts for sensitivities of the unsteady flow and unstructured, dynamic, over set grids, while sensitivities of structural responses to aerodynamic loads are computed using a complex-variable method. The methodology is applied to optimize the shape of UH-60A Blackhawk helicopter blades for hover and forward flight conditions. The objective of the multipoint design is to simultaneously increase the rotorcraft figure of merit in a hover flight and reduce the rotor power in a forward flight. Trimmed loose-coupling solutions for the baseline configuration are used to initiate the tight-coupling multidisciplinary analysis. Target thrust and rolling and pitching moments are enforced as optimization constraints. The optimized configuration improves the optimization metrics at both design points. The improved performance and all constraints are maintained over many revolutions beyond the optimization interval, satisfying the required flight conditions. Computational cost of the optimization cycle is assessed in a high-performance computing environment and found affordable for design of rotorcraft in general level-flight conditions.

Wang, Li↗

Toward Adjoint-Based Aeroacoustic Optimization for Propeller and Rotorcraft Applications

The goal of the present project is to build a multidisciplinary, rapid, robust, and accurate computational tool to optimize wing-mounted propeller designs. The full Farassat’s formulation F1Afor aeroacoustic analysis is implemented in the open-source software SU2.This extension enables the prediction of far-field noise generated by moving sources. The formulation is verified, for a stationary and rotating sphere in a wind tunnel and for a tiltrotor in forward flight, by comparing the acoustic predictions ofSU2 with the predictions computed by NASA’s aeroacoustics code ANOPP2. The algorithmic differentiation capability of SU2 provides discretely consistent, adjoint-based sensitivity analysis for this formulation. The adjoint-based sensitivities are verified through comparison with complex-step sensitivities

R Omur Icke↗

On the nature of models in remote sensing

An explicit framework can provide a better understanding of remote sensing models and their interrelationships. This framework distinguishes between the scene, which is real and exists on the ground, and the image, which is a collection of spatially arranged masurements drawn from the scene. The scene model generalizes and parameterizes the essential qualities of the scene. Scene models may be discrete, in which the scene model consists of discrete elements with boundaries, or continuous, in which matter and energy flows are taken to be continuous and there are no clear or sharp boundaries in the scene. In the discrete case, there are two possibilities for models: H- and L-resolution. In the H-resolution case, the resolution cells of the image are smaller than the elements, and thus the elements may be individually resolved. In the L-resolution case, the resolution cells are larger than the elements and cannot be resolved. Most canopy models are L-resolution, deterministic, and noninvertible in nature; image processing models, however, tend to be H-resolution, empirical, and invertible. This taxonomy helps add insight to the development of remote sensing theory and point the way to new, productive areas of research.

Strahler, A. H.↗

Direct and inverse problems in radiation of sound from discrete random sources on two coaxial rings

An analytical model consisting of two ring sources of sound is developed to study the direct radiation in terms of correlation, coherence, and phase and also to aid in solving the inverse-radiation problem of determining the noise source in terms of farfield measurements. The rings consist of discrete sources which are either monopoles or quadrupoles with Gaussian autocorrelations. Only adjacent sources, both within and between the rings, are correlated. Results show that from the farfield information one can determine when the sources are compact or noncompact with respect to the acoustic wavelength and distinguish between the types of sources. In addition, from the inverse-radiation approach one can recover the center of mass, the location and separation distance of the ring, and the respective diameters.

Maestrello, L.↗

A New Discrete Wavelength BUV Algorithm for Consistent Volcanic SO2 Retrievals from Multiple Satellite Missions

This paper describes a new discrete wavelength algorithm developed for retrieving volcanic sulfur dioxide (SO2) vertical column density (VCD) from UV observing satellites. The Multi-Satellite SO2 algorithm (MS_SO2) simultaneously retrieves column densities of sulfur dioxide, ozone, and Lambertian effective reflectivity (LER) and its spectral dependence. It is used operationally to process measurements from the heritage Total Ozone Mapping Spectrometer (TOMS) onboard NASA's Nimbus-7 satellite (N7/TOMS: 1978-1993) and from the current Earth Polychromatic Imaging Camera (EPIC) onboard Deep Space Climate Observatory (DSCOVR: 2015-) from the Earth-Sun Lagrange (L1) orbit. Results from MS_SO2 algorithm for several volcanic cases were assessed using the more sensitive principal component analysis (PCA) algorithm. The PCA is an operational algorithm used by NASA to retrieve SO2 from hyperspectral UV spectrometers, such as the Ozone Monitoring Instrument (OMI) onboard NASA's Earth Observing System Aura satellite and Ozone Mapping and Profiling Suite (OMPS) onboard NASA-NOAA Suomi National Polar Partnership (SNPP) satellite. For this comparative study, the PCA algorithm was modified to use the discrete wavelengths of the Nimbus-7/TOMS instrument, described in Sect. S1 of the Supplement. Our results demonstrate good agreement between the two retrievals for the largest volcanic eruptions of the satellite era, such as the 1991 Pinatubo eruption. To estimate SO2 retrieval systematic uncertainties, we use radiative transfer simulations explicitly accounting for volcanic sulfate and ash aerosols. Our results suggest that the discrete-wavelength MS_SO2 algorithm, although less sensitive than hyperspectral PCA algorithm, can be adapted to retrieve volcanic SO2 VCDs from contemporary hyperspectral UV instruments, such as OMI and OMPS, to create consistent, multi-satellite, long-term volcanic SO2 climate data records.

Bradford L Fisher↗

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

Second moment methods (SMMs) are developed that are consistent with the discontinuous Galerkin spatial discretization of the discrete ordinates (or S\(_N\)) transport equations. The low-order (LO) diffusion system of equations is discretized with fully consistent P\(_1\), local discontinuous Galerkin (LDG), and interior penalty (IP) methods. A discrete residual approach is used to derive SMM correction terms that make each of the LO systems consistent with the high-order discretization. We show that the consistent methods are more accurate and have better solution quality than independently discretized LO systems, that they preserve the diffusion limit, and that the LDG and IP consistent SMMs can be scalably solved in parallel on a challenging, multimaterial benchmark problem.

97 MATHEMATICS AND COMPUTING↗

Discrete modes and continuous spectra in supersonic boundary layers

The disturbance field induced due to a harmonic point source consists of discrete eigenmodes and a continuous spectrum; these are studied by using generalized Fourier transform techniques. For a supersonic boundary layer, there exist seven branches of the continuous spectrum in the complex wavenumber space, four of which (two acoustic waves, one vorticity wave and one entropy wave) contribute to the flowfield downstream of the source. The discrete eigenmodes spring off from these branches at some critical Reynolds numbers. The results for Mach 2 and 4.5 boundary layers show that the receptivity coefficients for the stable discrete modes are much larger than that for the unstable mode. Therefore, the flow very near the source is dominated by the continuous spectrum and the stable discrete modes. However, the unstable mode takes over sufficiently far away from the source. It is shown that it is only necessary to consider the first few discrete modes to construct the solution. Calculations also show that, in a supersonic boundary layer, upstream influence from a localized disturbance is minimal.

Balakumar, P.↗

DEKFIS user's guide: Discrete Extended Kalman Filter/Smoother program for aircraft and rotorcraft data consistency

The computer program DEKFIS (discrete extended Kalman filter/smoother), formulated for aircraft and helicopter state estimation and data consistency, is described. DEKFIS is set up to pre-process raw test data by removing biases, correcting scale factor errors and providing consistency with the aircraft inertial kinematic equations. The program implements an extended Kalman filter/smoother using the Friedland-Duffy formulation.

Source record↗

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The structure of random discrete spacetime

The usual picture of spacetime consists of a continuous manifold, together with a metric of Lorentzian signature which imposes a causal structure on the spacetime. A model, first suggested by Bombelli et al., is considered in which spacetime consists of a discrete set of points taken at random from a manifold, with only the causal structure on this set remaining. This structure constitutes a partially ordered set (or poset). Working from the poset alone, it is shown how to construct a metric on the space which closely approximates the metric on the original spacetime manifold, how to define the effective dimension of the spacetime, and how such quantities may depend on the scale of measurement. Possible desirable features of the model are discussed.

Brightwell, Graham↗

Structure of random discrete spacetime

The usual picture of spacetime consists of a continuous manifold, together with a metric of Lorentzian signature which imposes a causal structure on the spacetime. A model, first suggested by Bombelli et al., is considered in which spacetime consists of a discrete set of points taken at random from a manifold, with only the causal structure on this set remaining. This structure constitutes a partially ordered set (or poset). Working from the poset alone, it is shown how to construct a metric on the space which closely approximates the metric on the original spacetime manifold, how to define the effective dimension of the spacetime, and how such quantities may depend on the scale of measurement. Possible desirable features of the model are discussed.

Brightwell, Graham↗

Image reconstruction from multiple 1-D scans using filtered localized projection

The spatial resolution that can be attained using scanning linear arrays (consisting of discrete IR solid-state detectors) for image acquisition is considered, and a filtered local projection (FLP) method is described which efficiently combines all available scan information into one rectangular grid without the need for explicit interpolation. Mathematically, the FLP algorithm consists of a localized summation followed by an inverse-filter operation, and it has application to nonlinear restoring techniques. The present method is applied, using a linear array, to simulated data for staggered parallel scans and to multiple scan directions. Noise effects and limitations of the technique are also considered.

Frieden, B. Roy↗

Spectroscopic Determination of Trace Contaminants in High-Purity Oxygen

Oxygen used for extravehicular activities (EVAs) must be free of contaminants because a difference in a few tenths of a percent of argon or nitrogen content can mean significant reduction in available EVA time. These inert gases build up in the extravehicular mobility unit because they are not metabolized or scrubbed from the atmosphere. A prototype optical emission technique capable of detecting argon and nitrogen below 0.1% in oxygen has been developed. This instrument uses a glow discharge in reduced-pressure gas to produce atomic emission from the species present. Because the atomic emission lines from oxygen, nitrogen, and argon are discrete, and in many cases well-separated, trace amounts of argon and nitrogen can be detected in the ultraviolet and visible spectrum. This is a straightforward, direct measurement of the target contaminants, and may lend itself to a device capable of on-orbit verification of oxygen purity. A glow discharge is a plasma formed in a low-pressure (1 to 10 Torr) gas cell between two electrodes. Depending on the configuration, voltages ranging from 200 V and above are required to sustain the discharge. In the discharge region, the gas is ionized and a certain population is in the excited state. Light is produced by the transitions from the excited states formed in the plasma to the ground state. The spectrum consists of discrete, narrow emission lines for the atomic species, and broader peaks that may appear as a manifold for molecular species such as O2 and N2, the wavelengths and intensities of which are a characteristic of each atom. The oxygen emission is dominated by two peaks at 777 and 844 nm.

Hornung, Steven↗

Automated propellant leak detection

An automated hydrogen leak detection system is being developed for earth-to-orbit rocket engine applications. The system consists of three elements, a sensor array, a signal processing unit, and a diagnostic processor. The sensor array consists of discrete solid state sensors which are located at specific potential leak sites and in potential leak zones. The signal processing unit provides excitations power for the sensors and provides analog to digital data conversion of the sensor signals. The diagnostic computer analyzes the sensor outputs to determine leak sources and magnitude. Leak data from the sensor network is interpreted using knowledge based software and displayed on-line through a graphical user interface including 3-D leak visualization. The system requirements have been developed assuming eventual application to the Space Shuttle Main Engine which requires approximately 72 measurement locations. A prototype system has been constructed to demonstrate the operational features of the system. This system includes both prototype electronics and data processing software algorithms. Experiments are in progress to evaluate system operation at conditions which simulate prelaunch and flight. The prototype system consists of a network using 16 sensors within a testbed which simulates engine components. Sensor response, orientation, and data analysis algorithms are being evaluated using calibrated leaks produced within the testbed. A prototype flight system is also under development consisting of 8 sensors and flight capable electronics with autonomous control and data recording.

Makel, D. B.↗

Characterizing the microstructures of mammalian enamel by synchrotron phase contrast microCT

The enamel of mammalian teeth is a highly mineralized tissue that must endure a lifetime of cyclic contact and is inspiring the development of next-generation engineering materials. Attempts to implement enamel-inspired structures in synthetic materials have had limited success, largely due to the absence of a detailed understanding of its microstructure. Here, the present work used synchrotron phase-contrast microCT imaging to evaluate the three-dimensional microstructure of enamel from four mammals including Lion, Gray Wolf, Snow Leopard, and Black Bear. Quantitative results of image analysis revealed that the decussation pattern of enamel consists of discrete diazone (D) and parazone (P) bands of rods organized with stacking arrangement of D+/P/D-/P in all mammals evaluated; the D+ and D- refer to distinct diazone bands with juxtaposed rod orientations from the reference plane. Furthermore, the rod orientations in the bands can be described in terms of two principal angles, defined here as the pitch and yaw. While the pitch angle increases from the outer enamel to a maximum (up to ≈ 40°) near the dentin enamel junction, minimal spatial variations are observed in yaw across the enamel thickness. There are clear differences in the decussation parameters of enamel across species that are interpreted here with respect to the structural demands placed on their teeth. The rod pitch and band width of enamel are identified as important design parameters and appear to be correlated with the bite force quotient of the four mammals evaluated.

36 MATERIALS SCIENCE↗