Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameterization method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗

Probability of satellite collision

A method is presented for computing the probability of a collision between a particular artificial earth satellite and any one of the total population of earth satellites. The collision hazard incurred by the proposed modular Space Station is assessed using the technique presented. The results of a parametric study to determine what type of satellite orbits produce the greatest contribution to the total collision probability are presented. Collision probability for the Space Station is given as a function of Space Station altitude and inclination. Collision probability was also parameterized over miss distance and mission duration.

Mccarter, J. W.↗

Charge-transfer energy in closed-shell ion-atom interactions

The importance of charge-transfer energy in the interactions between closed-shell ions and atoms is investigated. Ab initio calculations on H(plus)-He and Li(plus)-He are used as a guide for the construction of approximate methods for the estimation of the charge-transfer energy for more complicated systems. For many alkali ion-rate gas systems the charge-transfer energy is comparable to the induction energy in the region of the potential minimum, although for doubly charged alkaline-earth ions in rare gases the induction energy always dominates. Surprisingly, an empirical combination of repulsion energy plus asymptotic induction energy plus asymptotic dispersion energy seems to give a fair representation of the total interaction, especially if the repulsion energy is parameterized, despite the omission of any explicit charge-transfer contribution. More refined interaction models should consider the charge-transfer energy contribution.

Alvarez-Rizzatti, M.↗

Viscous stability theory for thermally stratified flow - Discontinuous jets and shear layers.

The method used by Drazin to investigate the stability of unbounded, viscous, homogeneous, parallel shear flow to small wavenumber disturbances is extended to study the effect of thermal stratification on the stability of unbounded jets and shear layers. By this method the stability characteristics of continuous profiles are inferred from the stability characteristics of discontinuous profiles. The characteristic value problem for discontinuous jet and shear layers is posed by the requirement that the solutions of the governing differential equation satisfy the matching conditions and boundedness conditions for layers that extend to infinity. The analysis leads to a characteristic determinant which is required to vanish for the characteristic values of the parameters: the Reynolds number, the wavenumber, and the wave speed. The stabilizing effect of the thermal stratification as parameterized by the Richardson number was found to be most stabilizing for small wavenumber (large-scale) disturbances.

Gage, K. S.↗

Scattering by non-spherical particles of size comparable to a wavelength - A new semi-empirical theory

An approximate method is proposed for evaluating the interaction of randomly oriented, nonspherical particles with the total intensity component of electromagnetic radiation. When the particle size parameter, x, the ratio of particle circumference to wavelength, is less than some upper bound x(o) (about 5), Mie theory is used. For x greater than x(o), the interaction is divided into three components: diffraction, external reflection, and transmission. Physical optics theory is used to obtain the first of these components; geometrical optics theory is applied to the second; and a simple parameterization is employed for the third. The predictions of this theory are found to be in very good agreement with laboratory measurements for a wide variety of particle shapes, sizes, and refractive indexes. Limitations of the theory are also noted.

Pollack, J. B.↗

Comparison of tropical pacific temperature and current simulations with two vertical mixing schemes embedded in an ocean general circulation model and reference to observations

The Pacanowski-Philander (PP) and Mellor-Yamada (MY) parameterization models of vertical mixing by turbulent processes were embedded in the Geophysical Fluid Dynamics Laboratory high-resolution ocean general circulation model of the tropical Pacific Ocean. All other facets of the numerical simulations were the same. Simulations were made for the 1987-1988 period. At the equator the MY simulation produced near-surface temperatures more uniform with depth, a deeper thermocline, a deeper core speed of the Equatorial Undercurrent, and a South Equatorial Current with greater vertical thickness compared with that computed with the PP method. Along 140 deg W, between 5 deg N and 10 deg N, both simulations were the same. Moored buoy current and temperature observations had been recorded by the Pacific Marine Environmental Laboratory at three sites (165 deg E, 140 deg W, 110 deg W) along the equator and at three sites (5 deg N, 7 deg N, 9 deg N) along 140 deg W. Simulated temperatures were lower than those observed in the near-surface layer and higher than those observed in the thermocline. Temperature simulations were in better agreement with observations compared to current simulations. At the equator, PP current and temperature simulations were more representative of the observations than MY simulations.

Halpern, David↗

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗

Comparison of GOES-Retrieved and In Situ Measurements of Deep Convective Anvil Cloud Microphysical Properties During the Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4))

One of the main goals of the Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) during July and August 2007 was to gain a better understanding of the formation and life cycle of cirrus clouds in the upper troposphere and lower stratosphere and how their presence affects the exchange of water vapor between these layers. Additionally, it is important to compare in situ measurements taken by aircraft instruments with products derived from satellite observations and find a meaningful way to interpret the results. In this study, cloud properties derived using radiance measurements from the Geostationary Operational Environmental Satellite (GOES) imagers are compared to similar quantities from aircraft in situ observations and are examined for meaningful relationships. A new method using dual \angle satellite measurements is used to derive the ice water content (IWC) for the top portion of deep convective clouds and anvils. The results show the in situ and remotely sensed mean microphysical properties agree to within approx.10 microns in the top few kilometers of thick anvils despite the vastly different temporal and spatial resolutions of the aircraft and satellite instruments. Mean particle size and IWC are shown to increase with decreasing altitude in the top few kilometers of the cloud. Given these relationships, it may be possible to derive parameterizations for effective particle size and IWC as a function of altitude from satellite observations

Yost, Christoher R.↗

Curved Finite Elements and Curve Approximation

The approximation of parameterized curves by segments of parabolas that pass through the endpoints of each curve segment arises naturally in all quadratic isoparametric transformations. While not as popular as cubics in curve design problems, the use of parabolas allows the introduction of a geometric measure of the discrepancy between given and approximating curves. The free parameters of the parabola may be used to optimize the fit, and constraints that prevent overspill and curve degeneracy are introduced. This leads to a constrained optimization problem in two varibles that can be solved quickly and reliably by a simple method that takes advantage of the special structure of the problem. For applications in the field of computer-aided design, the given curves are often cubic polynomials, and the coefficient may be calculated in closed form in terms of polynomial coefficients by using a symbolic machine language so that families of curves can be approximated with no further integration. For general curves, numerical quadrature may be used, as in the implementation where the Romberg quadrature is applied. The coefficient functions C sub 1 (gamma) and C sub 2 (gamma) are expanded as polynomials in gamma, so that for given A(s) and B(s) the integrations need only be done once. The method was used to find optimal constrained parabolic approximation to a wide variety of given curves.

Baart, M. L.↗

A general rough-surface inversion algorithm: Theory and application to SAR data

Rough-surface inversion has significant applications in interpretation of SAR data obtained over bare soil surfaces and agricultural lands. Due to the sparsity of data and the large pixel size in SAR applications, it is not feasible to carry out inversions based on numerical scattering models. The alternative is to use parameter estimation techniques based on approximate analytical or empirical models. Hence, there are two issues to be addressed, namely, what model to choose and what estimation algorithm to apply. Here, a small perturbation model (SPM) is used to express the backscattering coefficients of the rough surface in terms of three surface parameters. The algorithm used to estimate these parameters is based on a nonlinear least-squares criterion. The least-squares optimization methods are widely used in estimation theory, but the distinguishing factor for SAR applications is incorporating the stochastic nature of both the unknown parameters and the data into formulation, which will be discussed in detail. The algorithm is tested with synthetic data, and several Newton-type least-squares minimization methods are discussed to compare their convergence characteristics. Finally, the algorithm is applied to multifrequency polarimetric SAR data obtained over some bare soil and agricultural fields. Results will be shown and compared to ground-truth measurements obtained from these areas. The strength of this general approach to inversion of SAR data is that it can be easily modified for use with any scattering model without changing any of the inversion steps. Note also that, for the same reason it is not limited to inversion of rough surfaces, and can be applied to any parameterized scattering process.

Moghaddam, M.↗

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗

A Model-Based, Bayesian Solution for Characterization of Complex Damage Scenarios in Aerospace Composite Structures

Ultrasonic damage detection and characterization is commonly used in nondestructive evaluation (NDE) of aerospace composite components. In recent years there has been an increased development of guided wave based methods. In real materials and structures, these dispersive waves result in complicated behavior in the presence of complex damage scenarios. Model-based characterization methods utilize accurate three dimensional finite element models (FEMs) of guided wave interaction with realistic damage scenarios to aid in defect identification and classification. This work describes an inverse solution for realistic composite damage characterization by comparing the wavenumber-frequency spectra of experimental and simulated ultrasonic inspections. The composite laminate material properties are first verified through a Bayesian solution (Markov chain Monte Carlo), enabling uncertainty quantification surrounding the characterization. A study is undertaken to assess the efficacy of the proposed damage model and comparative metrics between the experimental and simulated output. The FEM is then parameterized with a damage model capable of describing the typical complex damage created by impact events in composites. The damage is characterized through a transdimensional Markov chain Monte Carlo solution, enabling a flexible damage model capable of adapting to the complex damage geometry investigated here. The posterior probability distributions of the individual delamination petals as well as the overall envelope of the damage site are determined.

H. Reed↗

An Efficient GPU-Accelerated Multi-Source Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog ∼12000 Galactic Binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA global fit↗

Krylov methods preconditioned with incompletely factored matrices on the CM-2

The performance is measured of the components of the key interative kernel of a preconditioned Krylov space interative linear system solver. In some sense, these numbers can be regarded as best case timings for these kernels. Sweeps were timed over meshes, sparse triangular solves, and inner products on a large 3-D model problem over a cube shaped domain discretized with a seven point template. The performance of the CM-2 is highly dependent on the use of very specialized programs. These programs mapped a regular problem domain onto the processor topology in a careful manner and used the optimized local NEWS communications network. The rather dramatic deterioration in performance was documented when these ideal conditions no longer apply. A synthetic workload generator was developed to produce and solve a parameterized family of increasingly irregular problems.

Berryman, Harry↗

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign and the 1998 First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). We investigate two related approaches to prediction of GS cloud properties from SGS variability: the critical saturation ratio method, and the critical relative humidity method. We find significant correlation between GS cloud water and SGS supersaturation when the critical saturation ratio is set at 100%, as well as a notable relationship between GS cloud water and the width of the SGS total water distribution. Critical relative humidity also compares well with GS cloud water. However, we also find that the assumptions of a static critical saturation ratio of 100% to be unrealistic, as well as a fixed SGS distribution width. Empirical calculations from the ARISE data show a large sensitivity of these SGS variables to GS relative humidity, and so a SGS parameterization allowing them to vary according to GS thermodynamic properties may result in more realistic GS cloud water values.

J. Brant Dodson↗

Shortcomings with Tree-Structured Edge Encodings for Neural Networks

In evolutionary algorithms a common method for encoding neural networks is to use a tree structured assembly procedure for constructing them. Since node operators have difficulties in specifying edge weights and these operators are execution-order dependent, an alternative is to use edge operators. Here we identify three problems with edge operators: in the initialization phase most randomly created genotypes produce an incorrect number of inputs and outputs; variation operators can easily change the number of input/output (I/O) units; and units have a connectivity bias based on their order of creation. Instead of creating I/O nodes as part of the construction process we propose using parameterized operators to connect to preexisting I/O units. Results from experiments show that these parameterized operators greatly improve the probability of creating and maintaining networks with the correct number of I/O units, remove the connectivity bias with I/O units and produce better controllers for a goal-scoring task.

Hornby, Gregory S.↗

On the Ice Nucleation Spectrum

This work presents a novel formulation of the ice nucleation spectrum, i.e. the function relating the ice crystal concentration to cloud formation conditions and aerosol properties. The new formulation is physically-based and explicitly accounts for the dependency of the ice crystal concentration on temperature, supersaturation, cooling rate, and particle size, surface area and composition. This is achieved by introducing the concepts of ice nucleation coefficient (the number of ice germs present in a particle) and nucleation probability dispersion function (the distribution of ice nucleation coefficients within the aerosol population). The new formulation is used to generate ice nucleation parameterizations for the homogeneous freezing of cloud droplets and the heterogeneous deposition ice nucleation on dust and soot ice nuclei. For homogeneous freezing, it was found that by increasing the dispersion in the droplet volume distribution the fraction of supercooled droplets in the population increases. For heterogeneous ice nucleation the new formulation consistently describes singular and stochastic behavior within a single framework. Using a fundamentally stochastic approach, both cooling rate independence and constancy of the ice nucleation fraction over time, features typically associated with singular behavior, were reproduced. Analysis of the temporal dependency of the ice nucleation spectrum suggested that experimental methods that measure the ice nucleation fraction over few seconds would tend to underestimate the ice nuclei concentration. It is shown that inferring the aerosol heterogeneous ice nucleation properties from measurements of the onset supersaturation and temperature may carry significant error as the variability in ice nucleation properties within the aerosol population is not accounted for. This work provides a simple and rigorous ice nucleation framework where theoretical predictions, laboratory measurements and field campaign data can be reconciled, and that is suitable for application in atmospheric modeling studies.

Barahona, D.↗