Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational discretization”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

Discrete ordinate theory of radiative transfer. 2: Scattering from maritime haze

Discrete ordinate theory was used to calculate the reflected and transmitted radiance of photons which have interacted with plane parallel maritime haze layers. The results are presented for three solar zenith angles, three values of the surface albedo, and a range of optical thicknesses from very thin to very thick. The diffuse flux at the lower boundary and the cloud albedo were tabulated. The forward peak and other features in the single scattered phase function caused the radiance in many cases to be very different from that for Rayleigh scattering. The variation of the radiance with both the zenith or nadir angle and the azimuthal angle is more marked, and the relative limb darkening under very thick layers is greater, for haze than for Rayleigh scattering. The downward diffuse flux at the lower boundary for A = O is always greater and the cloud albedo is always less for haze than for Rayleigh layers.

Kattawar, G. W.↗
NASA NTRS1997

X-Ray Flare Characteristics in the B2e Star Lambda Eridani (ROSAT)

We document the results of a simultaneous wavelength monitoring on the B2e star (lambda) Eri. This campaign was carried out from ground stations and with the ROSAT, ASCA, IUE, and Voyager 2 space platforms during a week in February-March 1995; a smaller follow-up was conducted in September 1995. During the first of these intervals (lambda) Eri exhibited extraordinary wind and disk-ejection activity. The ROSAT/HRI X-ray light curves showed no large flares such as the one the ROSAT/PSCA observed in 1991. However, possible low level fluctuations in the February-March ROSAT data occurred at the same time as unusual activity in H(alpha) He I (lambda)6678, He II (lambda)1640, and the C IV doublet. For example, the hydrogen and helium lines exhibited an emission in the blue half of their profiles, probably lasting several hours. The C IV lines showed a strong high-velocity Discrete Absorption Component (DAC) accompanied by unusually strong absorption at lower velocities. The helium line activity suggests that a mass ejection occurred at the base of the wind while the strong C III (Voyager) and C IV (IUE) lines implies that shock interactions occurred in the wind flow. It is not clear that the X-ray elevations are directly related to the strong C IV absorptions because the former changed on a much more rapid timescale than absorptions in the C IV lines. Within hours of the mild X-ray flux variations found by ROSAT on February 28, the Voyager UVS observed a "ringing" that decayed over three 3-hr. cycles. The amplitude of these fluctuations was strong (50%) at (lambda)(lambda)950-1100, decreased rapidly with wavelength, and faded to nondetection longward of (lambda)1300. Various considerations indicate that these continuum variations were not due to an instrumental pathology in the UVS. Rather, they appear to be due to a time-dependent flux deficit in the (lambda)(lambda)950-1250 region. We outline a scenario in which a dense plasma structure over the star's surface is heated and cooled quasi-periodically to produce such flux changes. Observations of new examples of this phenomenon are badly needed. Amateur astronomers can make a significant contribution to its understanding by searching for ringing in light curves of Be stars during their outburst phases. Finally we draw attention to an increase in the emission of the H(alpha) line that occurred at about the time the FUV ringing started. This increased emission hints that approximately 50,000K plasma near the star's surface can influence the circumstellar disc at approximately 12R. by its increased Lyman continuum flux.

Smith, Myron A.↗
NASA NTRS1976

Numerical comparison of discrete Kalman filter algorithms - Orbit determination case study

Numerical characteristics of various Kalman filter algorithms are illustrated with a realistic orbit determination study. The case study of this paper highlights the numerical deficiencies of the conventional and stabilized Kalman algorithms. Computational errors associated with these algorithms are found to be so large as to obscure important mismodeling effects and thus cause misleading estimates of filter accuracy. The positive result of this study is that the U-D covariance factorization algorithm has excellent numerical properties and is computationally efficient, having CPU costs that differ negligibly from the conventional Kalman costs. Accuracies of the U-D filter using single precision arithmetic consistently match the double precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity to variations in the a priori statistics.

Bierman, G. J.↗
NASA NTRS1980

Newton's method: A link between continuous and discrete solutions of nonlinear problems

Newton's method for nonlinear mechanics problems replaces the governing nonlinear equations by an iterative sequence of linear equations. When the linear equations are linear differential equations, the equations are usually solved by numerical methods. The iterative sequence in Newton's method can exhibit poor convergence properties when the nonlinear problem has multiple solutions for a fixed set of parameters, unless the iterative sequences are aimed at solving for each solution separately. The theory of the linear differential operators is often a better guide for solution strategies in applying Newton's method than the theory of linear algebra associated with the numerical analogs of the differential operators. In fact, the theory for the differential operators can suggest the choice of numerical linear operators. In this paper the method of variation of parameters from the theory of linear ordinary differential equations is examined in detail in the context of Newton's method to demonstrate how it might be used as a guide for numerical solutions.

Thurston, G. A.↗
NASA NTRS1988

Radio flares from AE Aquarii - A low-power analog to Cygnus X-37

The magnetic cataclysmic variable AE Aquarii has been observed at 1.5, 4.9, 15, and 22.5 GHz. The source shows temporal variations in the flux density, often with a flarelike morphology. The flux variations are at times extreme, with values ranging from 1-12 mJy at 15 GHz; the degree of variability increases with frequency. The spectrum shows a variety of short-term behavior, at times increasing as nu, at other times decreasing with nu, and at still others flat. The long-term average of the spectrum increases as nu exp 0.3-0.4. No significant degree of circular polarization was observed. The source properties can be explained as a superposition of flare events in which electrons are accelerated to relativistic energies which then emit synchrotron radiation in discrete, expanding sources. The origin of the flare events may be connected to the disruption of the accretion disk some ditance from the white dwarf primary. The radio emission of AE Aqr shows similarities to the low-level flares from Cygnus X-3.

Bastian, T. S.↗
NASA NTRS2010

Statistical and Probabilistic Extensions to Ground Operations' Discrete Event Simulation Modeling

NASA's human exploration initiatives will invest in technologies, public/private partnerships, and infrastructure, paving the way for the expansion of human civilization into the solar system and beyond. As it is has been for the past half century, the Kennedy Space Center will be the embarkation point for humankind's journey into the cosmos. Functioning as a next generation space launch complex, Kennedy's launch pads, integration facilities, processing areas, launch and recovery ranges will bustle with the activities of the world's space transportation providers. In developing this complex, KSC teams work through the potential operational scenarios: conducting trade studies, planning and budgeting for expensive and limited resources, and simulating alternative operational schemes. Numerous tools, among them discrete event simulation (DES), were matured during the Constellation Program to conduct such analyses with the purpose of optimizing the launch complex for maximum efficiency, safety, and flexibility while minimizing life cycle costs. Discrete event simulation is a computer-based modeling technique for complex and dynamic systems where the state of the system changes at discrete points in time and whose inputs may include random variables. DES is used to assess timelines and throughput, and to support operability studies and contingency analyses. It is applicable to any space launch campaign and informs decision-makers of the effects of varying numbers of expensive resources and the impact of off nominal scenarios on measures of performance. In order to develop representative DES models, methods were adopted, exploited, or created to extend traditional uses of DES. The Delphi method was adopted and utilized for task duration estimation. DES software was exploited for probabilistic event variation. A roll-up process was used, which was developed to reuse models and model elements in other less - detailed models. The DES team continues to innovate and expand DES capabilities to address KSC's planning needs.

Trocine, Linda↗
NASA NTRS1975

Long-range energy-state maneuvers for minimum time to specified terminal conditions

The calculus of variations and energy-state approximation are used to determine the optimum maneuvers in three-dimensional minimum-time aircraft flight paths to a specified final point or line. Constraints on thrust, Mach number, angle of attack, dynamic pressure, and load factor are included. It is shown that when the initial range is sufficiently large that the maximum velocity constraint is encountered en route, the calculation of minimum-time maneuvers can be greatly simplified by a separation of arcs into two two-parameter problems. Suboptimal paths along which the bank angle is restricted to three discrete values (negative maximum angle, 0, positive maximum angle) compare favorably with the optimum, continuous bank angle solutions for transonic speeds and below.

Parsons, M. G.↗
NASA NTRS1978

An efficient, direct finite difference method for computing sound propagation in arbitrarily shaped two-dimensional and axisymmetric ducts without flow

An efficient, direct finite difference method is presented for computing sound propagation in non-stepped two-dimensional and axisymmetric ducts of arbitrarily varying cross section without mean flow. The method is not restricted by axial variation of acoustic impedance of the duct wall linings. The non-uniform two-dimensional or axisymmetric duct is conformally mapped numerically into a rectangular or cylindrical computational domain using a new procedure based on a method of fast direct solution of the Cauchy-Riemann equations. The resulting Helmholtz equation in the computational domain is separable. The solution to the governing equation and boundary conditions is expressed as a linear combination of fundamental solutions. The fundamental solutions are computed only once for each duct shape by means of the fast direct cyclic reduction method for the discrete solution of separable elliptic equations. Numerical results for several examples are presented to show the applicability and efficiency of the method.

Chakravarthy, S.↗
NASA NTRS1988

On the convergence of difference approximations to scalar conservation laws

A unified treatment is given for time-explicit, two-level, second-order-resolution (SOR), total-variation-diminishing (TVD) approximations to scalar conservation laws. The schemes are assumed only to have conservation form and incremental form. A modified flux and a viscosity coefficient are introduced to obtain results in terms of the latter. The existence of a cell entropy inequality is discussed, and such an equality for all entropies is shown to imply that the scheme is an E scheme on monotone (actually more general) data, hence at most only first-order accurate in general. Convergence for TVD-SOR schemes approximating convex or concave conservation laws is shown by enforcing a single discrete entropy inequality.

Osher, Stanley↗
NASA NTRS2021

Solving Two-Equation Turbulence Models With a Perspective on Solving Transport Equations

There are three principal objectives of this report. The first objective is to investigate frequently used two-equation models when solving the RANS equations. In this study, we consider the 2006Wilcox and the 2003 Menter Shear Stress Transport (SST) models. Also, a simple change of the Wilcox model is introduced to improve computational predictions for transonic flows. Computations for flows over two different airfoils are examined to compare these turbulence models. Effects on modeling the physics of each flow due to variations in the models, such as using either strainrate or vorticity in the turbulence production term, are considered and discussed. The second objective of this report is to not only explore but also evaluate the performance of the solution algorithm for the mean flow and the transport equations. The RANS and turbulence modeling equations aresolved in a weakly coupled manner with a diagonal implicit Runge-Kutta (DIRK) solution algorithm. Throughout this report, emphasis is given to reducing the residuals to machine zero in all flow calculations, so as to eliminate the error due to numerical integration of the discrete governing equations. The final objective is to provide a perspective on solving transport equations for turbulence modeling. Aspects of solving such stiff systems of equations, as well as various numerical difficulties and possible techniques to overcome them, are considered. Discussion is also provided concerning what is called ’numerical compatibility’, which is an essential requirement when designing a solution algorithm for solving transport equations.

turbulence modeling↗
NASA NTRS2005

Discrete element simulations of gravitational volcanic deformation: Deformation structures and geometries - 1

We have carried out two-dimensional particle dynamics simulations of granular piles subject to frictional Coulomb failure criteria to gain a first-order understanding of different modes of gravitational deformation within volcanoes. Under uniform basal and internal strength conditions, granular piles grow self-similarly, developing distinctive stratigraphies, morphologies, and structures. Piles constructed upon cohesive substrates exhibit particle avalanching, forming outward dipping strata and angle of repose slopes. Systematic decreases in basal strength lead to progressively deeper and steeper internal detachment faults and slip along a basal decollement; landslide forms grade from shallow slumps to deep-seated landslide and, finally, to axial subsidence and outward flank displacements, or volcanic spreading. Surface slopes decrease and develop concave up morphologies with decreasing decollement strength; depositional layers tilt progressively inward. Spatial variations in basal strength cause lateral transitions in pile structure, stratigraphy, and morphology. This approximation of volcanoes as Coulomb granular piles reproduces the richness of deformational structures and surface morphologies in many volcanic settings. The gentle slopes of Hawaiian volcanoes and Olympus Mons on Mars suggest weak basal decollements that enable volcanic spreading. High-angle normal faults, favored above weak decollements, are interpreted in both settings and may explain catastrophic sector collapse in Hawaii and broad aureole deposits surrounding Olympus Mons. In contrast, steeper slopes and shallow detachment faults predominate in the Canary Islands, thought to lack a weak decollement, favoring smaller, more frequent slope failures than predicted for Hawaii. The numerical results provide a useful predictive tool for interpreting dynamic behavior and associated geologic hazards of active volcanoes.

Morgan, Julia K.↗
NASA NTRS1976

A numerical comparison of discrete Kalman filtering algorithms: An orbit determination case study

The numerical stability and accuracy of various Kalman filter algorithms are thoroughly studied. Numerical results and conclusions are based on a realistic planetary approach orbit determination study. The case study results of this report highlight the numerical instability of the conventional and stabilized Kalman algorithms. Numerical errors associated with these algorithms can be so large as to obscure important mismodeling effects and thus give misleading estimates of filter accuracy. The positive result of this study is that the Bierman-Thornton U-D covariance factorization algorithm is computationally efficient, with CPU costs that differ negligibly from the conventional Kalman costs. In addition, accuracy of the U-D filter using single-precision arithmetic consistently matches the double-precision reference results. Numerical stability of the U-D filter is further demonstrated by its insensitivity of variations in the a priori statistics.

Thornton, C. L.↗
NASA NTRS2004

Boundary Layer Control for Hypersonic Airbreathing Vehicles

Active and passive methods for tripping hypersonic boundary layers have been examined in NASA Langley Research Center wind tunnels using a Hyper-X model. This investigation assessed several concepts for forcing transition, including passive discrete roughness elements and active mass addition (or blowing), in the 20-Inch Mach 6 Air and the 31-Inch Mach 10 Air Tunnels. Heat transfer distributions obtained via phosphor thermography, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. The comparisons between the active and passive methods for boundary layer control were conducted at test conditions that nearly match the Hyper-X nominal Mach 7 flight test-point of an angle-of-attack of 2-deg and length Reynolds number of 5.6 million. For passive roughness, the primary parametric variation was a range of trip heights within the calculated boundary layer thickness for several trip concepts. The passive roughness study resulted in a swept ramp configuration, scaled to be roughly 0.6 of the calculated boundary layer thickness, being selected for the Mach 7 flight vehicle. For the active blowing study, the manifold pressure was systematically varied (while monitoring the mass flow) for each configuration to determine the jet penetration height, with schlieren, and transition movement, with the phosphor system, for comparison to the passive results. All the blowing concepts tested, which included various rows of sonic orifices (holes), two- and three-dimensional slots, and random porosity, provided transition onset near the trip location with manifold stagnation pressures on the order of 40 times the model surface static pressure, which is adequate to ensure sonic jets. The present results indicate that the jet penetration height for blowing was roughly half the height required with passive roughness elements for an equivalent amount of transition movement.

Berry, Scott A.↗
NASA NTRS2019

Verification and Validation of a Three-Dimensional Composite Impact Model with Tabulated Input

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in the commercial transient dynamic finite element code LS-DYNA® has been developed. The material model utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure, which allows an arbitrarily shaped failure surface to be defined. A systematic series of validation and verification studies, at a variety of length scales ranging from single element simulations to simulations of a flat panel impact test, have been performed to fully exercise and evaluate the capabilities of the developed model.

Impact↗
NASA NTRS1973

Design, development, and fabrication of a electronic analog microminiaturized electronic analog signal to discrete time interval converter

The microminiaturization of an electronic analog signal to discrete time interval converter is presented. Discrete components and integrated circuits comprising the converter were assembled on a thin-film ceramic substrate containing nichrome resistors with gold interconnections. The finished assembly is enclosed in a flat package measuring 3.30 by 4.57 centimeters. The module can be used whenever conversion of analog to digital signals is required, in particular for the purpose of regulation by means of pulse modulation. In conjunction with a precision voltage reference, the module was applied to control the duty cycle of a switching regulator within a temperature range of -55 C to +125 C, and an input voltage range of 10V to 35V. The output-voltage variation was less than + or - 300 parts per million, i.e., less than + or - 3mV for a 10V output.

Schoenfeld, A. D.↗
NASA NTRS1991

Lake Hoare, Antarctica: sedimentation through a thick perennial ice cover

Lake Hoare in the Dry Valleys of Antarctica is covered with a perennial ice cover more than 3 m thick, yet there is a complex record of sedimentation and of growth of microbial mats on the lake bottom. Rough topography on the ice covering the lake surface traps sand that is transported by the wind. In late summer, vertical conduits form by melting and fracturing, making the ice permeable to both liquid water and gases. Cross-sections of the ice cover show that sand is able to penetrate into and apparently through it by descending through these conduits. This is the primary sedimentation mechanism in the lake. Sediment traps retrieved from the lake bottom indicate that rates of deposition can vary by large amounts over lateral scales as small as 1 m. This conclusion is supported by cores taken in a 3 x 3 grid with a spacing of 1.5 m. Despite the close spacing of the cores, the poor stratigraphic correlation that is observed indicates substantial lateral variability in sedimentation rate. Apparently, sand descends into the lake from discrete, highly localized sources in the ice that may in some cases deposit a large amount of sand into the lake in a very short time. In some locations on the lake bottom, distinctive sand mounds have been formed by this process. They are primary sedimentary structures and appear unique to the perennially ice-covered lacustrine environment. In some locations they are tens of centimetres high and gently rounded with stable slopes; in others they reach approximately 1 m in height and have a conical shape with slopes at angle of repose. A simple formation model suggests that these differences can be explained by local variations in water depth and sedimentation rate. Rapid colonization of fresh sand surfaces by microbial mats composed of cyanobacteria, eukaryotic algae, and heterotrophic bacteria produces a complex intercalation of organic and sandy layers that are a distinctive form of modern stromatolites.

NASA Discipline Exobiology↗
NASA NTRS1972

Adaptive stochastic control for a class of linear systems.

The problem considered in this paper deals with the control of linear discrete-time stochastic systems with unknown (possibly time-varying and random) gain parameters. The philosophy of control is based on the use of an open-loop feedback optimal (OLFO) control using a quadratic index of performance. It is shown that the OLFO system consists of (1) an identifier that estimates the system state variables and gain parameters and (2) a controller described by an 'adaptive' gain and correction term. Several qualitative properties and asymptotic properties of the OLFO adaptive system are discussed. Simulation results dealing with the control of stable and unstable third-order plants are presented. The key quantitative result is the precise variation of the control system adaptive gains as a function of the future expected uncertainty of the parameters; thus, in this problem the ordinary 'separation theorem' does not hold.

Tse, E.↗
NASA NTRS1994

Auroral activity associated with Kelvin-Helmholtz instability at the inner edge of the low-latitude boundary layer

Auroral activity occurred in the late afternoon sector (approx. 16 MLT) in the northern hemisphere during the passage at Earth of an interplanetary magnetic cloud on January 14, 1988. The auroral activity consisted of a very dynamic display which was preceded and followed by quiet auroral displays. During the quiet displays, discrete rayed arcs aligned along the geomagnetic L shells were observed. In the active stage, rapidly evolving spiral forms centered on magnetic zenith were evident. The activity persisted for many minutes and was characterized by the absence of directed motion. They were strongly suggestive of intense filaments of upward field-aligned currents embedded in the large-scale region 1 current system. Distortions of the flux ropes as they connect from the equatorial magnetosphere to the ionosphere were witnessed. We assess as possible generating mechanisms three nonlocal sources known to be associated with field-aligned currents. Of these, partial compressions of the magnetosphere due to variations of solar wind dynamic pressure seem an unlikely source. The possibility that the auroral forms are due to reconnection is investigated but is excluded because the active aurora were observed on the closed field line region just equatorward of the convection reversal boundary. To support this conclusion further, we apply recent results on the mapping of ionospheric regions to the equatorial plane based on the Tsyganenko 1989 model (Kaufmann et al., 1993). We find that for comparable magnetic activity the aurora map to the equatorial plane at X(sub GSM) = approx. 3 R(sub E) and approx. 2 R(sub E) inward of the magnetopause, that is, the inner edge of the boundary layer close to dusk. Since the auroral forms are manifestly associated with magnetic field shear, a vortical motion at the equatorial end of the flux rope is indicated, making the Kelvin-Helmholtz instability acting at the inner edge of the low-latitude boundary layer the most probable generating source.

Farrugia, C. J.↗