Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗
NASA NTRS2006

Spectral Measurements of PMCs from SBUV/2 Instruments

The SBUV/2 (Solar Backscattered Ultraviolet, model 2) instrument is designed to monitor ozone stratospheric profile and total column ozone using measurements of the Earth's backscattered ultraviolet albedo. We have previously demonstrated that the normal radiance measurements from SBUV/2 instruments, which sample 12 discrete wavelengths between 252 and 340 nm during each scan, can be used to identify polar mesospheric clouds (PMCs). Some SBUV/2 instruments also periodically view the earth in continuous scan mode, covering the wavelength range 160-400 nm with 0.15 nm sampling. Analysis of these data show PMC occurrence rates similar to the normal discrete scan results, although the observation technique reduces the number of daily measurements by a factor of six. PMC observed by SBUV/2 instruments show a monotonic variation in the residual spectral albedo over the wavelength range 250 300 nm, with maximum enhancements of 10 15% at 250 nm. This result is consistent with microphysical model predictions from Jensen [1989. A numerical model of polar mesospheric cloud formation and evolution, Ph. D. Thesis, University of Colorado]. We find no evidence for a systematic localized increase in PMC residual albedo for wavelengths near 260 nm, in contrast to the recently reported results from the MSX UVISI instrument [Carbary J.F., et al., 2004. Evidence for bimodal particle distribution from the spectra of polar mesospheric clouds. Geophysics Research. Letters 31, L13108]. This result is observed for three different SBUV/2 instruments in both Northern and Southern Hemisphere data over a 13-year span. Our Mie scattering calculations show that the location and magnitude of the 260 nm hump feature is dependent upon the specific scattering angles appropriate to the MSX measurements. Although it explains the MSX spectrum, the bimodal size distribution proposed by Carbary et al. (2004), cannot explain the lack of scattering angle dependence of the SBUV/2 spectral shapes. The spectral signature of the SBUV/2 continuous scan PMC data is thus inconsistent with the bimodal particle size distribution suggested by Carbary et al. (2004).

DeLand, Matthew T.↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Attosecond spectroscopy of size-resolved water clusters

Electron dynamics in water are of fundamental importance for a broad range of phenomena, but their real-time study faces numerous conceptual and methodological challenges. Here, in this paper, we introduce attosecond size-resolved cluster spectroscopy and build up a molecular-level understanding of the attosecond electron dynamics in water. We measure the effect that the addition of single water molecules has on the photoionization time delays of water clusters. We find a continuous increase of the delay for clusters containing up to four to five molecules and little change towards larger clusters. We show that these delays are proportional to the spatial extension of the created electron hole, which first increases with cluster size and then partially localizes through the onset of structural disorder that is characteristic of large clusters and bulk liquid water. These results indicate a previously unknown sensitivity of photoionization delays to electron-hole delocalization and indicate a direct link between electronic structure and attosecond photoionization dynamics. Our results offer new perspectives for studying electron-hole delocalization and its attosecond dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Retention of Long-period Gas Giant Planets: Type II Migration Revisited

During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Nonplanar Method for Predicting Incompressible Aerodynamic Coefficients of Rectangular Wings with Circular-Arc Camber

The development of a nonplanar lifting surface method having a continuous distribution of singularities and satisfying the tangent flow boundary condition on the mean camber surface is given. The method predicts some incompressible longitudinal aerodynamic coefficients of rectangular wings which have circular-arc camber. The solution method is of the integral-equation type and the resulting surface integrals are evaluated by either using numerical or analytical techniques, as are appropriate. Applications are made and the results compared with those from an exact two-dimensional circular-arc camber solution, a three-dimensional flat-wing solution which represents the camber by a projected slope onto the flat surface, and a flat-wing experiment. From these comparisons, the present method is found to predict well the flat-wing experiment and limiting values, in addition to the center of pressure variation at an angle of attack of zero for any camber. For wings having camber ratios larger than about 1.25% and moderate to high aspect ratios, the results deterioriate due to the inadequacy of lifting pressure modes employed.

Lamar, J. E.↗
NASA NTRS1973

Transonic flow past lifting wings.

Work conducted by Norstrud (1971) has been extended to lifting flows with the inclusion of embedded continuous supercritical regions. The approach taken follows some fundamental steps proposed for two-dimensional flows by Oswatitsch (1950). The governing integral equation is replaced by a system of nonlinear algebraic equations. The method of parametric differentiation is applied to the solution of this system of equations. The analytical analysis is discussed together with a numerical analysis in the case of a wing configuration with arbitrary thickness distribution.

Norstrud, H.↗
NASA NTRS1982

Numerical simulation of the interaction of the plasma sheet with the lobes of the Earth's magnetotail

Codes involving one and two spatial dimensions and three velocity dimensions were used to model the Earth's magnetotail. It was shown that the magnetotail can become inflated as a consequence of low energy plasma convection toward the neutral plane. The computer study exhibits a conversion of both magnetic field energy and of energy supplied by the convection electric field into particle energy. The numerical simulations suggest that much of the magnetotail substorm morphology may be a simple consequence of an increase, followed by a decrease, in the convection electric field, without the requirement of any magnetospheric size scale plasma instability or other disruptive processes. It is also concluded that the presence of the convection electric field and a continuing replenishment of low energy particles in the magnetotail are both necessary for maintenance of the magnetotail.

Swift, D. W.↗
NASA NTRS1988

The emergence of multi-user expert systems

There exists a set of problems for which the single user workstation is not a viable solution. In some cases, it may be necessary for the system to support numerous experts working on different aspects of a single logical activity concurrently. In others, the activity may extend over a prolonged period of time during which different experts will work on different phases of the project to produce a single integrated result. Another possibility is that the results produced by the expert system must be continuously available to some group other than that actively exercising the expert system. These are the types of applications which require multiuser expert systems. The nature of NASA tends to produce expert system domains which are inherently multiuser. It is these domains that this new class of expert systems has first been encountered. As expert systems move into the main stream of large scale engineering projects, these encounters will become more frequent. A need for expert system building tools which can accommodate the unique properties of these multiuser domains as they emerge is a problem which needs to be addressed as the current set of tools migrate to time sharing environments, and as new tools are developed.

Freeman, Michael S.↗
NASA NTRS1988

Photometry of Pluto-Charon mutual events and Hirayama family asteroids

Once every 124 years, nature provides earth-bound astronomers with the opportunity to observe occultation and transit phenomena between Pluto and its satellite, Charon. Ground-based observations of these events will allow precise physical parameters for the Pluto-Charon system to be derived which are unlikely to be improved upon until in situ spacecraft observations are obtained. The proposed program will continue to support photometry observations from McDonald Observatory, a critical location in an international Pluto Campaign network. Knowledge of the diameters, masses, densities, and compositions derived from these observations will augment our understanding of Pluto's origin and its context within the problem of solar system formation. A second task will continue to research the evolutionary processes which have occurred in the asteroid belt by measuring the physical properties of specific Hirayama family members. Photoelectric lightcurve observations of Koronis and Themis family members will be used to investigate the individual catastrophic collision events which formed each family. By comparing these properties with results of laboratory and numerical experiments, the outcomes of catastrophic disruptions and collisional evolution may be more precisely determined.

Binzel, Richard P.↗
NASA NTRS1990

Investigation of air transportation technology at Princeton University, 1988-1989

The Air Transportation Technology Program at Princeton University, a program emphasizing graduate and undergraduate student research, proceeded along several avenues during the past year. A study of optimal trajectories for penetration of microbursts when encounter is unavoidable was conducted. The emphasis of current wind shear research is on developing an expert system for wind shear avoidance. A knowledge-based reconfigurable flight control system that is implemented with the Pascal programming language using parallel microprocessors was developed. This expert system could be considered a prototype for a failure-tolerant control system that can be constructed using existing hardware. Development of a real-time cockpit simulator continued during the year. The simulator provides a single-person crew station with both conventional and advanced control devices; it currently is programmed to simulate the Navion single-engine general aviation airplane. Alternatives for the air traffic control system giving particular attention to the institutional structure of the FAA are analyzed. A simple numerical procedure for estimating the stochastic robustness of control systems is being investigated. The revitalization of the general aviation industry is also discussed.

Stengel, Robert F.↗
NASA NTRS1993

The determination of crustal motions using laser ranging to Lageos-2 and Lageos-1

During the period October 1992 through June 1993, Center for Space Research efforts continued on four areas: (1) production of a solution for terrestrial reference frame including site positions, velocities, and Earth orientation and rotation spanning the entire Lageos-1 mission, (2) production of a solution for terrestrial reference frame including site positions and Earth orientation spanning the Lageos-2 mission, independently from Lageos-1, (3) evaluation of numerical characteristics of the integration of the Lageos-2 orbit, (4) evaluation of contributions of Lageos-1 and Lageos-2 to temporal variations in the geopotential and 18.6 year tidal response of the Earth, and (5) attendence and participation in Lageos-2 Investigator Working Group Meetings and Committees.

Tapley, Byron D.↗
NASA NTRS1999

Radiation Measurements on the Russian MIR Orbital Station

The Russian Mir orbital station was launched into an orbit with a 51.65 deg inclination in March 1986. It has operated continuously in the altitude range of 380 to 460 km. Almost immediately, cosmonauts began to carry out missions to the station. To date, there have been some 25 Russian missions. As part of the NASA Mir Program, a comprehensive set of radiation measurements were made to map the radiation in all of the Mir module. Numerous radiation measurements have been made on the Mir station throughout its lifetime. However, the comparison of these measurements have been difficult because of different sensitivities of detectors, some active and some passive, differing self shielding and in most cases unknown location shielding. In spite of these complications, very significant progress in the knowledge of the radiation environment onboard the Mir station has been made. These results are directly applicable to expected radiation environment on the International Space Station. In this paper, we describe the combined results from all seven NASA Mir missions. We show: (1) the absorbed dose rate from trapped particles is well correlated with the atmospheric density computed nearly 400 days earlier than the time of observation, (2) developed a relationship between the absorbed dose rate from galactic cosmic rays to the deceleration potential derived using the Climax neutron monitor rate, giving a tool to predict GCR dose rates to +/- 15% nearly 90 days prior to observations, (3) describe the drift of the South Atlantic Anomaly (SAA) with time, (4) compare the predictions of the dose rates as a function of time from the November 6-8, 1998 solar particle event with observations, (5) compare measurements made with NASA Tissue Equivalent Proportional Counter (TEPC) with the ESA DOSTEL device, the Hungarian Pille system, and the Russian R-16 dosimeter. Implications of these measurements for the ISS will be discussed.

Badhwar, Gautam D.↗
NASA NTRS2001

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗
NASA NTRS2001

A Hybrid Numerical Method for Turbulent Mixing Layers

A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental data.

Georgiadis, Nicholas J.↗
NASA NTRS2012

The Advanced Technology Microwave Sounder (ATMS): A New Operational Sensor Series

ATMS is a new satellite microwave sounding sensor designed to provide operational weather agencies with atmospheric temperature and moisture profile information for global weather forecasting and climate applications. ATMS will continue the microwave sounding capabilities first provided by its predecessors, the Microwave Sounding Unit (MSU) and Advanced Microwave Sounding Unit (AMSU). The first ATMS was launched October 28, 2011 on board the Suomi National Polar-orbiting Partnership (S-NPP) satellite. Microwave soundings by themselves are the highest-impact input data used by Numerical Weather Prediction (NWP) models; and ATMS, when combined with the Cross-track Infrared Sounder (CrIS), forms the Cross-track Infrared and Microwave Sounding Suite (CrIMSS). The microwave soundings help meet NWP sounding requirements under cloudy sky conditions and provide key profile information near the surface

Kim, Edward↗
NASA NTRS2021

Assessment of the Advanced Very High Resolution Radiometer (AVHRR) for Snowfall Retrieval in High Latitudes Using CloudSat and Machine Learning

Precipitation retrieval is a challenging topic, especially in high latitudes (HL), and current precipitation products face ample challenges over these regions. This study investigates the potential of the Advanced Very High Resolution Radiometer (AVHRR) for snowfall retrieval in HL using CloudSat radar information and machine learning (ML). With all the known limitations, AVHRR observations should be considered for HL snowfall retrieval because 1) AVHRR data have been continuously collected for about four decades on multiple platforms with global coverage, and similar observations will likely continue in the future; 2) current passive microwave satellite precipitation products have several issues over snow and ice surfaces; and 3) good coincident observations between AVHRR and CloudSat are available for training ML algorithms. Using ML, snowfall rate was retrieved from AVHRR’s brightness temperature and cloud probability, as well as auxiliary information provided by numerical reanalysis. The results indicate that the ML-based retrieval algorithm is capable of detection and estimation of snowfall with comparable or better statistical scores than those obtained from the Atmospheric Infrared Sounder (AIRS) and two passive microwave sensors contributing to the Global Precipitation Measurement (GPM) mission constellation. The outcomes also suggest that AVHRR-based snowfall retrievals are spatially and temporally reasonable and can be considered as a quantitatively useful input to the merged precipitation products that require frequent sampling or long-term records.

Mohammad Reza↗
NASA NTRS1990

Numerical modeling of the atmosphere with an isentropic vertical coordinate

A theta-coordinate model simulating the nonlinear evolution of a baroclinic wave is presented. In the model, vertical discretization maintains important integral constraints such as conservation of the angular momentum and total energy. A massless-layer approach is used in the treatment of the intersections of coordinate surfaces with the lower boundary. This formally eliminates the intersection problem, but raises other computational problems. Horizontal discretization of the continuity and momentum equations in the model are designed to overcome these problems. Selected results from a 10-day integration with the 25-layer, beta-plane version of the model are presented. It is concluded that the model can simulate the nonlinear evolution of a baroclinic wave and associated dynamical processes without major computational difficulties.

Hsu, Yueh-Jiuan G.↗