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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Application of a flux-split algorithm to chemically relaxing, hypervelocity blunt-body flows

Viscous, nonequilibrium, hypervelocity flow fields over two axisymmetric configurations are numerically simulated using a factored, implicit, flux-split algorithm. The governing gas-dynamic and species-continuity equations for laminar flow are presented. The gas-dynamics/nonequilibrium-chemistry coupling procedure is developed as part of the solution procedure and is described in detail. Numerical solutions are presented for hypervelocity flows over a hemisphere and over an axisymmetric aeroassisted orbital transfer vehicle using three different chemistry models. The gas models considered are those for an ideal gas, for a frozen gas, and for chemically relaxing air consisting of five species. The calculated results are compared with existing numerical solutions in the literature along the stagnation line of the hemisphere. The effects of free-stream Reynolds number on the nonequilibrium flow field are discussed.

Balakrishnan, A.↗
DOE OSTI2022

Ensuring planetary survival: the centrality of organic carbon in balancing the multifunctional nature of soils

Not only do soils provide 98.7% of the calories consumed by humans, they also provide numerous other functions upon which planetary survivability closely depends. However, our continuously increasing focus on soils for biomass provision (food, fiber, and energy) through intensive agriculture is rapidly degrading soils and diminishing their capacity to deliver other vital functions. These tradeoffs in soil functionality – the increased provision of one function at the expense of other critical planetary functions – are the focus of this review. Here, we examine how land-use change for biomass provision has decreased the ability of soils to regulate the carbon pool and thereby contribute profoundly to climate change, to cycle the nutrients that sustain plant growth and ecosystem health, to protect the soil biodiversity upon which many other functions depend, and to cycle the Earth’s freshwater supplies. We also examine how this decreasing ability of soil to provide these other functions can be halted and reversed. Despite the complexity and the interconnectedness of soil functions, we show that soil organic carbon plays a central role and is a master indicator for soil functioning and that we require a better understanding of the factors controlling the behavior and persistence of C in soils. Given the threats facing humanity and their economies, it is imperative that we recognize that Soil Security is itself an existential challenge and that we need to increase our focus on the multiple functions of soils for long-term human welfare and survivability of the planet.

Carboncarbon sequestration↗
DOE OSTI2020

The look-elsewhere effect from a unified Bayesian and frequentist perspective

When searching over a large parameter space for anomalies such as events, peaks, objects, or particles, there is a large probability that spurious signals with seemingly high signi ficance will be found. This is known as the look-elsewhere effect and is prevalent throughout cosmology, (astro)particle physics, and beyond. To avoid making false claims of detection, one must account for this effect when assigning the statistical significance of an anomaly. This is typically accomplished by considering the trials factor, which is generally computed numerically via potentially expensive simulations. In this paper we develop a continuous generalization of the Bonferroni and Sidak corrections by applying the Laplace approximation to evaluate the Bayes factor, and in turn relating the trials factor to the prior-to-posterior volume ratio. Here, we use this to define a test statistic whose frequentist properties have a simple interpretation in terms of the global p-value, or statistical significance. We apply this method to various physics-based examples and show it to work well for the full range of p-values, i.e. in both the asymptotic and non-asymptotic regimes. We also show that this method naturally accounts for other model complexities such as additional degrees of freedom, generalizing Wilks' theorem. This provides a fast way to quantify statistical significance in light of the look-elsewhere effect, without resorting to expensive simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

NUM-DAT File Format Specification: Used in M-9 Gun Experiment Data Archiving

The M-9 Shock and Detonation Physics group executes experiments on gun and explosive platforms with large numbers of oscilloscopes used for data acquisition. The data acquisition from these oscilloscopes was automated many years ago using a custom piece of software called RunDig . The default save format from this software is a custom structure referred to as "NUM-DAT" format. This file format includes a text ".DAT" file which is a header file used to interpret the binary ".NUM" file which contains the oscilloscope data. The data save format was originally developed by John Vorthman and has been in use by M-9 personnel for over 20 years. This data format has been used for archiving data from experiments performed by M-9 personnel at TA-40, TA-39, and the TA-55 Impact Test Facility. Numerous custom analysis and visualization programs have also been developed, and continue to be used, that utilize this data format. This document describes the NUM-DAT format and provides code examples for reading the format and converting it to other formats.

47 OTHER INSTRUMENTATION↗
NASA NTRS1961

Theory of Electrostatic Fields in the Ionosphere at Equatorial Latitudes

The properties of the elongated electrostatic fields that are required to provide the coupling mechanism in the dynamo-motor concept of the E and F regions of the ionosphere are examined theoretically for the conditions that prevail in equatorial latitudes. The analysis is developed for an electrostatic field of arbitrary horizontal scale in a horizontally stratified partly ionized gas subject to an imposed magnetic field having the form of a parabolic arch over the equator. The anisotropic character and continuous variation with height of the conductivity are retained throughout , and numerical solutions are determined for the attenuation of the electric field with distance along the field line. The results are similar qualitatively to those found previously upon analysis of the corresponding problem for middle latitudes, but the attenuation of the electrostatic field with height is considerably greater. It develops, in particular, that the coupling between E and F regions is very small for fields having horizontal wavelengths of a few kilometers. It follows that the dynamo-motor concept could not be used to account for the presence of irregularities of this scale. On the other hand, it is found that almost all the attenuation occurs at heights near that of the E region. If an electrostatic field having a horizontal wavelength of a few kilometers could be produced at heights of, say 200 km or greater, it follows that it would be very elongated and extend with little change in amplitude from hemisphere to hemisphere. The results also indicate that effective coupling could be achieved for fields having horizontal wavelengths of the order of tens of kilometers, but the assumptions introduced to simplify the analysis may impair the quantitative reliability of the results for fields of this scale.

Briggs, Benjamin R.↗
NASA NTRS1982

Computations and turbulent flow modeling in support of helicopter rotor technology

Efforts continue to investigate the applicability of the Deiwert's time dependent numerical airfoil calculation code to the simulation of two dimensional airfoil flows with large amounts of separation present. A modified logarithmic region mixing length term produced substantial differences in the calculated flow fields for the large separation of the M=0.4, high angle of attack cases. Reasons are given for interest in applying that model to previously investigated cases of the 64A010 airfoil section at Mach number 0.8 for angles of attack of 4 deg and 6.2 deg. The effects of these turbulence models at alterations and the time dependent behavior from impulsively started conditions are discussed.

Rose, W. C.↗
NASA NTRS1989

Aeroassist flight experiment aerodynamics and aerothermodynamics

The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.

Brewer, Edwin B.↗
NASA NTRS1989

Low frequency electric and magnetic fields

Following preliminary investigations of the low frequency electric and magnetic fields that may exists in the Earth-ionospheric cavity, measurements were taken with state-of-the art spectrum analyzers. As a follow up to this activity, an investigation was initiated to determine sources and values for possible low frequency signal that would appear in the cavity. The lowest cavity resonance is estimated at about 8 Hz, but lower frequencies may be an important component of our electromagnetic environment. The potential field frequencies produced by the electron were investigated by a classical model that included possible cross coupling of the electric and gravitation fields. During this work, an interesting relationship was found that related the high frequency charge field with the extremely low frequency of the gravitation field. The results of numerical calculations were surprisingly accurate and this area of investigation is continuing. The work toward continued development of a standardized monitoring facility is continuing with the potential of installing the prototype at West Virginia State College early in 1990. This installation would be capable of real time monitoring of ELF signals in the Earth-ionoshpere cavity and would provide some directional information. A high gain, low noise, 1/f frequency corrected preamplifier was designed and tested for the ferrite core magnetic sensor. The potential application of a super conducting sensor for the ELF magnetic field detection is under investigation. It is hoped that a fully operational monitoring network could pinpoint the location of ELF signal sources and provide new information on where these signals originate and what causes them, assuming that they are natural in origin.

Spaniol, Craig↗
NASA NTRS1991

Numerical methods for the analysis of sampled-data systems and for the computation of system zeros

MARSYAS is a computer aided control system design package for the simulation and analysis of open loop and closed loop dynamic systems. Outlined here is the numerical and theoretical basis of the MARSYAS functions developed during the summer of 1991. In particular, the numerical computation of the matrix exponential e(exp A) = I + A + A(exp 2)/2! + A(exp 3)/3! + ... and the numerical computation of the finite system zeros of a dynamic system (continuous time or discrete time) are presented.

Hodel, A. Scottedward↗
NASA NTRS1992

Lidar Atmospheric Observatory in the Canadian Arctic

During the last decade there have been growing concerns about a broad variety of atmospheric properties. Among these, a depletion of the stratospheric ozone layer has attracted considerable attention from the general public, politicians and scientists due to its vital impact for the entire global biosphere. One of the major warning signs was the discovery of the 'ozone hole' in the Antarctic region where the concentration of the ozone in the stratosphere was significantly reduced. At present the stratospheric ozone layer in this region is being continuously monitored by groups of scientists from around the world and numerous observations of the ozone layer on the global scale have clearly demonstrated the process of ozone depletion. Recent observations by NASA have shown significant depletion in the Arctic region. This paper provides an initial description of two lidars that are planned to be installed in a new observatory for atmospheric studies in the Canadian Arctic. This observatory is being constructed under the supervision of the Atmospheric Environment Services (AES) of Canada as a part of Green Plan - an initiative of the Federal Government of Canada. The station is located at Eureka on Ellesmere Island at a latitude of 80 degrees N and a longitude of 86 degrees W.

Ulitsky, Arkady↗
NASA NTRS1992

A boundary integral method for an inverse problem in thermal imaging

An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.

Bryan, Kurt↗
NASA NTRS1987

Needs and status of CFD code validation

The two types of Computational Fluid Dynamics code validations, solution-to-solution comparison and solution-to-experiment comparison, are discussed. It is suggested that to develop more detailed experiments the following things are necessary: (1) further development of turbulence models; (2) better methods for numerical validation of CFD codes; (3) evaluation of disagreements; and (4) continued determination of experimental scatter. All data and results are presented in viewgraph form.

Holst, Terry L.↗
NASA NTRS1994

Structure and Dynamics of the Solar Corona

Advanced computational techniques were used to study solar coronal heating and coronal mass ejections. A three dimensional, time dependent resistive magnetohydrodynamic code was used to study the dynamic response of a model corona to continuous, slow, random magnetic footpoint displacements in the photosphere. Three dimensional numerical simulations of the response of the corona to simple smooth braiding flows in the photosphere were calculated to illustrate and understand the spontaneous formation of current filaments. Two dimensional steady state helmet streamer configurations were obtained by determining the time asymptotic state of the interaction of an initially one dimensinal transponic solar wind with a spherical potential dipole field. The disruption of the steady state helmet streamer configuration was studied as a response to shearing of the magnetic footpoints of the closed field lines under the helmet.

Schnack, D. D.↗
NASA NTRS1994

The birthplace of compact groups of galaxies

We use complete redshift surveys to study the redshift neighborhoods of 38 Hickson compact groups (HCGs). Twenty-nine of these HCGs (76%) are embedded in rich looser systems which we call HCG associations. Analysis of the redshift neighborhood of HCGs outside the CfA survey suggests that most HCGs are embedded in more extended physical systems. Rich loose groups extracted from the CfA survey (Ramella et al. (1994)) have physical properties similar to those of the HCG associations. These rich loose groups often contain compact configurations. N-body experiments (Diaferio (1994)) suggest that compact configurations analogous to HCGs form continually during the collapse of rich loose groups. These observational and numerical results suggest that rich loose groups are the birthplace of HCGs.

Ramella, Massimo↗
NASA NTRS1994

Parallel inhomogeneity and the Alfven resonance. 1: Open field lines

In light of a recent demonstration of the general nonexistence of a singularity at the Alfven resonance in cold, ideal, linearized magnetohydrodynamics, we examine the effect of a small density gradient parallel to uniform, open ambient magnetic field lines. To lowest order, energy deposition is quantitatively unaffected but occurs continuously over a thickened layer. This effect is illustrated in a numerical analysis of a plasma sheet boundary layer model with perfectly absorbing boundary conditions. Consequences of the results are discussed, both for the open field line approximation and for the ensuing closed field line analysis.

Hansen, P. J.↗
NASA NTRS1997

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel↗
NASA NTRS2006

Metrics for Linear Kinematic Features in Sea Ice

The treatment of leads as cracks or discontinuities (see Coon et al. presentation) requires some shift in the procedure of evaluation and comparison of lead-resolving models and their validation against observations. Common metrics used to evaluate ice model skills are by and large an adaptation of a least square "metric" adopted from operational numerical weather prediction data assimilation systems and are most appropriate for continuous fields and Eilerian systems where the observations and predictions are commensurate. However, this class of metrics suffers from some flaws in areas of sharp gradients and discontinuities (e.g., leads) and when Lagrangian treatments are more natural. After a brief review of these metrics and their performance in areas of sharp gradients, we present two new metrics specifically designed to measure model accuracy in representing linear features (e.g., leads). The indices developed circumvent the requirement that both the observations and model variables be commensurate (i.e., measured with the same units) by considering the frequencies of the features of interest/importance. We illustrate the metrics by scoring several hypothetical "simulated" discontinuity fields against the lead interpreted from RGPS observations.

Levy, G.↗