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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗

NASA NTRS2023

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Chris Ostrom↗

NASA NTRS1973

A field investigation and numerical simulation of coastal fog

A field investigation of the microphysical and micrometeorological features of fogs occurring near Los Angeles and Vandenberg, California was conducted. Observations of wind speed and direction, temperature, dew point, vertical wind velocity, dew deposition, drop-size distribution, liquid water content, and haze and cloud nucleus concentration were obtained. These observations were initiated in late evening prior to fog formation and continued until the time of dissipation in both advection and radiation fogs. Data were also acquired in one valley fog and several dense haze situations. The behavior of these parameters prior to and during fog are discussed in detail. A two-dimensional numerical model was developed to investigate the formation and dissipation of advection fogs under the influence of horizontal variations in surface temperature. The model predicts the evolution of potential temperature, water vapor content, and liquid water content in a vertical plane as determined by vertical turbulent transfer and horizontal advection. Results are discussed from preliminary numerical experiments on the formation of warm-air advection fog and dissipation by natural and artificial heating from the surface.

Mack, E. J.↗

NASA NTRS2011

Recent Changes in the Greenland Ice Sheet as Seen from Space

Many changes in the Greenland Ice Sheet have been reported in the recent scientific literature and have been attributed to various responses of the ice sheet due to regional (and global) warming. Because melting of the ice sheet would contribute approximately 7 m to sea-level rise, the lives and habitat of hundreds of millions of people worldwide would be directly and indirectly affected if continued ice-sheet melting occurs. As mean-annual global temperatures have increased, there has been an increasing focus on studying the Greenland Ice Sheet using available satellite data, and numerous expeditions have been undertaken. Regional "clear-sky" surface temperature increases since the early 1980s in the Arctic, measured using Advanced Very High Resolution Radiometer (AVHRR) infrared data, range from 0.57+/-0.02 C to 0.72+/-0.10 C per decade. Arctic warming has important implications for ice-sheet mass balance because much of the periphery of the Greenland Ice Sheet is already near O C during the melt season, and is thus vulnerable to more extensive melting if temperatures continue to increase. An increase in melting of the ice sheet would accelerate sea-level rise, an issue of increasing concern to billions of people worldwide. The surface temperature of the ice sheet has been studied in even greater detail using Moderate-Resolution Imaging Spectroradiometer (MODIS) data in the six individual drainage basins as well as for the ice sheet as a whole. Surface temperature trends in the decade of the 2000s have not been strong, according to the MODIS measurements. In addition to surface-temperature increases over the last few decades as measured by AVHRR, other changes have been observed such as accelerated movement of many of Greenland's outlet glaciers and sudden draining of supraglacial lakes. Decreasing mass of the ice sheet since (at least) 2002 has been measured using Gravity Recovery and Climate Experiment (GRACE) data, along with an build-up of ice at the higher elevations and a decrease of ice at the lower elevations as measured using airborne Lidar and Ice, Cloud and Land Elevation Satellite (ICESat) data. The seminar will address the above issues using a variety of NASA satellite data and ground observations.

Hall, Dorothy K.↗

NASA NTRS2000

Numerical Calculation of the Drag Force Acting on an Insoluble Particle Moving in Front of a Solidifying Interface

The distribution of insoluble particles in a metal casting depends primarily on the interaction of the particles with the solid/liquid interface (SLI) during the solidification process. Whether a particle will be engulfed or continuously pushed by SLI is essentially determined by the balance of forces acting on the particle. An important component of this force balance is the drag force generated by the particle motion in front of the SLI. Previously developed mathematical models for particle/SLI interaction made use of steady-state solutions of this force provided by the lubrication theory. However, our numerical model based on the SLI tracking approach shows that not only the steady-state approach is inappropriate to model the interaction process but also that even at steady-state the theoretical solution underestimates the drag force. A regression analysis of steady-state numerical solutions for cylindrical particles moving normal to a flat SLI gave a relationship of the form F(sub d) (exp num) = (square root of 3) pi eta V(sub p)((R(sub p) / d) (exp 10 gamma/3)) as compared to the theoretical solution F(sub d) (exp theor) = 3 X (square root of 2) pi eta V(sub p)((R(sub p) / d) (exp 3/2)) where F(sub d) is the drag force, eta is the dynamic viscosity of the fluid, V(sub p) is the particle velocity, R(sub p) is the particle radius, gamma is Euler's constant (gamma is approximately equal to 0.577), and d is the width of the gap between the particle and the SLI. Model validation as well as the influence of the interface shape on the value of F(sub d) will also be discussed.

Catalina, Adrian V.↗

NASA NTRS1989

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗

NASA NTRS2002

Quadratic Optimization in the Problems of Active Control of Sound

We analyze the problem of suppressing the unwanted component of a time-harmonic acoustic field (noise) on a predetermined region of interest. The suppression is rendered by active means, i.e., by introducing the additional acoustic sources called controls that generate the appropriate anti-sound. Previously, we have obtained general solutions for active controls in both continuous and discrete formulations of the problem. We have also obtained optimal solutions that minimize the overall absolute acoustic source strength of active control sources. These optimal solutions happen to be particular layers of monopoles on the perimeter of the protected region. Mathematically, minimization of acoustic source strength is equivalent to minimization in the sense of L(sub 1). By contrast. in the current paper we formulate and study optimization problems that involve quadratic functions of merit. Specifically, we minimize the L(sub 2) norm of the control sources, and we consider both the unconstrained and constrained minimization. The unconstrained L(sub 2) minimization is certainly the easiest problem to address numerically. On the other hand, the constrained approach allows one to analyze sophisticated geometries. In a special case, we call compare our finite-difference optimal solutions to the continuous optimal solutions obtained previously using a semi-analytic technique. We also show that the optima obtained in the sense of L(sub 2) differ drastically from those obtained in the sense of L(sub 1).

Loncaric, J.↗

NASA NTRS1982

An atmospheric exposure chamber for small animals

The purpose of this project was to design a long-term environmental exposure chamber for small animals. This chamber is capable of producing hypoxic, normoxic and hyperoxic atmospheres which are closely regulated. The chamber, which is of the recycling type, is fashioned after clear plastic germ-free isolators. Oxygen concentration is set and controlled by a paramagnetic O2 analyzer and a 3-way solenoid valve. In this way either O2 or N2 may be provided to the system by way of negative O2 feedback. Relative humidity is maintained at 40-50 percent by a refrigeration type dryer. Carbon dioxide is absorbed by indicating soda lime. A diaphragm pump continuously circulates chamber gas at a high enough flow rate to prevent buildup of CO2 and humidity. This chamber has been used for numerous studies which involve prolonged exposure of small animals to various O2 concentrations.

Glaser, R. M.↗

NASA NTRS1988

Simulation tests of the optimization method of Hopfield and Tank using neural networks

The method proposed by Hopfield and Tank for using the Hopfield neural network with continuous valued neurons to solve the traveling salesman problem is tested by simulation. Several researchers have apparently been unable to successfully repeat the numerical simulation documented by Hopfield and Tank. However, as suggested to the author by Adams, it appears that the reason for those difficulties is that a key parameter value is reported erroneously (by four orders of magnitude) in the original paper. When a reasonable value is used for that parameter, the network performs generally as claimed. Additionally, a new method of using feedback to control the input bias currents to the amplifiers is proposed and successfully tested. This eliminates the need to set the input currents by trial and error.

Paielli, Russell A.↗

NASA NTRS1989

Oscillatory flow with heat transfer in a square cavity

A computational study is presented for the flow inside an oscillatory cavity. The numerical scheme employs a semi-implicit, time-splitting method to integrate the two-dimensional full Navier-Stokes equations satisfying continuity to machine accuracy. The efficient use of direct solvers for the uncoupled momentum and pressure equations is demonstrated. The oscillatory cavity flow is studied considering the effects of heat transfer, Reynolds number, and oscillatory Stokes number.

Biringen, S.↗

NASA NTRS1994

The development and validation of command schedules for SeaWiFS

An automated method for developing and assessing spacecraft and instrument command schedules is presented for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) project. SeaWiFS is to be carried on the polar-orbiting SeaStar satellite in 1995. The primary goal of the SeaWiFS mission is to provide global ocean chlorophyll concentrations every four days by employing onboard recorders and a twice-a-day data downlink schedule. Global Area Coverage (GAC) data with about 4.5 km resolution will be used to produce the global coverage. Higher resolution (1.1 km resolution) Local Area Coverage (LAC) data will also be recorded to calibrate the sensor. In addition, LAC will be continuously transmitted from the satellite and received by High Resolution Picture Transmission (HRPT) stations. The methods used to generate commands for SeaWiFS employ numerous hierarchical checks as a means of maximizing coverage of the Earth's surface and fulfilling the LAC data requirements. The software code is modularized and written in Fortran with constructs to mirror the pre-defined mission rules. The overall method is specifically developed for low orbit Earth-observing satellites with finite onboard recording capabilities and regularly scheduled data downlinks. Two software packages using the Interactive Data Language (IDL) for graphically displaying and verifying the resultant command decisions are presented. Displays can be generated which show portions of the Earth viewed by the sensor and spacecraft sub-orbital locations during onboard calibration activities. An IDL-based interactive method of selecting and testing LAC targets and calibration activities for command generation is also discussed.

Woodward, Robert H.↗

NASA NTRS1994

High-voltage electron collection by a spherical satellite

The numerical solution to the cold electron flow equations is presented for the region surrounding an ionospheric satellite at high positive potential. The momentum, continuity, and Poisson equations are solved simultaneously for satellite potentials of 1000, 10,000, 30,000, and 100,000 kT(kT = 0.17 eV). The electron velocity vectors in the sheath region are presented. The electron current density to the satellite surface is also reported. Electrons impact the satellite surface with average velocity vectors having angles of incidence (relative to the surface normal) that increase from the satellite's magnetic pole to equator. A toroidal flow develops in the vicinity of the magnetic equator and expands away from the satellite surface as the surface potential is increased. The current density to the satellite surface increases with increasing satellite potential and decreases with increasing polar angle.

Sheldon, J. W.↗

NASA NTRS1997

Super-optimal CO2 reduces seed yield but not vegetative growth in wheat

Although terrestrial atmospheric CO2 levels will not reach 1000 micromoles mol-1 (0.1%) for decades, CO2 levels in growth chambers and greenhouses routinely exceed that concentration. CO2 levels in life support systems in space can exceed 10000 micromoles mol-1(1%). Numerous studies have examined CO2 effects up to 1000 micromoles mol-1, but biochemical measurements indicate that the beneficial effects of CO2 can continue beyond this concentration. We studied the effects of near-optimal (approximately 1200 micromoles mol-1) and super-optimal CO2 levels (2400 micromoles mol-1) on yield of two cultivars of hydroponically grown wheat (Triticum aestivum L.) in 12 trials in growth chambers. Increasing CO2 from sub-optimal to near-optimal (350-1200 micromoles mol-1) increased vegetative growth by 25% and seed yield by 15% in both cultivars. Yield increases were primarily the result of an increased number of heads per square meter. Further elevation of CO2 to 2500 micromoles mol-1 reduced seed yield by 22% (P < 0.001) in cv. Veery-10 and by 15% (P < 0.001) in cv. USU-Apogee. Super-optimal CO2 did not decrease the number of heads per square meter, but reduced seeds per head by 10% and mass per seed by 11%. The toxic effect of CO2 was similar over a range of light levels from half to full sunlight. Subsequent trials revealed that super-optimal CO2 during the interval between 2 wk before and after anthesis mimicked the effect of constant super-optimal CO2. Furthermore, near-optimal CO2 during the same interval mimicked the effect of constant near-optimal CO2. Nutrient concentration of leaves and heads was not affected by CO2. These results suggest that super-optimal CO2 inhibits some process that occurs near the time of seed set resulting in decreased seed set, seed mass, and yield.

NASA Discipline Life Support Systems↗

NASA NTRS2005

ON THE ROLE OF INVOLUTIONS IN THE DISCONTINUOUS GALERKIN DISCRETIZATION OF MAXWELL AND MAGNETOHYDRODYNAMIC SYSTEMS

The role of involutions in energy stability of the discontinuous Galerkin (DG) discretization of Maxwell and magnetohydrodynamic (MHD) systems is examined. Important differences are identified in the symmetrization of the Maxwell and MHD systems that impact the construction of energy stable discretizations using the DG method. Specifically, general sufficient conditions to be imposed on the DG numerical flux and approximation space are given so that energy stability is retained These sufficient conditions reveal the favorable energy consequence of imposing continuity in the normal component of the magnetic induction field at interelement boundaries for MHD discretizations. Counterintuitively, this condition is not required for stability of Maxwell discretizations using the discontinuous Galerkin method.

Barth, Timothy↗

NASA NTRS2011

Ground-Based Network and Supersite Observations to Complement and Enrich EOS Research

Since 1997 NASA has been successfully launching a series of satellites - the Earth Observing System (EOS) - to intensively study, and gain a better understanding of, the Earth as an integrated system. Space-borne remote sensing observations, however, are often plagued by contamination of surface signatures. Thus, ground-based in-situ and remote-sensing measurements, where signals come directly from atmospheric constituents, the sun, and/or the Earth-atmosphere interactions, provide additional information content for comparisons that confirm quantitatively the usefulness of the integrated surface, aircraft, and satellite datasets. Through numerous participations, particularly but not limited to the EOS remote-sensing/retrieval and validation projects over the years, NASA/GSFC has developed and continuously refined ground-based networks and mobile observatories that proved to be vital in providing high temporal measurements, which complement and enrich the satellite observations. These are: the AERO NET (AErosol RObotic NETwork) a federation of ground-based globally distributed network of spectral sun-sky photometers; the MPLNET (Micro-Pulse Lidar NETwork, a similarly organized network of micro-pulse lidar systems measuring aerosol and cloud vertical structure continuously; and the SMART-COMMIT (Surface-sensing Measurements for Atmospheric Radiative Transfer - Chemical, Optical & Microphysical Measurements of In-situ Troposphere, mobile observatories, a suite of spectral radiometers and in-situ probes acquiring supersite measurements. Most MPLNET sites are collocated with those of AERONET, and both networks always support the deployment of SMART-COMMIT worldwide. These data products follow the data structure of EOS conventions: Level-0, instrument archived raw data; Level-1 (or 1.5), real-time data with no (or limited) quality assurance; Level-2, not real high temporal and spectral resolutions. In this talk, we will present NASA/GSFC groundbased facilities, serving as network or supersite observations, which have been playing key roles in major international research projects over diverse aerosol regimes to complement and enrich the EOS scientific research.

Tsay, Si-Chee↗

NASA NTRS1991

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗

NASA NTRS2003

Optimized Shapes of Ocsillating Resonators for Generating High-Amplitude Pressure Waves

It is well known that the resonator geometry strongly influences the resonant frequencies of an acoustical resonator and the generated nonlinear standing pressure waveform. Maximizing the ratio of maximum to minimum gas pressure at an end of an oscillating resonator by optimizing the cavity contour is investigated numerically. A quasi-Newton type scheme is used to find optimized axisymmetric resonator shapes to achieve the maximum pressure compression ratio. The acoustical field is solved using a one-dimensional model, and the resonance frequency shift and hysteresis effects are obtained through an automation scheme based on continuation methods. Results are presented from optimizing cone, horn-cone, and cosine resonator geometries. Significant performance improvement is found in the optimized shapes over others previously published. Different optimized shapes are found when starting with different initial guesses, indicating multiple local extrema. The numerical model is validated by comparing with the experimental results of a horn-cone shaped resonator.

Li, Xiao-Fan↗

NASA NTRS1966

Numerical integration of periodic orbits in the elliptic restricted three-body problem

The problem of finding periodic orbits in the circular restricted three-body problem has been very extensively studied in celestial mechanics. It is well known that continuous families of periodic orbits exist, for which the period varies in a continuous way. However, all the applications which are found in the solar system correspond to cases with non-zero eccentricities, and the elliptic restricted three-body problem is thus a better approximation than the circular one. For instance, for the motion of a satellite in the Earth-Moon system, as a first approximation, we may assume that the Moon moves around the Earth in circular motion; but as a much better approximation, we can also assume that the Moon moves in an elliptic orbit around the Earth.

NUMERICAL INTEGRATION↗