Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical and computational techniques”

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 631 records · Page 35

Modeling, estimation and identification methods for static shape determination of flexible structures

This paper outlines methods for modeling, identification and estimation for static determination of flexible structures. The shape estimation schemes are based on structural models specified by (possibly interconnected) elliptic partial differential equations. The identification techniques provide approximate knowledge of parameters in elliptic systems. The techniques are based on the method of maximum-likelihood that finds parameter values such that the likelihood functional associated with the system model is maximized. The estimation methods are obtained by means of a function-space approach that seeks to obtain the conditional mean of the state given the data and a white noise characterization of model errors. The solutions are obtained in a batch-processing mode in which all the data is processed simultaneously. After methods for computing the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the related estimation error is conducted. In addition to outlining the above theoretical results, the paper presents typical flexible structure simulations illustrating performance of the shape determination methods.

Rodriguez, G.↗

Advanced Methods for Aircraft Engine Thrust and Noise Benefits: Nozzle-Inlet Flow Analysis

The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. Results obtained are based on analytical methods, numerical simulations and experimental tests at the NASA LaRC and Hampton University computer complexes and experimental facilities. The main objective of this research is injection, mixing and combustion enhancement in propulsion systems. The sub-projects in the reporting period are: (A) Aero-performance and acoustics of Telescope-shaped designs. The work included a pylon set application for SCRAMJET. (B) An analysis of sharp-edged nozzle exit designs for effective fuel injection into the flow stream in air-breathing engines: triangular-round and diamond-round nozzles. (C) Measurement technique improvements for the HU Low Speed Wind Tunnel (HU LSWT) including an automatic data acquisition system and a two component (drag-lift) balance system. In addition, a course in the field of aerodynamics was developed for the teaching and training of HU students.

Morgan, Morris H.↗

Hydration Mechanisms in Nanoparticle Interaction and Surface Energetics

Work on the project advanced molecular level understanding, prediction and control of nanoscale hydration in salt solutions under electric stimuli or ionic patterning. Permeation of nanoporous electrodes at preset voltage underlies the function of ultracapacitors. Transitory regulation of wetting in nanoporous media by electric field spans an array of applications in materials, energy storage, and separation sciences. The physically related modulation of nanoparticle solubility by surface charges can significantly extend the range of the nanomaterial applications, improve processing techniques, and potentially alleviate environmental concerns. Pore/solution equilibria, and activated kinetics of liquid gating by nanoconfined electrolytes are challenging problems at the forefront of experimental and theoretical research. Addressing these problems from a molecular perspective required the development of state-of-the-art simulation algorithms in statistical mechanics to capture complex processes in open systems under electric control. Parallel studies of wetting and dispersibility of polar and ionizing particles aim to uncover predictive relations between electrowetting, chemical functionalization, and geometry of nanomaterial particles. By nonequilibrium dynamic modeling, we elucidated electrolyte flow in nanochanels and associated electrokinetic energy conversion. Research on the program provided training opportunities for the next generation of scientists in computational chemistry. Insights, and methods from the project implicate broad segment of researchers in energy and nanosciences, materials and surface chemistry.

74 ATOMIC AND MOLECULAR PHYSICS↗

Flow field predictions for a slab delta wing at incidence

Theoretical results are presented for the structure of the hypersonic flow field of a blunt slab delta wing at moderately high angle of attack. Special attention is devoted to the interaction between the boundary layer and the inviscid entropy layer. The results are compared with experimental data. The three-dimensional inviscid flow is computed numerically by a marching finite difference method. Attention is concentrated on the windward side of the delta wing, where detailed comparisons are made with the data for shock shape and surface pressure distributions. Surface streamlines are generated, and used in the boundary layer analysis. The three-dimensional laminar boundary layer is computed numerically using a specially-developed technique based on small cross-flow in streamline coordinates. In the rear sections of the wing the boundary layer decreases drastically in the spanwise direction, so that it is still submerged in the entropy layer at the centerline, but surpasses it near the leading edge. Predicted heat transfer distributions are compared with experimental data.

Conti, R. J.↗

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗

Excited State Phenomena in Correlated Nanostructures: Transition Metal Oxide Clusters and Nanocrystals

The primary goal of our research program was to apply advanced computational techniques to predict electronic properties of molecular species. These include transition metal monoxide and dioxide clusters and negatively charged copper oxide clusters which are important systems for their utilization in a variety of technological applications. The enhanced electron-electron interactions inherent in these systems make it challenging to model their excited state properties within a parameter-free approach. The excited state properties in question include the interaction of light with these molecules which results in ejection of electrons from the system, or the absorption of light by the molecular system resulting in a molecule in a higher-energy state. Our overall goal was to investigate the impact of various computational and theoretical approaches on such electronic and optical excitations in these systems. Comparing several variants of theoretical approaches with available experimental data, and high-level quantum chemistry computations helped us identify the level of theory and approximations that are needed to achieve a reasonable balance between accuracy and computational demand.

74 ATOMIC AND MOLECULAR PHYSICS↗

A technique for global monitoring of net solar irradiance at the ocean surface. I - Model

An accurate long-term (84-month) climatology of net surface solar irradiance over the global oceans from Nimbus-7 earth radiation budget (ERB) wide-field-of-view planetary-albedo data is generated via an algorithm based on radiative transfer theory. Net surface solar irradiance is computed as the difference between the top-of-atmosphere incident solar irradiance (known) and the sum of the solar irradiance reflected back to space by the earth-atmosphere system (observed) and the solar irradiance absorbed by atmospheric constituents (modeled). It is shown that the effects of clouds and clear-atmosphere constituents can be decoupled on a monthly time scale, which makes it possible to directly apply the algorithm with monthly averages of ERB planetary-albedo data. Compared theoretically with the algorithm of Gautier et al. (1980), the present algorithm yields higher solar irradiance values in clear and thin cloud conditions and lower values in thick cloud conditions.

Frouin, Robert↗

Constraining chaos: Enforcing dynamical invariants in the training of reservoir computers

Drawing on ergodic theory, we introduce a novel training method for machine learning based forecasting methods for chaotic dynamical systems. The training enforces dynamical invariants—such as the Lyapunov exponent spectrum and the fractal dimension—in the systems of interest, enabling longer and more stable forecasts when operating with limited data. The technique is demonstrated in detail using reservoir computing, a specific kind of recurrent neural network. Finally, results are given for the Lorenz 1996 chaotic dynamical system and a spectral quasi-geostrophic model of the atmosphere, both typical test cases for numerical weather prediction.

97 MATHEMATICS AND COMPUTING↗

Uncertainties in Atomic Data and Their Propagation Through Spectral Models. I.

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data.We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of Oiii and Fe ii and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe ii]. Key words: atomic data - atomic processes - line: formation - methods: data analysis - molecular data - molecular processes - techniques: spectroscopic

linear sets↗

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

Automated matching of pairs of SIR-B images for elevation mapping

During the SIR-B mission in October 1984, a significant number of overlapping synthetic aperture radar (SAR) images of various ground areas was collected. This has offered the first opportunity to perform stereo analyses on images from space that cover large ground areas to determine elevation information. This paper presents the preliminary results of an investigation to obtain elevation data from stereo pairs of SIR-B images. First, the accuracy with which elevation information can be derived from SIR-B image pairs is evaluated theoretically. It is shown that elevation accuracy is a function of the slant range resolution, the incidence angles with which the stereo pair is obtained, the accuracies in spacecraft state estimation, and determination of corresponding pixels in the stereo pair. Next, a hierarchical method is developed to match the corresponding pixels. This method involves iterative removal of local distortions and correlations of pairs of local neighborhoods in the two images. Since it is necessary to perform the matching at every pixel in the image, it is very computationally intensive. Therefore, it has been implemented on the Massively Parallel Processor (MPP) at the Goddard Space Flight Center (GSFC). The MPP's speed permits two iterations of this technique to operate on a pair of 512 x 512 images within 7 s. Results of applying this algorithm of SIR-B images of Mount Shasta, CA, are shown. The matching algorithm performs well in regions of the image with significant features. An approximate elevation image derived from the matching process corresponds to published topographic map data, except for certain obvious discontinuities.

Ramapriyan, H. K.↗

Simulating a pulsed-power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the magnetohydrodynamics (MHD) domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This much more affordable ideal MHD representation is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Fast Time-Stepping Strategy for Dynamical Systems Equipped with a Surrogate Model

Simulation of complex dynamical systems arising in many applications is computationally challenging due to their size and complexity. Model order reduction, machine learning, and other types of surrogate modeling techniques offer cheaper and simpler ways to describe the dynamics of these systems but are inexact and introduce additional approximation errors. In order to overcome the computational difficulties of the full complex models, on one hand, and the limitations of surrogate models, on the other, this work proposes a new accelerated time-stepping strategy that combines information from both. This approach is based on the multirate infinitesimal general-structure additive Runge--Kutta framework. The inexpensive surrogate model is integrated with a small time step to guide the solution trajectory, and the full model is treated with a large time step to occasionally correct for the surrogate model error and ensure convergence. Here, we provide a theoretical error analysis, and several numerical experiments, to show that this approach can be significantly more efficient than using only the full or only the surrogate model for the integration.

Surrogate models↗

Acting to gain information

This report is concerned with agents that act to gain information. In previous work, we developed agent models combining qualitative modeling with real-time control. That work, however, focused primarily on actions that affect physical states of the environment. The current study extends that work by explicitly considering problems of active information-gathering and by exploring specialized aspects of information-gathering in computational perception, learning, and language. In our theoretical investigations, we analyzed agents into their perceptual and action components and identified these with elements of a state-machine model of control. The mathematical properties of each was developed in isolation and interactions were then studied. We considered the complexity dimension and the uncertainty dimension and related these to intelligent-agent design issues. We also explored active information gathering in visual processing. Working within the active vision paradigm, we developed a concept of 'minimal meaningful measurements' suitable for demand-driven vision. We then developed and tested an architecture for ongoing recognition and interpretation of visual information. In the area of information gathering through learning, we explored techniques for coping with combinatorial complexity. We also explored information gathering through explicit linguistic action by considering the nature of conversational rules, coordination, and situated communication behavior.

Rosenchein, Stanley J.↗

Determining the Molecular Packing Arrangements on Protein Crystal Faces by Atomic Force Microscopy

Periodic Bond Chain (PBC) analysis of the packing of tetragonal lysozyme crystals have revealed that there are two possible molecular packing arrangements for the crystal faces. The analysis also predicted that only one of these, involving the formation of helices about the 4(sub 3) axes, would prevail during crystal growth. In this study high resolution atomic force microscopy (AFM) was employed to verify these predictions for the (110) crystal face. A computer program was developed which constructs the expected AFM image for a given tip shape for each possible molecular packing arrangement. By comparing the actual AFM image with the predicted images the correct packing arrangement was determined. The prediction of an arrangement involving 4(sub 3) helices was confirmed in this manner,"while the alternate arrangement was not observed. The investigation also showed the protein molecules were packed slightly closer about the 4(sub 3) axes than in the crystallographic arrangement of the crystal interior. This study demonstrates a new approach for determining the molecular packing arrangements on protein crystal faces. It also shows the power of combining a theoretical PBC analysis with experimental high resolution AFM techniques in probing protein crystal growth processes at the molecular level.

Li, Huayu↗

A Search for Supernova in Star-Burst Galaxies

The rate at which massive supernovae occur in galaxies can be related directly to the rate of formation of massive stars this rate in turn can be computed from measurements and theoretical models of stellar mss distributions. Other measurements such as chemical composition and the dynamical structure of galaxies can also be used to infer the rate of massive supernovae. Various lines of reasoning using these measurements and theories of star formation have predicted a supernova rate which might be as much as ten times the rate of supernovae that had been observed using optical telescopes. The unseen supernovae were thought to be undetected due to a variety of circumstances, with the primary problem being extinction imposed by dust along the line of observation. This research used new infrared observations of the Infrared Space Observatory (ISO) space craft to explore active galaxies with a technique unhindered by dust extinction. New image processing algorithms for ISO point source analysis and cosmic ray deglitching were developed in order to search through a fairly complete set of the ISO data.

Rank, David M.↗

Shear localization in metallic materials at high strain rates

Three factors govern adiabatic shear localization: strain hardening (or softening), strain-rate hardening, and thermal softening. It is typically associated with large shear strains (>1), high strain rates (10 3 –10 7 s –1 ), and high temperatures (40–100% of melting point), all of which happen within narrow regions with widths of about 1–200 μm. It is often an undesirable phenomenon, leading to failure, but there are situations where it is desirable, e. g., the generation of machining chips. Here, we review the development of both theoretical and experimental achievements, from the initiation of shear bands to their propagation with emphasis on three aspects: novel experimental techniques, novel materials, and nano/microstructural effects. The principal characteristics of adiabatic shear bands in metallic materials at the nano- and micro-scale are described. Bands that were formerly identified as transformed actually consist of nanocrystalline/ultrafine grains. These grains result from the breakup of the microstructure by a rotational recrystallization process. The evolution of the microstructure inside shear bands and their interactions for hcp, bcc, and fcc alloys, high-entropy alloys, nanocrystalline alloys, and metallic glasses are analyzed mechanistically. In conclusion, the gaps in the field and opportunities for future research are identified. Modern experimental characterization and computational techniques enable a more profound and predictive understanding of adiabatic shear localization and its avoidance in advanced materials.

36 MATERIALS SCIENCE↗

The leaking mode problem in atmospheric acoustic-gravity wave propagation

The problem of predicting the transient acoustic pressure pulse at long horizontal distances from large explosions in the atmosphere is examined. Account is taken of poles off the real axis and of branch line integrals in the general integral governing the transient waveform. Perturbation techniques are described for the computation of the imaginary ordinate of the poles and numerical studies are described for a model atmosphere terminated by a halfspace with c = 478 m/sec above 125 km. For frequencies less than 0.0125 rad/sec, the GR sub 1 mode, for example, is found to have a frequency dependent amplitude decay of the order of 0.0001 nepers/km. Examples of numerically synthesized transient waveforms are exhibited with and without the inclusion of leaking modes. The inclusion of leaking modes results in waveforms with a more marked beginning rather than a low frequency oscillating precursor of gradually increasing amplitude. Also, the revised computations indicate that waveforms invariably begin with a pressure rise, a result supported by other theoretical considerations and by experimental data.

Kinney, W. A.↗