Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 1,315 records · Page 73

Mesoscale Assimilation of TMI Rainfall Data with 4DVAR: Sensitivity Studies

Sensitivity studies are performed on the assimilation of TRMM (Tropical Rainfall Measurement Mission) Microwave Imager (TMI) derived rainfall data into a mesoscale model using a four-dimensional variational data assimilation (4DVAR) technique. A series of numerical experiments is conducted to evaluate the impact of TMI rainfall data on the numerical simulation of Hurricane Bonnie (1998). The results indicate that rainfall data assimilation is sensitive to the error characteristics of the data and the inclusion of physics in the adjoint and forward models. In addition, assimilating the rainfall data alone is helpful for producing a more realistic eye and rain bands in the hurricane but does not ensure improvements in hurricane intensity forecasts. Further study indicated that it is necessary to incorporate TMI rainfall data together with other types of data such as wind data into the model, in which case the inclusion of the rainfall data further improves the intensity forecast of the hurricane. This implies that proper constraints may be needed for rainfall assimilation.

Tao, Wei-Kuo↗

Electric organ discharges and electric images during electrolocation

Weakly electric fish use active electrolocation - the generation and detection of electric currents - to explore their surroundings. Although electrosensory systems include some of the most extensively understood circuits in the vertebrate central nervous system, relatively little is known quantitatively about how fish electrolocate objects. We believe a prerequisite to understanding electrolocation and its underlying neural substrates is to quantify and visualize the peripheral electrosensory information measured by the electroreceptors. We have therefore focused on reconstructing both the electric organ discharges (EODs) and the electric images resulting from nearby objects and the fish's exploratory behaviors. Here, we review results from a combination of techniques, including field measurements, numerical and semi-analytical simulations, and video imaging of behaviors. EOD maps are presented and interpreted for six gymnotiform species. They reveal diverse electric field patterns that have significant implications for both the electrosensory and electromotor systems. Our simulations generated predictions of the electric images from nearby objects as well as sequences of electric images during exploratory behaviors. These methods are leading to the identification of image features and computational algorithms that could reliably encode electrosensory information and may help guide electrophysiological experiments exploring the neural basis of electrolocation.

Review, Tutorial↗

Realistic Covariance Prediction For the Earth Science Constellations

Routine satellite operations for the Earth Science Constellations (ESC) include collision risk assessment between members of the constellations and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed via Monte Carlo techniques as well as numerically integrating relative probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by NASA Goddard's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the ESC satellites: Aqua, Aura, and Terra

Duncan, Matthew↗

The Impact of Microphysics on Intensity and Structure of Hurricanes

During the past decade, both research and operational numerical weather prediction models, e.g. Weather Research and Forecast (WRF) model, have started using more complex microphysical schemes originally developed for high-resolution cloud resolving models (CRMs) with a 1-2 km or less horizontal resolutions. WFW is a next-generation mesoscale forecast model and assimilation system that has incorporated modern software framework, advanced dynamics, numeric and data assimilation techniques, a multiple moveable nesting capability, and improved physical packages. WFW model can be used for a wide range of applications, from idealized research to operational forecasting, with an emphasis on horizontal grid sizes in the range of 1-10 km. The current WRF includes several different microphysics options such as Lin et al. (1983), WSM 6-class and Thompson microphysics schemes. We have recently implemented three sophisticated cloud microphysics schemes into WRF. The cloud microphysics schemes have been extensively tested and applied for different mesoscale systems in different geographical locations. The performances of these schemes have been compared to those from other WRF microphysics options. We are performing sensitivity tests in using WW to examine the impact of six different cloud microphysical schemes on hurricane track, intensity and rainfall forecast. We are also performing the inline tracer calculation to comprehend the physical processes @e., boundary layer and each quadrant in the boundary layer) related to the development and structure of hurricanes.

Tao, Wei-Kuo↗

The Impact of Microphysics on Intensity and Structure of Hurricanes and Mesoscale Convective Systems

During the past decade, both research and operational numerical weather prediction models, e.g. Weather Research and Forecast (WRF) model, have started using more complex microphysical schemes originally developed for high-resolution cloud resolving models (CRMs) with a 1-2 km or less horizontal resolutions. WRF is a next-generation mesoscale forecast model and assimilation system that has incorporated modern software framework, advanced dynamics, numeric and data assimilation techniques, a multiple moveable nesting capability, and improved physical packages. WRF model can be used for a wide range of applications, from idealized research to operational forecasting, with an emphasis on horizontal grid sizes in the range of 1-10 km. The current WRF includes several different microphysics options such as Purdue Lin et al. (1983), WSM 6-class and Thompson microphysics schemes. We have recently implemented three sophisticated cloud microphysics schemes into WRF. The cloud microphysics schemes have been extensively tested and applied for different mesoscale systems in different geographical locations. The performances of these schemes have been compared to those from other WRF microphysics options. We are performing sensitivity tests in using WRF to examine the impact of six different cloud microphysical schemes on precipitation processes associated hurricanes and mesoscale convective systems developed at different geographic locations [Oklahoma (IHOP), Louisiana (Hurricane Katrina), Canada (C3VP - snow events), Washington (fire storm), India (Monsoon), Taiwan (TiMREX - terrain)]. We will determine the microphysical schemes for good simulated convective systems in these geographic locations. We are also performing the inline tracer calculation to comprehend the physical processes (i.e., boundary layer and each quadrant in the boundary layer) related to the development and structure of hurricanes and mesoscale convective systems.

Tao, Wei-Kuo↗

Explosion/Blast Dynamics for Constellation Launch Vehicles Assessment

An assessment methodology is developed to guide quantitative predictions of adverse physical environments and the subsequent effects on the Ares-1 crew launch vehicle associated with the loss of containment of cryogenic liquid propellants from the upper stage during ascent. Development of the methodology is led by a team at Sandia National Laboratories (SNL) with guidance and support from a number of National Aeronautics and Space Administration (NASA) personnel. The methodology is based on the current Ares-1 design and feasible accident scenarios. These scenarios address containment failure from debris impact or structural response to pressure or blast loading from an external source. Once containment is breached, the envisioned assessment methodology includes predictions for the sequence of physical processes stemming from cryogenic tank failure. The investigative techniques, analysis paths, and numerical simulations that comprise the proposed methodology are summarized and appropriate simulation software is identified in this report.

Baer, Mel↗

Advanced X-Ray Inspection of Reinforced Carbon Composite Materials on the Orbiter Leading Edge Structural Subsystem (LESS)

The post return-to-flight (RTF) inspection methodology for the Orbiter Leading Edge Structural Subsystem (LESS) is currently being defined. Numerous NDT modalities and techniques are being explored to perform the flight-to-flight inspections of the reinforced carbon/carbon (RCC) composite material for impact damage, general loss of mass in the bulk layers, or other anomalous conditions that would pose risk to safe return upon re-entry. It is possible to have an impact upon ascent that is not visually observable on the surface, yet causes internal damage. Radiographic testing may be a useful NDT technique for such occurrences. The authors have performed radiographic tests on full-sized mock samples of LESS hardware with embedded image quality phantoms. Digitized radiographic film, computed radiography and flat panel digital real-time radiography was acquired using a GE Eresco 200 x-ray tube, and Se-75 and Yb-169 radioisotopes.

Hernandez, Jose M.↗

Software Testing and Verification in Climate Model Development

Over the past 30 years most climate models have grown from relatively simple representations of a few atmospheric processes to a complex multi-disciplinary system. Computer infrastructure over that period has gone from punch card mainframes to modem parallel clusters. Model implementations have become complex, brittle, and increasingly difficult to extend and maintain. Existing verification processes for model implementations rely almost exclusively upon some combination of detailed analysis of output from full climate simulations and system-level regression tests. In additional to being quite costly in terms of developer time and computing resources, these testing methodologies are limited in terms of the types of defects that can be detected, isolated and diagnosed. Mitigating these weaknesses of coarse-grained testing with finer-grained "unit" tests has been perceived as cumbersome and counter-productive. In the commercial software sector, recent advances in tools and methodology have led to a renaissance for systematic fine-grained testing. We discuss the availability of analogous tools for scientific software and examine benefits that similar testing methodologies could bring to climate modeling software. We describe the unique challenges faced when testing complex numerical algorithms and suggest techniques to minimize and/or eliminate the difficulties.

Clune, Thomas L.↗

The Impact of Microphysical Schemes on Hurricane Intensity and Track

During the past decade, both research and operational numerical weather prediction models [e.g. the Weather Research and Forecasting Model (WRF)] have started using more complex microphysical schemes originally developed for high-resolution cloud resolving models (CRMs) with 1-2 km or less horizontal resolutions. WRF is a next-generation meso-scale forecast model and assimilation system. It incorporates a modern software framework, advanced dynamics, numerics and data assimilation techniques, a multiple moveable nesting capability, and improved physical packages. WRF can be used for a wide range of applications, from idealized research to operational forecasting, with an emphasis on horizontal grid sizes in the range of 1-10 km. The current WRF includes several different microphysics options. At NASA Goddard, four different cloud microphysics options have been implemented into WRF. The performance of these schemes is compared to those of the other microphysics schemes available in WRF for an Atlantic hurricane case (Katrina). In addition, a brief review of previous modeling studies on the impact of microphysics schemes and processes on the intensity and track of hurricanes is presented and compared against the current Katrina study. In general, all of the studies show that microphysics schemes do not have a major impact on track forecasts but do have more of an effect on the simulated intensity. Also, nearly all of the previous studies found that simulated hurricanes had the strongest deepening or intensification when using only warm rain physics. This is because all of the simulated precipitating hydrometeors are large raindrops that quickly fall out near the eye-wall region, which would hydrostatically produce the lowest pressure. In addition, these studies suggested that intensities become unrealistically strong when evaporative cooling from cloud droplets and melting from ice particles are removed as this results in much weaker downdrafts in the simulated storms. However, there are many differences between the different modeling studies, which are identified and discussed.

Tao, Wei-Kuo↗

The Impact of Microphysical Schemes on Intensity and Track of Hurricane

During the past decade, both research and operational numerical weather prediction models [e.g. Weather Research and Forecasting Model (WRF)] have started using more complex microphysical schemes originally developed for high-resolution cloud resolving models (CRMs) with a 1-2 km or less horizontal resolutions. The WRF is a next-generation meso-scale forecast model and assimilation system that has incorporated a modern software framework, advanced dynamics, numeric and data assimilation techniques, a multiple moveable nesting capability, and improved physical packages. The WRF model can be used for a wide range of applications, from idealized research to operational forecasting, with an emphasis on horizontal grid sizes in the range of 1-10 km. The current WRF includes several different microphysics options. At Goddard, four different cloud microphysics schemes (warm rain only, two-class of ice, two three-class of ice with either graupel or hail) are implemented into the WRF. The performances of these schemes have been compared to those from other WRF microphysics scheme options for an Atlantic hurricane case. In addition, a brief review and comparison on the previous modeling studies on the impact of microphysics schemes and microphysical processes on intensity and track of hurricane will be presented. Generally, almost all modeling studies found that the microphysics schemes did not have major impacts on track forecast, but did have more effect on the intensity. All modeling studies found that the simulated hurricane has rapid deepening and/or intensification for the warm rain-only case. It is because all hydrometeors were very large raindrops, and they fell out quickly at and near the eye-wall region. This would hydrostatically produce the lowest pressure. In addition, these modeling studies suggested that the simulated hurricane becomes unrealistically strong by removing the evaporative cooling of cloud droplets and melting of ice particles. This is due to the much weaker downdraft simulated. However, there are many differences between different modeling studies and these differences were identified and discussed.

Tao, W. K.↗

Comprehensive T-Matrix Reference Database: A 2012 - 2013 Update

The T-matrix method is one of the most versatile, efficient, and accurate theoretical techniques widely used for numerically exact computer calculations of electromagnetic scattering by single and composite particles, discrete random media, and particles imbedded in complex environments. This paper presents the fifth update to the comprehensive database of peer-reviewed T-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2012. It also lists several earlier publications not incorporated in the original database, including Peter Waterman's reports from the 1960s illustrating the history of the T-matrix approach and demonstrating that John Fikioris and Peter Waterman were the true pioneers of the multi-sphere method otherwise known as the generalized Lorenz - Mie theory.

Complex scattering↗

Vibroacoustic Characterization of Corrugated-Core and Honeycomb-Core Sandwich Panels

The vibroacoustic characteristics of two candidate launch vehicle fairing structures, corrugated- core and honeycomb-core sandwich designs, were studied. The study of these structures has been motivated by recent risk reduction efforts focused on mitigating high noise levels within the payload bays of large launch vehicles during launch. The corrugated-core sandwich concept is of particular interest as a dual purpose structure due to its ability to harbor resonant noise control systems without appreciably adding mass or taking up additional volume. Specifically, modal information, wavelength dispersion, and damping were determined from a series of vibrometer measurements and subsequent analysis procedures carried out on two test panels. Numerical and analytical modeling techniques were also used to assess assumed material properties and to further illuminate underlying structural dynamic aspects. Results from the tests and analyses described herein may serve as a reference for additional vibroacoustic studies involving these or similar structures.

Allen, Albert↗

Efficient Trajectory Propagation for Orbit Determination Problems

Regularized formulations of orbital motion apply a series of techniques to improve the numerical integration of the orbit. Despite their advantages and potential applications little attention has been paid to the propagation of the partial derivatives of the corresponding set of elements or coordinates, required in many orbit-determination scenarios and optimization problems. This paper fills this gap by presenting the general procedure for integrating the state-transition matrix of the system together with the nominal trajectory using regularized formulations and different sets of elements. The main difficulty comes from introducing an independent variable different from time, because the solution needs to be synchronized. The correction of the time delay is treated from a generic perspective not focused on any particular formulation. The synchronization using time-elements is also discussed. Numerical examples include strongly-perturbed orbits in the Pluto system, motivated by the recent flyby of the New Horizons spacecraft, together with a geocentric flyby of the NEAR spacecraft.

numerical methods↗

Coilgun Acceleration Model Containing Multiple Interacting Coils

A coilgun operates by pulsing current through an axially-arranged series of independently-controlled coils inductively interacting with a small, electrically-conductive, azimuthally-symmetric projectile to accelerate it to high velocities. The electrical circuits are programmed to pulse current through the coils in such a way so as to impart further electromagnetic acceleration in each stage. A method is developed to calculate the mutual inductance between the coils and between each coil and the projectile. These terms are used to write a system of first-order ordinary differential equations governing the projectile velocity and the current flow in each coil. While the inclusion of the electromagnetic interactions between coils significantly complicates the equation set as more coil sets are included in the problem, casting the problem symbolically in mass matrix form permits solution using standard numerical Runge-Kutta techniques. Comparing a projectile with a single-turn to that comprised of nine-turns, the inductance of the former is much smaller, but this leads to a greater induced projectile current. The lower inductance and greater current appear to offset each other with little difference in the acceleration profile for the two cases. For the limited cases studied, coils with a discharge half-cycle equal to the time for a projectile to transit from one coil to the next yield increased efficiency.

Polzin, Kurt A.↗

X-ray Evaluation of the Marshall Grazing Incidence X-Ray Spectrometer (MaGIXS) Nickel-Replicated Mirrors

X-ray observations of astronomical objects provides diagnostics not available in any other wavelength regime, however the capability of making these observation at a high spatial resolution has proven challenging. Recently, NASA Marshall Space Flight Center (MSFC) has made good progress in employing computer numerical control (CNC) polishing techniques on electroless nickel mandrels as part of our replicated grazing incidence optics program. CNC polishing has afforded the ability to deterministically refine mandrel figure, thereby improving mirror performance. The Marshall Grazing Incidence X-ray Spectrometer (MaGIXS) is a MSFC-led sounding rocket instrument that is designed to make the first ever soft x-ray spectral observations of the Sun spatially resolved along a narrow slit. MaGIXS incorporates some of the first mirrors produced at MSFC using this polishing technique. Here we present the predicted mirror performance obtained from metrology, after completion of CNC polishing, as well as the results of X-ray tests performed on the MaGIXS telescope mirror before and after mounting.

MaGIXS X-ray Tests↗

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in a Mach 4.5 Boundary Layer

We use deep learning, an ensemble variational technique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (DeepONets), which have the known ability to learn complex nonlinear operators within dynamical systems, are used for machine learning. For the baseline configuration of a smooth flat plate, the second-mode waves at the DNS inflow cause a quick nonlinear breakdown of the high-speed boundary layer within the computational domain. Results reported in the present study validate the ability of DeepONets to model the transition delay via a given topography of the roughness element. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substantially lowered by the DeepONets-based reduced-order model. In comparison to the baseline method of EnVar optimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition past the outflow boundary of the computational domain while utilizing almost 5–6 times fewer DNS.

Machine Learning↗

A New Search for Low-mass Dark Matter and an Examination and Reduction of the Uncertainty due to the Photoelectric Absorption Cross Section using a Cryogenic Silicon Detector with Single-charge Sensitivity

Many astrophysical observations point to the abundant existence of dark matter (DM) in the universe composed of particles beyond the Standard Model. However despite numerous experiments using different techniques, DM particles have yet to be directly observed. The Super Cryogenic Dark Matter Search (SuperCDMS) SNOLAB experiment will employ silicon and germanium crystal detectors operated at temperatures as low as 15~mK to probe DM interactions using phonon and ionization signals. Recently, R\&D; facilities have developed gram-sized high-voltage eV-scale (HVeV) silicon detectors that achieve single-electron-hole-pair resolution. By probing effective DM-electron interactions, these devices can be used for searches of low-mass DM candidates.This dissertation presents a DM search experiment known as HVeV Run~2 that employs a second-generation HVeV detector operated in an above-ground laboratory at Northwestern University (IL, USA). Energy spectra are obtained from a blind analysis with 0.39 and 1.2\,g-days of exposure with the detector biased at 60 and 100\,V, respectively. The 0.93~gram detector achieves a 3\,eV phonon energy resolution, corresponding to a world-leading charge resolution of 3\,\% of a single electron-hole pair for a detector bias of 100\,V. With charge carrier trapping and impact ionization effects incorporated into the DM signal models, the resulting exclusion limits are reported for inelastic DM-electron scattering for DM masses from 0.5--$10^{4}$\,MeV$/c^{2}$; in the mass range from 1.2--50\,eV$/c^{2}$ the limits for dark photon and axion-like particle absorption are reported.Several DM search experiments, including HVeV Run~2, are sensitive to low-mass DM candidates that rely on the temperature-dependent photoelectric absorption cross section of silicon. However discrepancies in the underlying literature data result in dominating systematic uncertainties on the DM exclusion limits. In order to reduce these systematic uncertainties, this dissertation presents a novel method of making a direct, low-temperature measurement of the photoelectric absorption cross section of silicon at energies near the band gap (1.2--2.8\,eV).

Wilson, Matthew James↗

Mean Flow from Phase Averages in the 2D Boussinesq Equations

The atmosphere and ocean are described by highly oscillatory PDEs that challenge both our understanding of their dynamics and their numerical approximation. This paper presents a preliminary numerical study of one type of phase averaging applied to mean flows in the 2D Boussinesq equations that also has application to numerical methods. The phase averaging technique, well-known in dynamical systems theory, relies on a mapping using the exponential operator, and then an averaging over the phase. The exponential operator has connections to the Craya–Herring basis pioneered by Jack Herring to study the fluid dynamics of oscillatory, nonlinear fluid dynamics. In this paper, we perform numerical experiments to study the effect of this averaging technique on the time evolution of the solution. We explore its potential as a definition for mean flows. We also show that, as expected from theory, the phase-averaging method can reduce the magnitude of the time rate of change in the PDEs, making them potentially suitable for time stepping methods.

Wingate, Beth A.↗