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 343 records · Page 19

Variational Algorithms for Nonlinear Smoothing Applications

A variational approach is presented for solving a nonlinear, fixed-interval smoothing problem with application to offline processing of noisy data for trajectory reconstruction and parameter estimation. The nonlinear problem is solved as a sequence of linear two-point boundary value problems. Second-order convergence properties are demonstrated. Algorithms for both continuous and discrete versions of the problem are given, and example solutions are provided.

Bach, Ralph E. Jr.↗

Near-infrared linear array spectrometer for space applications

A 128-element InSb linear detector array read out with a field effect transistor (FET) multiplexing scheme has been developed at the Jet Propulsion Laboratory. State-of-the-art NEP, uniform response, on-detector charge integration, and very low focal plane heat load compared to that of discrete detector arrays allows this type of array to open new possibilities for spaceborne near-infrared imaging spectrometers. Large number of cooled detector elements are feasible with the new technology, making it practical to design sensitive spectrometers even for spinning spacecraft where simultaneous wavelength coverage is necessary to avoid confusion of spatial with spectral variations.

Apt, J.↗

Initialization of mesoscale models - The possible impact of remotely sensed data

Little or no improvement appears to have been achieved in mesoscale numerical prediction. This has been blamed on inherent unpredictability, insufficient spatial resolution, inadequate parameterization of important physical processes, and errors in the numerical discretization of the governing hydrodynamic equations. Attention is presently given to the more fundamental difficulty of forecasting when present conditions are not known with sufficient accuracy, and to the matching of remotely sensed sounding-capability satellite and all-weather Doppler radar data with mesoscale numerical models. This matching is not straightforward; simple and variational forms of four-dimensional assimilation, and Kalman filtering, are suggested approaches but their practical testing remains to be accomplished. The most severe problem faced is that of the temperature and humidity profiles' vertical resolution.

Gal-Chen, T.↗

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗

Large Sulphur Isotope Fractionations in Martian Sediments at Gale Crater

Variability in the sulfur isotopic composition in sediments can reflect atmospheric, geologic and biological processes. Evidence for ancient fluvio-lacustrine environments at Gale crater on Mars and a lack of efficient crustal recycling mechanisms on the planet suggests a surface environment that was once warm enough to allow the presence of liquid water, at least for discrete periods of time, and implies a greenhouse effect that may have been influenced by sulfur-bearing volcanic gases. Here we report in situ analyses of the sulfur isotopic compositions of SO2 volatilized from ten sediment samples acquired by NASA's Curiosity rover along a 13 km traverse of Gale crater. We find large variations in sulfur isotopic composition that exceed those measured for Martian meteorites and show both depletion and enrichment in S-34. Measured values of δS-34 range from -47 +/- 14% to 28 +/- 7%, similar to the range typical of terrestrial environments. Although limited geochronological constraints on the stratigraphy traversed by Curiosity are available, we propose that the observed sulfur isotopic signatures at Gale crater can be explained by equilibrium fractionation between sulfate and sulfide in an impact-driven hydrothermal system and atmospheric processing of sulfur-bearing gases during transient warm periods.

Curiosity↗

Space and Time Variations and Turbopause Dynamical Structure

Time spectra of the critical frequency f sub o E sub s and blanketing frequency f sub b E sub s for the records of a spaced chain of ionosphere vertical sounding stations were studied by means of a maximum entropy method, while coherence spectra were analyzed with the Blackman and Tukey method. It was found out that the frequency parameters spectrum of E sub s within the range of 1-10 cycles/h are discrete with one or two maxima. The first, a low frequency maximum with T=40 min, is of larger amplitude, stable and exists permanently; the second one with T = 10 min is less stable, occurs irregularly and more often within f sub o E sub s. Further, the first maximum corresponds to a higher level of coherence which decreases as the distance to the stations increases. It is concluded that the low-frequency maximum is induced by cellular eddies, usually interpreted as wind shears with horizontal dimensions of not more than 300 km. Irregularity, small amplitude, and low coherence in the range of the second maximum are indicative of the fact that turbulence must be its only source.

Ovezgeldiyev, O. G.↗

The onset of Wolf-Rayet wind outflow and the nature of the hot component in the symbiotic nova PU Vulpecula

We have analyzed temporal variations in the far ultraviolet He II (1640), Si IV (1393, 1402), and C IV (1548, 1550) line profiles in eight high dispersion, International Ultraviolet Explorer Short Wavelength Prime spectra of the symbiotic nova PU Vul by comparatively examining these profiles on a common velocity scale. We see clear evidence of the onset of a Wolf-Rayet-like wind outflow from the bloated, contracting white dwarf hot component with terminal velocity of approximately equals -550 to -600 km/s. We have quantitatively analyzed the complicated He II (1640) emission region for the first time and show that the discrete absorption features seen in the He II region occur at precisely the same velocites in each spectrum, thus demonstrating that the absorbing source is steady and not affected by any orbital motion. We demonstrate that there is an underlying He II wind emission feature whose true shape is hidden by superposed absorption due to the foreground red giant wind flowing in front of the white dwarf and abscuring the white dwarf's wind outflow. We present synthetic spectra of He II emission behind an absorbing slab with u = 20 km/s, T = 5000 K, and column densities in the range N = 1 x 10(exp 22) and 1 x 10(exp 23)/sq cm which explain these absorptions. Our analysis of the Si IV and C IV resonance doublets, in velocity space, reveal temporal variations in the profile between 1987 and 1991 with the emergence of clear P Cygni profiles in Si IV by 1990. A nebular emission feature in C III 1909 also appears in the most recent spectra (e.g., SW42538H) while it was absent or extremely weak in the earliest spectra (e.g., SW36332H), thus strengthening evidence that the nebular emission, as seen in permitted and semiforbidden lines, intensities in step with the onset of the hot, fast, wind outflow. We also report the first detection of narrow interstellar (circumbinary shell?) absorption lines near -1 km/s, most strongly in Al III (1854, 1862) and Si IV (1392, 1402). We have carried out a rough quantitative analysis of the He II wind emission by using the theoretical He II Wolf-Rayet profiles of Hamann & Schmutz (1987). We obtain a lower limit to the He II net emission equivalent width of approximately 1 A, a hot component temperature in the range 25,000 to 35,000 K, a hot component radius in the range 5 solar radius to 30 solar radius, a maximum wind velocity of approximately equals -600 km/s and a rough upper limit to the mass-loss rate of less than 1 x 10(exp -5) solar mass/yr. To our knowledge, this is the first quantitative wind analysis, albeit crude, to be carried out for the hot component of a symbiotic nova or symbiotic variable.

Sion, Edward M.↗

Active flutter control using discrete optimal constrained dynamic compensators

A method for synthesizing digital active flutter suppression controllers using the concept of optimal output feedback is presented. A recently developd convergent algorithm is employed to determine constrained control law parameters that minimize an infinite-time discrete quadratic performance index. Low-order compensator dynamics are included in the control law and the compensator parameters are computed along with the output feedback gain as part of the optimization process. An input noise adjustment procedure is used to improve the stability margins of the digital active flutter controller. Results from investigations into sample rate variation, prefilter pole variation, and effects of varying flight condtions are discussed. The study indicates that a digital control law which accommodates computation delay can stabilize the wing with reasonable rms performance and adequate stability margins.

Broussard, J. R.↗

Signal processing for airborne doppler radar detection of hazardous wind shear as applied to NASA 1991 radar flight experiment data

Radar data collected during the 1991 NASA flight tests have been selectively analyzed to support research directed at developing both improved as well as new algorithms for detecting hazardous low-altitude windshear. Analysis of aircraft attitude data from several flights indicated that platform stability bandwidths were small compared to the data rate bandwidths which should support an assumption that radar returns can be treated as short time stationary. Various approaches at detection of weather returns in the presence of ground clutter are being investigated. Non-coventional clutter rejection through spectrum mode tracking and classification algorithms is a subject of continuing research. Based upon autoregressive modeling of the radar return time sequence, this approach may offer an alternative to overcome errors in conventional pulse-pair estimates. Adaptive filtering is being evaluated as a means of rejecting clutter with emphasis on low signal-to-clutter ratio situations, particularly in the presence of discrete clutter interference. An analysis of out-of-range clutter returns is included to illustrate effects of ground clutter interference due to range aliasing for aircraft on final approach. Data are presented to indicate how aircraft groundspeed might be corrected from the radar data as well as point to an observed problem of groundspeed estimate bias variation with radar antenna scan angle. A description of how recorded clutter return data are mixed with simulated weather returns is included. This enables the researcher to run controlled experiments to test signal processing algorithms. In the summary research efforts involving improved modelling of radar ground clutter returns and a Bayesian approach at hazard factor estimation are mentioned.

Baxa, Ernest G., Jr.↗

Numerical analysis of base flowfield at high altitude for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a pressure based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution and turbulence model are two primary factors that influence the accuracy of the base flowfield prediction.

Wang, Ten-See↗

Numerical study of base pressure characteristic curve for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a three-dimensional pressure-based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution had a strong influence on the accuracy of the base flowfield prediction.

Wang, Ten-See↗

The discrete and localized nature of the variable emission from active regions

Using data from the Extreme Ultraviolet (EUV) Spectroheliometer on Skylab, we study the empirical characteristics of the variable emission in active regions. These simultaneous multi-wavelength observations clearly confirm that active regions consist of a complex of loops at different temperatures. The variable emission from this complex has very well-defined properties that can be quantitatively summarized as follows: (1) It is localized predominantly around the footpoints where it occurs at discrete locations. (2) The strongest variability does not necessarily coincide with the most intense emission. (3) The fraction of the area of the footpoints, (delta n)/N, that exhibits variable emission, varies by +/- 15% as a function of time, at any of the wavelengths measured. It also varies very little from footpoint to footpoint. (4) This fractional variation is temperature dependent with a maximum around 10(exp 5) K. (5) The ratio of the intensity of the variable to the average background emission, (delta I)/(bar-I), also changes with temperature. In addition, we find that these distinctive characteristics persist even when flares occur within the active region.

Arndt, Martina Belz↗

On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach

The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows.

Gastaldi, Fabio↗

Newly found evidence of Sun-climate relationships

Solar radiation cycles drive climatic changes intercyclically. These interdecadal changes were detected as variations in solar total irradiances over the time period of recorded global surface-air-temperature (SAT) and have been restored utilizing Earth Radiation Budget Channel 10C measurements (1978-1990), Greenwich Observatory faculae data (1874-1975), and Taipei Observatory Active Region data (1964-1991). Analysis of the two separate events was carried out by treating each as a discrete time series determined by the length of each solar cycle. The results show that the global SAT responded closely to the input of solar cyclical activities, S, with a quantitative relation of T = 1.62 * S with a correlation coefficient of 0.61. This correlation peaks at 0.71 with a built-in time lag of 32 months in temperature response. Solar forcing in interannual time scale was also detected and the derived relationship of T = 0.17 * S with a correlation coefficient of 0.66 was observed. Our analysis shows derived climate sensitivities approximately fit the theoretical feedback slope, 4T(sup 3).

Kim, Hongsuk H.↗

Experimental Feedback Control of Flow Induced Cavity Tones

Discrete-time, linear quadratic methods were used to design feedback controllers for reducing tones generated by flow over a cavity. The dynamics of a synthetic jet actuator mounted at the leading edge of the cavity as observed by two microphones in the cavity were modeled over a broad frequency range using state space models computed from experimental data. Variations in closed loop performance as a function of model order, control order, control bandwidth, and state estimator design were studied using a cavity in the Probe Calibration Tunnel at NASA Langley. The controller successfully reduced the levels of multiple cavity tones at the tested flow speeds of Mach 0.275, 0.35, and 0.45. In some cases, the closed loop results were limited by excitation of sidebands of the cavity tones, or the creation of new tones at frequencies away from the cavity tones. Nonetheless, the results validate the combination of optimal control and experimentally-generated state space models, and suggest this approach may be useful for other flow control problems. The models were not able to account for non-linear dynamics, such as interactions between tones at different frequencies.

Cabell, Randolph H.↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Aeroheating Predictions for a Hypersonic, Turbulent Near-Wake

The accuracy of heating predictions using various turbulence models is examined for an axisymmetric near-wake at Mach 6. The CFD predictions are compared with experimental data collected under AGARD Working Group 18 on the wake of a 70-degree sphere-cone. The impact of grid resolution and discretization error is estimated, which allows a comparison of stacked-block and conventional structured meshes. The accuracy of steady Reynolds-averaged Navier-Stokes (RANS) models is contrasted with that of a hybrid RANS/Large-Eddy Simulation model. The predictions are made with three different CFD codes (LAURA, FUN3D, and HyperSolve), to demonstrate the code-to-code variation in the results. Steady SST models capture the qualitative nature of the heating in the wake, including the increasing peak heating with increasing Reynolds number. Spalart-Allmaras models, including SA-Catris, under-predicted the peak heating in the wake. Hybrid RANS/LES models improve upon the SA results but have their own modeling difficulties near the shear layer impingement. These results are generally consistent across solvers and grid topologies.

RANS↗

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗