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 73 records · Page 4

Interactions Between Thresholds and Spatial Discretizations of Snow: Insights From Estimates of Wolverine Denning Habitat in the Colorado Rocky Mountains

Thresholds can be used to interpret environmental data in a way that is easily communicated and useful for decision making purposes. However, thresholds are often developed for specific data products and time periods, changing findings when the same threshold is applied to datasets or periods with different characteristics. Here, we test the impact of different spatial discretizations of snow on annual estimates of wolverine denning opportunities in the Colorado Rocky Mountains, defined using a snow water equivalent (SWE) threshold (0.20 m) and threshold date (15 May) from previous habitat assessments. Annual potential wolverine denning area (PWDA) was thresholded from a 36-year (1985 – 2020) snow reanalysis model with three different spatial discretizations: 1) 480 m grid cells (D480), 2) 90 m grid cells (D90), and 3) 480 m grid cells with implicit representations of subgrid snow spatial heterogeneity (S480). Relative to the D480 and S480 discretizations, D90 resolved shallower snow deposits on slopes between 3050 and 3350 m elevation, decreasing PWDA by 10%, on average. In years with warmer and/or drier winters, S480 discretizations with subgrid representations of snow heterogeneity increased PWDA, even within grid cells where mean 15 May SWE was less than the SWE threshold. These simulations increased PWDA by upwards of 30% in low snow years, as compared to the D480 and D90 simulations without subgrid snow heterogeneity. Despite PWDA sensitivity to different snow spatial discretizations, PWDA was controlled more by annual variations in winter precipitation and temperature. However, small changes to the SWE threshold (± 0.07 m) and threshold date (± 2 weeks) also affected PWDA by as much as 82%. Across these threshold ranges, PWDA was approximately 18% more sensitive to the SWE threshold than the threshold date. However, the sensitivity to the threshold date was larger in years with late spring snowfall, when PWDA depended on whether modeled SWE was thresholded before, during, or after spring snow accumulation. Our results demonstrate that snow thresholds are useful but may not always provide a complete picture of the annual variability in snow-adapted wildlife denning opportunities. Studies thresholding spatiotemporal datasets could be improved by including 1) information about the fidelity of thresholds across multiple spatial discretizations, and 2) uncertainties related to ranges of realistic thresholds.

Justin M. Pflug↗

Images in the rocket ultraviolet - Photometry of M101

It is pointed out that ultraviolet imagery of spiral galaxies is a sensitive probe of the large-scale morphology of the star-forming regions and the associated interstellar dust. A rocket-borne Ritchey-Chretien telescope was employed to obtain three ultraviolet images of the Sc 1 galaxy M101 (NGC 5457). The present investigation is concerned with the procedures for reduction and absolute calibration of the obtained images and an analysis in terms of discrete sources. Attention is also given to the radial variation and the azimuthal variation of disk and arms. It is found that the emission from M101 in the UV bandpass is dominated by extreme Population I, in the form of discrete OB/H II complexes.

Hill, J. K.↗

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗

Shock modelling in transonic and supersonic flow

Numerical simulation studies of problems in transonic flow are conducted in connection with the design of aircraft flying near the speed of sound, while recently there has also been a great deal of interest in supersonic and hypersonic flow calculations. Three important models for inviscid compressible gasdynamics are utilized in the study of these flow problems, taking into account the transonic small disturbance equation (TSD), the transonic full potential equation (FP), and the full Euler equations for compressible gasdynamics (EU). The present paper has the objective to provide a description of shock capturing methods recently developed for each model. The design principles behind the algorithms are also examined, giving attention to consistency, conservative form, monotone and sharp discrete shocks, and a bound on the variation, which prevents spurious oscillations.

Osher, S.↗

The influence of sampling rate on the dynamics and control of aeronautical systems

It is shown in the design of aeronautical control systems that sampling rate variations give rise to migrations of the discrete eigenvalues for both the plant and the controllers. In general, these migrations tend to degrade simulation fidelity and the overall performance of Ultrareliable Fault Tolerant Control Systems (UFTCS). These systems consist of a redundant asynchronous network of microprocessors which use sensor information to generate the servo commands that drive actuators, thrust linkages, etc. A change in the sampling period of a single microprocessor will, in addition to causing an eigenvalue shift, alter the magnitude of the servo command issued by the controller. The influence of these changes and the effect of initial conditions and skew are considered.

Howard, J. C.↗

Acceleration techniques for dependability simulation

As computer systems increase in complexity, the need to project system performance from the earliest design and development stages increases. We have to employ simulation for detailed dependability studies of large systems. However, as the complexity of the simulation model increases, the time required to obtain statistically significant results also increases. This paper discusses an approach that is application independent and can be readily applied to any process-based simulation model. Topics include background on classical discrete event simulation and techniques for random variate generation and statistics gathering to support simulation.

Barnette, James David↗

Transition in a Supersonic Boundary-Layer Due to Roughness and Acoustic Disturbances

The transition process induced by the interaction of an isolated roughness with acoustic disturbances in the free stream is numerically investigated for a boundary layer over a flat plate with a blunted leading edge at a free stream Mach number of 3.5. The roughness is assumed to be of Gaussian shape and the acoustic disturbances are introduced as boundary condition at the outer field. The governing equations are solved using the 5'h-~rder accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third- order total-variation-diminishing (TVD) Runge- Kutta scheme for time integration. The steady field induced by the two and three-dimensional roughness is also computed. The flow field induced by two-dimensional roughness exhibits different characteristics depending on the roughness heights. At small roughness heights the flow passes smoothly over the roughness, at moderate heights the flow separates downstream of the roughness and at larger roughness heights the flow separates upstream and downstream of the roughness. Computations also show that disturbances inside the boundary layer is due to the direct interaction of the acoustic waves and isolated roughness plays a minor role in generating instability waves.

Balakumar, P.↗

Effects of Wall Cooling on Hypersonic Boundary Layer Receptivity Over a Cone

Effects of wall cooling on the receptivity process induced by the interaction of slow acoustic disturbances in the free-stream are numerically investigated for a boundary layer flow over a 5-degrees straight cone. The free-stream Mach number is 6.0 and the Reynolds number is 7.8x10(exp 6)/ft. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using 3rd-order total variation diminishing (T VD) Runge-K utta scheme for time integration. Computations are performed for a cone with nose radius of 0.001 inch for adiabatic wall temperature (T(sub aw)), 0.75*T(sub aw), 0.5*T(sub aw), 0.40*T(sub aw), 0.30*T(sub aw), and 0.20*T(sub aw). Once the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations showed that wall cooling has strong stabilization effect on the first mode disturbances as was observed in the experiments. T ransition location moved to upstream when wall cooling was applied It is also found that the boundary layer is much more receptive to fast acoustic wave (by almost a factor of 50). When simulations performed using the same forcing frequency growth of the second mode disturbances are delayed with wall cooling and they attained values two times higher than that of adiabatic case. In 0.20*T(sub aw) case the transition Reynolds number is doubled compared to adiabatic conditions. The receptivity coefficient for adiabatic wall case (804 R) is 1.5225 and for highly cooled cones (241, and 161 R); they are in the order of 10(exp -3).

Kara, K.↗

Development of Logistics for Building Radiation Storm Shelters and Their Operational Evaluation

Over the past three years NASA has been studying the operational effectiveness and astronaut protection efficacy of numerous radiation protection shelters for use in space exploration activities outside of earth's magnetosphere. The work presented was part of NASA's Advanced Exploration Systems (AES) RadWorks Storm Shelter project. This paper is a summary of the concept development activities of this third year. Fabricated items were integrated into mock up deep space habitat vehicle sections for operational evaluations. Two full scale human-in-loop simulations were designed, fabricated, and implemented through an Institutional Review Board approved solicited participant assessment process. Fabricated items are described, along with usage scenarios of two protection approaches. Existing ISS type logistics along with proposed variations of those logistics were used. Preliminary Discrete Event Simulation (DES) work is noted to be useful in quantifying and documenting operational performance measures for the two primary shelter methods, including some characterization of radiation dose accumulation over a mission timeline. The project also performed correlation analyses between effective radiation dose and the Risk of Exposure Induced Death (REID) to show that concept level work may be able to include such a performance metric in early stages of mission scenario habitat design trade space investigation.

Cerro, Jeffrey A.↗

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

The Delta band dissociation of nitric oxide - A potential mechanism for coupling thermospheric variations to the mesosphere and stratosphere

Depletion of solar radiation at discrete wavelengths by nitric oxide has a significant impact on the dissociation rate of this gas in the delta(0, 0) and delta(1, 0) bands. Inclusion of the opacity provided by a typical NO profile reduces the dissociation rate in the upper stratosphere to 50-75 pct of that predicted when the optical depth calculation omits this contribution. A substantial fraction of the NO column abundance as measured near the stratopause resides in the lower thermosphere where correlations of NO with solar and magnetic activity are well documented. Variations in the thermospheric NO abundance therefore modulate the radiation field at the precise wavelengths absorbed by this molecule in the mesosphere and upper stratosphere. Predicted changes in the dissociation rate arising from a varying thermospheric opacity exceed 8 pct throughout the mesosphere and reach 15 pct between 65 km and 95 km.

Frederick, J. E.↗

The response of the dayside aurora to sharp northward and southward transitions of the interplanetary magnetic field and to magnetospheric substorms

Latitudinal shifts of the discrete dayside aurora are examined with respect to variations in the interplanetary magnetic field (IMF) and to magnetospheric substorms. Within 10-15 min after steplike southward (northward) transitions in the IMF the dayside auroral oval moves equatorward (poleward). In at least some cases the auroral shift is similar to an exponential relaxation in latitude from an initial to a final steady state value; for these cases the average exponential time constant is estimated to be 17 min. Substorm features in the dayside aurora include (1) an equatorward shift by 1-3 deg of the equatorward discrete auroral boundary, (2) a brightening of the discrete aurora near substorm onset, (3) the formation of multiple auroral bands, and (4) poleward motion of short-lived individual auroral arcs occurring nearly coincident in time with the equatorward boundary shifts, the average poleward velocity of these arcs being approximately 800 m/s.

Horwitz, J. L.↗

Existence and stability, and discrete BB and rank conditions, for general mixed-hybrid finite elements in elasticity

In this paper, all possible forms of mixed-hybrid finite element methods that are based on multi-field variational principles are examined as to the conditions for existence, stability, and uniqueness of their solutions. The reasons as to why certain 'simplified hybrid-mixed methods' in general, and the so-called 'simplified hybrid-displacement method' in particular (based on the so-called simplified variational principles), become unstable, are discussed. A comprehensive discussion of the 'discrete' BB-conditions, and the rank conditions, of the matrices arising in mixed-hybrid methods, is given. Some recent studies aimed at the assurance of such rank conditions, and the related problem of the avoidance of spurious kinematic modes, are presented.

Xue, W.-M.↗

A comparison of digital flight control design methods

Many variations in design methods for aircraft digital flight control have been proposed in the literature. In general, the methods fall into two categories: those where the design is done in the continuous domain (or s-plane), and those where the design is done in the discrete domain (or z-plane). This paper evaluates several variations of each category and compares them for various flight control modes of the Langley TCV Boeing 737 aircraft. Design method fidelity is evaluated by examining closed loop root movement and the frequency response of the discretely controlled continuous aircraft. It was found that all methods provided acceptable performance for sample rates greater than 10 cps except the 'uncompensated s-plane design' method which was acceptable above 20 cps. A design procedure based on optimal control methods was proposed that provided the best fidelity at very slow sample rates and required no design iterations for changing sample rates.

Powell, J. D.↗

Forced Boundary-Layer Transition on X-43 (Hyper-X) in NASA LaRC 20-Inch Mach 6 Air Tunnel

Aeroheating and boundary layer transition characteristics for the X-43 (Hyper-X) configuration have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel. Global surface heat transfer distributions, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 15.4 million; and inlet cowl door both open and closed. The effects of discrete roughness elements on the forebody boundary layer, which included variations in trip configuration and height, were investigated. This document is intended to serve as a release of preliminary data to the Hyper-X program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

Berry, Scott A.↗

Forced Boundary-Layer Transition on X-43 (Hyper-X) in NASA LaRC 31-Inch Mach 10 Air Tunnel

Aeroheating and boundary layer transition characteristics for the X-43 (Hyper-X) configuration have been experimentally examined in the Langley 31-Inch Mach 10 Air Tunnel. Global surface heat transfer distributions, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, 3-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 5.1 million; and inlet cowl door both open and closed. The effects of discrete roughness elements on the forebody boundary layer, which included variations in trip configuration and height, were investigated. This document is intended to serve as a release of preliminary data to the Hyper-X program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

Berry, Scott A.↗

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗