Results for “collocation techniques” — Page 12

NASA NTRS2006

Detection of Ice Polar Stratospheric Clouds from Assimilation of Atmospheric Infrared Sounder Data

A novel technique is presented for detection of ice polar stratospheric clouds (PSCs) that form at extremely low temperatures in the lower polar stratosphere during winter. Temperature is a major factor in determining abundance of PSCs, which in turn provide surfaces for heterogeneous chemical reactions leading to ozone loss and radiative cooling. The technique infers the presence of ice PSCs using radiances from the Atmospheric Infrared Sounder (AIRS) in the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. Brightness temperatures are computed from short-term GEOS-5 forecasts for several hundred AIRS channels, using a radiation transfer module. The differences between collocated AIRS observations and these computed values are the observed-minus-forecast (O-F) residuals in the assimilation system. Because the radiation model assumes clear-sky conditions, we hypothesize that these O-F residuals contain quantitative information about PSCs. This is confirmed using sparse data from the Polar Ozone and Aerosol Measurement (POAM) III occultation instrument. The analysis focuses on 0-F residuals for the 6.79pm AIRS moisture channel. At coincident locations, when POAM III detects ice clouds, the AIRS O-F residuals for this channel are lower than -2K. When no ice PSCs are evident in POAM III data, the AIRS 0-F residuals are larger. Given this relationship, the high spatial density of AIRS data is used to construct maps of regions where 0-F residuals are lower than -2K, as a proxy for ice PSCs. The spatial scales and spatio-temporal variations of these PSCs in the Antarctic and Arctic are discussed on the basis of these maps.

Stajner, Ivanka↗

DOE OSTI2024

S-OPT: A Points Selection Algorithm for Hyper-Reduction in Reduced Order Models

While projection-based reduced order models can reduce the dimension of full order solutions, the resulting reduced models may still contain terms that scale with the full order dimension. Hyper-reduction techniques are sampling-based methods that further reduce this computational complexity by approximating such terms with a much smaller dimension. The goal of this work is to introduce the points selection algorithm developed by Shin and Xiu as a hyper-reduction method. The selection algorithm was originally proposed as a stochastic collocation method for uncertainty quantification. Since the algorithm aims at maximizing a quantity $\mathcal{S}$ that measures both the column orthogonality and the determinant, we refer to the algorithm as S-OPT. Numerical examples are provided to demonstrate the performance of S-OPT and to compare its performance with a gappy proper orthogonal decomposition (POD) algorithm. Here, we found that using the S-OPT algorithm is shown to predict the full order solutions with higher accuracy than gappy POD especially when the number of sampling points is small, although we note that S-OPT shows slow asymptotic convergence with respect to the number of samples for some applications, e.g., Lagrangian hydrodynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Dual-frequency (Ka-band and G-band) radar estimates of liquid water content profiles in shallow clouds

The profile of the liquid water content (LWC) in clouds provides fundamental information for understanding the internal structure of clouds, their radiative effects, propensity to precipitate, and degree of entrainment and mixing with the surrounding environment. In principle, differential absorption techniques based on coincident dual-frequency radar reflectivity observations have the potential to provide the LWC profile. Previous differential frequency radar reflectivity (DFR) efforts were challenged by the fact that the measurable differential attenuation for small quantities of LWC is usually comparable to the system measurement error. This typically renders the retrieval impractical, as the uncertainty can become many times greater than the retrieved value itself. Theoretically, this drawback can be mitigated following two interconnected approaches: (1) increasing the frequency separation between the dual-frequency radar system to measure greater differential attenuation and (2) increasing the radar operating frequency to reduce the instrument measurement random error. Our recently developed 239 GHz radar was deployed during the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) along with a variety of collocated remote sensing and in situ instruments. We have combined Ka-band (35 GHz) and G-band (239 GHz) observations to retrieve the LWC from more than 100 vertical profiles of shallow clouds with typical amounts of LWC smaller than 1 g m -3 . We theoretically and experimentally demonstrate that the Ka-band and G-band pair of frequencies offers at least a 65 % relative improvement in the LWC retrieval sensitivity compared to previous works reported in the literature using lower-frequency radars. This new technique provides a missing capability to determine the LWC in the challenging low liquid water path (LWP) range (< 200 g m -2 ) and suggests a way forward to characterize microphysical and dynamical processes more precisely in shallow clouds.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2013

MULTILEVEL ACCELERATION OF STOCHASTIC COLLOCATION METHODS FOR PDE WITH RANDOM INPUT DATA

Stochastic Collocation (SC) methods for stochastic partial differential equa- tions (SPDEs) suffer from the curse of dimensionality, whereby increases in the stochastic dimension cause an explosion of computational effort. To combat these challenges, multilevel approximation methods seek to decrease computational complexity by balancing spatial and stochastic discretization errors. As a form of variance reduction, multilevel techniques have been successfully applied to Monte Carlo (MC) methods, but may be extended to accelerate other methods for SPDEs in which the stochastic and spatial degrees of freedom are de- coupled. This article presents general convergence and computational complexity analysis of a multilevel method for SPDEs, demonstrating its advantages with regard to standard, single level approximation. The numerical results will highlight conditions under which multilevel sparse grid SC is preferable to the more traditional MC and SC approaches.

Webster, Clayton G↗

NASA NTRS1977

Geodetic linear estimation technique and the norm choice problem

In this work the mathematical and probabilistic background of standard linear estimation techniques used in geodesy is clarified, and their interrelationship is revealed with the help of best approximation theory and the normal equations. Emphasis is given to the separation of the deterministic solution to the approximation problem from the probabilistic justification of the metric of the approximation. Least squares prediction has been related to deterministic (exact) collocation, and minimum error bound has been identified as a prediction optimality criterion in the latter. Criteria for the optimal choice of norm in Hilbert space collocation are proposed for gravimetric geodesy problems.

Dermanis, A.↗

NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

NASA NTRS1995

Efficient modeling of interconnects and capacitive discontinuities in high-speed digital circuits

Modeling of interconnects and associated discontinuities with the recent advances high-speed digital circuits has gained a considerable interest over the last decade although the theoretical bases for analyzing these structures were well-established as early as the 1960s. Ongoing research at the present time is focused on devising methods which can be applied to more general geometries than the ones considered in earlier days and, at the same time, improving the computational efficiency and accuracy of these methods. In this thesis, numerically efficient methods to compute the transmission line parameters of a multiconductor system and the equivalent capacitances of various strip discontinuities are presented based on the quasi-static approximation. The presented techniques are applicable to conductors embedded in an arbitrary number of dielectric layers with two possible locations of ground planes at the top and bottom of the dielectric layers. The cross-sections of conductors can be arbitrary as long as they can be described with polygons. An integral equation approach in conjunction with the collocation method is used in the presented methods. A closed-form Green's function is derived based on weighted real images thus avoiding nested infinite summations in the exact Green's function; therefore, this closed-form Green's function is numerically more efficient than the exact Green's function. All elements associated with the moment matrix are computed using the closed-form formulas. Various numerical examples are considered to verify the presented methods, and a comparison of the computed results with other published results showed good agreement.

Oh, K. S.↗

NASA NTRS1976

Higher-order numerical methods derived from three-point polynomial interpolation

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and Hermitian finite-difference discretization. The equations generally apply for both uniform and variable meshes. Hybrid schemes resulting from different polynomial approximations for first and second derivatives lead to the nonuniform mesh extension of the so-called compact or Pade difference techniques. A variety of fourth-order methods are described and this concept is extended to sixth-order. Solutions with these procedures are presented for the similar and non-similar boundary layer equations with and without mass transfer, the Burgers equation, and the incompressible viscous flow in a driven cavity. Finally, the interpolation procedure is used to derive higher-order temporal integration schemes and results are shown for the diffusion equation.

Rubin, S. G.↗

NASA NTRS1991

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

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

Pruett, C. David↗

DOE OSTI2020

A Polarimetric Radar Analysis of Ice Microphysical Processes in Melting Layers of Winter Storms Using S-Band Quasi-Vertical Profiles

Quasi-vertical profiles (QVPs) obtained from a database of U.S. WSR-88D data are used to document polarimetric characteristics of the melting layer (ML) in cold-season storms with high vertical resolution and accuracy. A polarimetric technique to define the top and bottom of the ML is first introduced. Using the QVPs, statistical relationships are developed to gain insight into the evolution of microphysical processes above, within, and below the ML, leading to a statistical polarimetric model of the ML that reveals characteristics that reflectivity data alone are not able to provide, particularly in regions of weak reflectivity factor at horizontal polarization Z H . QVP ML statistics are examined for two regimes in the ML data: Z H ≥ 20 dB Z and Z H < 20 dB Z . Regions of Z H ≥ 20 dB Z indicate locations of MLs collocated with enhanced differential reflectivity Z DR and reduced copolar correlation coefficient ρ hv , while for Z H < 20 dB Z a well-defined ML is difficult to discern using Z H alone. Evidence of large Z DR up to 4 dB, backscatter differential phase δ up to 8°, and low ρ hv down to 0.80 associated with lower Z H (from -10 to 20 dB Z ) in the ML is observed when pristine, nonaggregated ice falls through it. Positive correlation is documented between maximum specific differential phase K DP and maximum Z H in the ML; these are the first QVP observations of K DP in MLs documented at S band. Negative correlation occurs between minimum ρ hv in the ML and ML depth and between minimum ρ hv in the ML and the corresponding enhancement of Z H (Δ Z H = Z H max - Z H rain ).

Meteorology & Atmospheric Sciences↗

DOE OSTI2016

Simplex-stochastic collocation method with improved scalability

The Simplex-Stochastic Collocation (SSC) method is a robust tool used to propagate uncertain input distributions through a computer code. However, it becomes prohibitively expensive for problems with dimensions higher than 5. The main purpose of this paper is to identify bottlenecks, and to improve upon this bad scalability. In order to do so, we propose an alternative interpolation stencil technique based upon the Set-Covering problem, and we integrate the SSC method in the High-Dimensional Model-Reduction framework. In addition, we address the issue of ill-conditioned sample matrices, and we present an analytical map to facilitate uniformly-distributed simplex sampling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

Enhancing Radiative Transfer Models for Optimized Assimilation of Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

NASA NTRSDate not supplied

Optimizing Assimilation of Microwave and Radar Observations in the NWP Models

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

NASA NTRSDate not supplied

Assimilation of Spaceborne Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗

DOE OSTI2014

Parasol: An Architecture for Cross-Cloud Federated Graph Querying

Large scale data fusion of multiple datasets can often provide in- sights that examining datasets individually cannot. However, when these datasets reside in different data centers and cannot be collocated due to technical, administrative, or policy barriers, a unique set of problems arise that hamper querying and data fusion. To ad- dress these problems, a system and architecture named Parasol is presented that enables federated queries over graph databases residing in multiple clouds. Parasol’s design is flexible and requires only minimal assumptions for participant clouds. Query optimization techniques are also described that are compatible with Parasol’s lightweight architecture. Experiments on a prototype implementation of Parasol indicate its suitability for cross-cloud federated graph queries.

Graph database, federated database, data fusion, s↗

NASA NTRS1994

Separation of lithospheric, external, and core components of the south polar geomagnetic field at satellite altitudes

We present a new approach to producing scalar Magsat magnetic anomaly maps based on correlation coefficient filtering and the use of almost all of the available orbits. Our method differs from earlier techniques with respect to the following: (1) Passes are selected based on their variance properties rather than planetary indices such as Kp. (2) The core field model is least squares fit to individual passes and subsequently removed instead of substracting the model directly. This technique replaces band pass filtering or polynomial trend removal methods. (3) Each selected pass is sorted geographically and by local time, placed into one of four different altitude bands, and correlation coefficient filtered with the two adjacent passes. The filtering is the second step toward isolating the static lithospheric signal from the more dynamic external field signals. (4) Least squares collocation is used to grid the correlated passes; subsequently, the dawn and dusk maps are also correlation filtered providing another step toward removal of external fields. (5) The four resultant total field maps are continued to a common altitude and again correlation filtered for the static lithospheric anomalies. (6) The filtered results are then averaged together to provide a new total field map of the lithosphere south of 40 deg S latitude. Our total field map differs from previous efforts over the crustal blocks of West Antarctica. We obtained a positive anomaly over Edward VII Peninsula, extending into the Byrd subglacial basin and obtained a negative anomaly over the Ellsworth Mountains and parts of the Byrd subglacial basin. Also, a positive anomaly extending from the Ross Sea to offshore Wilkes Land is present in our map; however, this feature is absent in other maps. Positive anomalies marking the Weddell Sea in previous efforts are not present in our map. Prominent external field anomalies in the quadrant offshore of Wilkes Land are present in all previous efforts; however, these signatures are reduced in our total field map. This map significantly advances efforts for the identification and interpretation of magnetic anomalies of the south polar lithosphere.

Alsdorf, Douglas E.↗

NASA NTRS2007

An Analysis of Peak Wind Speed Data from Collocated Mechanical and Ultrasonic Anemometers

This study compared peak wind speeds reported by mechanical and ultrasonic anemometers at Cape Canaveral Air Force Station and Kennedy Space Center (CCAFS/KSC) on the east central coast of Florida and Vandenberg Air Force Base (VAFB) on the central coast of California. Launch Weather Officers, forecasters, and Range Safety analysts need to understand the performance of wind sensors at CCAFS/KSC and VAFB for weather warnings, watches, advisories, special ground processing operations, launch pad exposure forecasts, user Launch Commit Criteria (LCC) forecasts and evaluations, and toxic dispersion support. The legacy CCAFS/KSC and VAFB weather tower wind instruments are being changed from propeller-and-vane (CCAFS/KSC) and cup-and-vane (VAFB) sensors to ultrasonic sensors under the Range Standardization and Automation (RSA) program. Mechanical and ultrasonic wind measuring techniques are known to cause differences in the statistics of peak wind speed as shown in previous studies. The 45th Weather Squadron (45 WS) and the 30th Weather Squadron (30 WS) requested the Applied Meteorology Unit (AMU) to compare data between the RSA ultrasonic and legacy mechanical sensors to determine if there are significant differences. Note that the instruments were sited outdoors under naturally varying conditions and that this comparison was not designed to verify either technology. Approximately 3 weeks of mechanical and ultrasonic wind data from each range from May and June 2005 were used in this study. The CCAFS/KSC data spanned the full diurnal cycle, while the VAFB data were confined to 1000-1600 local time. The sample of 1-minute data from numerous levels on five different towers on each range totaled more than 500,000 minutes of data (482,979 minutes of data after quality control). The ten towers were instrumented at several levels, ranging from 12 ft to 492 ft above ground level. The ultrasonic sensors were collocated at the same vertical levels as the mechanical sensors and typically within 15 ft horizontally of each another. Data from a total of 53 RSA ultrasonic sensors, collocated with mechanical sensors were compared. The 1- minute average wind speed/direction and the 1-second peak wind speed/direction were compared.

Short, David A.↗

NASA NTRS1975

High-order numerical solutions using cubic splines

The cubic spline collocation procedure for the numerical solution of partial differential equations was reformulated so that the accuracy of the second-derivative approximation is improved and parallels that previously obtained for lower derivative terms. The final result is a numerical procedure having overall third-order accuracy for a nonuniform mesh and overall fourth-order accuracy for a uniform mesh. Application of the technique was made to the Burger's equation, to the flow around a linear corner, to the potential flow over a circular cylinder, and to boundary layer problems. The results confirmed the higher-order accuracy of the spline method and suggest that accurate solutions for more practical flow problems can be obtained with relatively coarse nonuniform meshes.

Rubin, S. G.↗