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 271 records · Page 15

Stellar convection 3: Convection at large Rayleigh numbers

A three dimensional study of convection in a self gravitating sphere of Boussinesq fluid with a Rayleigh number of 10 to the 10th power and a Prandtl of 1 is presented. The velocity and temperature of the fluid are computed at the largest wavelengths using spectral methods. A confirmation that the fluid is anisotropic and that the energy spectra are not smooth functions of wavelength but have a large amount of fine structure is discussed. The parameterization of the transport properties of the unresolvable inertial subrange with eddy viscosities and diffusivities is described. The time dependent fluctuations in the energy spectra and how they cascade from large to small wavelengths is examined.

Marcus, P. S.↗

The correction of infrasound signals for upper atmospheric winds

Infrasound waves propagate in the atmosphere by a well known mechanism produced by refraction of the waves, return to earth, and reflection at the surface into the atmosphere for subsequent bounces. A figure illustrates this phenomenon with results from a ray trace model. In this instance three rays are returned to earth from a region centered at about 50 kilometers in altitude and two from a region near 110 kilometers in altitude. The control of the wave refraction is largely dominated by the temperature-height profile and inversions; however, a major influence is also produced by the atmospheric wind profile. Another figure illustrates the considerable ray differences for rays moving in the wind direction (to the right) and in the counter direction (to the left). It obviously can be expected that infrasonic signal amplitudes will be greatly influenced by the winds in the atmosphere. The seasonal variation of the high altitude atmospheric winds is well documented. A third figure illustrates this with average statistics on the observed zonal wind in the region of 50 plus or minus 5 kilometers in altitude. The results are based upon a survey by Webb; Webb terms this parameterization the Stratospheric Circulation Index (SCI). The very strong seasonal variation has the ability to exert a major seasonal influence on infrasonic signals. The purpose here is to obtain a method for the correction of this effect.

Mutschlecner, J. Paul↗

Parameter dimension of turbulence-induced phase errors and its effects on estimation in phase diversity

The method of phase diversity has been used in the context of incoherent imaging to estimate jointly an object that is being imaged and phase aberrations induced by atmospheric turbulence. The method requires a parametric model for the phase-aberration function. Typically, the parameters are coefficients to a finite set of basis functions. Care must be taken in selecting a parameterization that properly balances accuracy in the representation of the phase-aberration function with stability in the estimates. It is well known that over parameterization can result in unstable estimates. Thus a certain amount of model mismatch is often desirable. We derive expressions that quantify the bias and variance in object and aberration estimates as a function of parameter dimension.

Thelen, Brian J.↗

A Unified Directional Spectrum for Long and Short Wind-Driven Waves

Review of several recent ocean surface wave models finds that while comprehensive in many regards, these spectral models do not satisfy certain additional, but fundamental, criteria. We propose that these criteria include the ability to properly describe diverse fetch conditions and to provide agreement with in situ observations of Cox and Munk [1954] and Jiihne and Riemer [1990] and Hara et al. [1994] data in the high-wavenumber regime. Moreover, we find numerous analytically undesirable aspects such as discontinuities across wavenumber limits, nonphysical tuning or adjustment parameters, and noncentrosymmetric directional spreading functions. This paper describes a two-dimensional wavenumber spectrum valid over all wavenumbers and analytically amenable to usage in electromagnetic models. The two regime model is formulated based on the Joint North Sea Wave Project (JONSWAP) in the long-wave regime and on the work of Phillips [1985] and Kitaigorodskii [1973] at the high wavenumbers. The omnidirectional and wind-dependent spectrum is constructed to agree with past and recent observations including the criteria mentioned above. The key feature of this model is the similarity of description for the high- and low-wavenumber regimes; both forms are posed to stress that the air-sea interaction process of friction between wind and waves (i.e., generalized wave age, u/c) is occurring at all wavelengths simultaneously. This wave age parameterization is the unifying feature of the spectrum. The spectrum's directional spreading function is symmetric about the wind direction and has both wavenumber and wind speed dependence. A ratio method is described that enables comparison of this spreading function with previous noncentrosymmetric forms. Radar data are purposefully excluded from this spectral development. Finally, a test of the spectrum is made by deriving roughness length using the boundary layer model of Kitaigorodskii. Our inference of drag coefficient versus wind speed and wave age shows encouraging agreement with Humidity Exchange Over the Sea (HEXOS) campaign results.

Elfouhaily, T.↗

Robust Conceptual Design of Transonic Airfoils

This paper describes an integrated, multi-fidelity analysis and heuristic design approach that can be used to derive initial airfoil designs for transonic flight. If successful, the final result is a geometry that can be expected to produce reasonable aerodynamic performance when used with higher order analysis methods. A key aspect of the methodology is the use of a “sonic-plateau” pressure distribution as the target distribution for inverse design. The sonic- plateau distribution is easily parameterized and has the advantage of automatically resulting in a smooth airfoil shape without any discontinuities built into the surface due to the presence of a shock in the target pressure distribution. Inverse design is performed on each airfoil using a parametrically defined pressure distribution at a reduced lift coefficient and Mach number from the design condition. The methodology is demonstrated by designing an airfoil at 38% of the wing semispan for a 737-200-like aircraft. The demonstration problem shows that the methodology is able to achieve rapid and robust convergence to the solution. The calculated designed airfoil was found to be sufficiently higher than Mach number, and the maximum thickness was close to the targeted value.

Olson, Erik D.↗

Rapid and Accurate Computation of Invariant Tori, Manifolds, and Connections Near Mean Motion Resonances in the Periodically Perturbed Planar Circular Restricted 3-Body Problem Models

When the planar circular restricted 3-body problem (RTBP) is periodically perturbed, most unstable resonant periodic orbits become invariant tori. In this study, we 1) develop a quasi-Newton method which simultaneously solves for the tori and their center, stable, and unstable directions; 2) implement continuation by both perturbation as well as rotation numbers; 3) compute Fourier-Taylor parameterizations of the stable and unstable manifolds; 4) globalize these manifolds; 5) compute homoclinic and heteroclinic connections. Our methodology improves on efficiency and accuracy compared to prior studies, and applies to a variety of periodic perturbations. We demonstrate the tools on the planar elliptic RTBP.

de la Llave, Rafael↗

A Comparison of Cumulus Parameterizations in Idealized Sea-Breeze Simulations

Four cumulus parameterizations in the Penn State-NCAR model MM5 are compared in idealized sea-breeze simulations, with the aim of discovering why they work as they do. The most realistic results appear to be those using the Kain-Fritsch scheme. Rainfall is significantly delayed with the Betts-Miller-Janjic scheme, due to the method of computing the reference sounding. This method can be corrected, but downdrafts should be added in a physically realistic manner. Even without downdrafts, a corrected version of the BMJ scheme produces nearly the same timing and location of deep convection as the KF scheme, despite the very different physics. In order to simulate the correct timing of the rainfall, a minimum amount of mass is required in the layer that is the source of a parameterized updraft. The Grell parameterization, in the present simulation, always derives the updraft from the top of the mixed layer, where vertical advection predominates over horizontal advection in increasing the moist static energy. This makes the application of the quasi-equilibrium closure more correct than it would be if the updrafts were always derived from the most unstable layer, but it evades the question of whether or not horizontal advection generates instability. Using different physics, the parameterizations produce significantly different cloud-top heights.

Cohen, Charles↗

Principal normal curvature of surfaces

Certain principal normal curvatures of differential geometry were developed for use in curvature matrices associated with the asymptotic solution of electromagnetic diffraction problems. The effort is directed toward microwave antenna simulations and high speed digital computer analysis of radiometric instruments used to obtain soil moisture, sea state, salinity and temperature data. It is shown that the methods used to develop the principal normal curvatures for paraboloid, hyperboloid, ellipsoid, sphere, and cone can be applied to other radiometer geometries such as the parabolic torus, even though the surface parameterizations are different. It is concluded that deployable offset geometries, distorted by rotational forces and solar loads may be analyzed by similar means given a suitable surface description.

Schmidt, R. F.↗

Optimization of linear, controlled structures

An approach to the integrated design of linear controlled structures, which is being investigated as part of NASA's controlled structures interaction (CSI) methodology development program, is presented. The approach is integrated in the sense that the structure and its controller are simultaneously designed. The design methodology uses constrained nonlinear optimization procedures based on analytically obtained gradients of the structural responses. Design of the controller is based on the so-called Q-parameterization theory, which parameterizes all closed-loop input/output maps achievable with stabiizing linear controllers. Very general objective and constraint functions are possible, and structural shape can be included in the design variables. This method has been partially implemented and demonstrated; early findings are reported.

Hunziker, K. Scott↗

The NASA/GISS Mars general circulation model: Preliminary experiments

The NASA/GISS Mars General Circulation Model (GCM) is an adapted version of the GISS Global Climate/Middle Atmosphere Model, specifically developed for the diagnostic validation and objective analysis of measured atmospheric temperatures from the Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) experiment. The GISS Mars GCM has 23 vertical layers extending from the surface to approximately 80 km altitude, representing a vertical resolution of about 0.3 scale heights. The primitive (vertically hydrostatic) equations are solved in finite difference form on the Krakawa B grid, with a horizontal resolution of 8 deg x 10 deg (latitude-longitude). The model includes a diurnal solar cycle, heat transport within a two-layer ground, and a high-order 'slopes-scheme' for the advection of heat in the upper atmosphere. The radiative transfer scheme is based on the correlated k distribution method for the treatment of nongray gaseous absorption thermal emission, and multiple scattering, including options for suspended dust. A special feature of the model of particular importance for Mars is a parameterization of gravity-wave-induced drag incorporating orographic forcing, wind shear, convection, and radiative damping. The implementation of the GISS Mars model includes global maps of topography, roughness, and albedo.

Allison, Michael↗

Metafitting: Weight optimization for least-squares fitting of PTTI data

For precise time intercomparisons between a master frequency standard and a slave time scale, we have found it useful to quantitatively compare different fitting strategies by examining the standard uncertainty in time or average frequency. It is particularly useful when designing procedures which use intermittent intercomparisons, with some parameterized fit used to interpolate or extrapolate from the calibrating intercomparisons. We use the term 'metafitting' for the choices that are made before a fitting procedure is operationally adopted. We present methods for calculating the standard uncertainty for general, weighted least-squares fits and a method for optimizing these weights for a general noise model suitable for many PTTI applications. We present the results of the metafitting of procedures for the use of a regular schedule of (hypothetical) high-accuracy frequency calibration of a maser time scale. We have identified a cumulative series of improvements that give a significant reduction of the expected standard uncertainty, compared to the simplest procedure of resetting the maser synthesizer after each calibration. The metafitting improvements presented include the optimum choice of weights for the calibration runs, optimized over a period of a week or 10 days.

Douglas, Rob J.↗

Mass Conservation in a Chemical Transport Model and its Effect on CO2 and SF6 Simulations

Chemical transport models (CTMs) must conserve mass to be useful for applications involving assessment of the effect of various pollutants on the troposphere and stratosphere. Furthermore, calculations of the evolution of constituents such as SF6 are used to evaluate overall model transport, and interpretation of such simulations is clouded if mass conservation is not assured. For realistic simulations or predictions, it is crucial that constituents are not produced or lost by transport or other processes in the CTMs. Analysis of CO2 and SF6 experiments using a CTM shows that problems with mass conservation can seriously degrade the simulations. Failure to conserve mass results from inconsistency of the surface pressure tendency and the divergence of horizontal mass flux when the model is forced by assimilated meteorological data. We have developed an effective method to eliminate the inconsistency by modifying the divergent part of the wind field. The changes in the wind fields are quite small but the impact on mass conservation is large. Parameterizations of physical processes such as convection or turbulent transport can also affect mass conservation. The lack of conservation is small but accumulates when integrations are lengthy such as required for SF6. This lack of conservation is found using winds from either a GCM or from an assimilation system. A simple adjustment removes much of the inaccuracy in the convective parameterization. A CO2 simulation using assimilated winds from the most recent version of the Goddard Earth Observing System Data Assimilation System will be used to illustrate the impact of these transport improvements.

Zhu, Z.↗

Synergy of Satellite Radiation, Precipitation, and Other Meteorological Variable Observations for Global Mean Sea Surface Turbulent Heat Flux Estimation

Sea surface turbulent heat flux is one of the key components in global freshwater and energy balances and plays an important role in atmospheric dynamics, thermodynamics and general circulation. This turbulent heat flux is dominantly decided by sea surface latent heat release with some contribution from sensible heat exchange. The flux and its anomaly could significantly affect ocean heat storage and ocean circulation. They are part of and have great impacts on climate variability. Though the turbulent heat flux is extremely important for the climate and weather systems, only very limited ship and buoy observations of the flux are available over open oceans. There are no global, operational, direct measurements of this crucial meteorological variable. Global estimates are, basically, indirectly calculated from bulk turbulent flux parameterization with a combination of satellite column water vapor, sea surface water temperature, and wind speed observations. Critical parameters such as sea surface air temperature and humidity are estimated from empirical relations of theses variables with the water vapor and sea surface water temperature, respectively. Lacking accurate knowledge on surface air temperature and humidity, uncertainties in the parameterization and potential changes in the non-linear turbulent processes with long-term climate variations could cause large errors in estimated long-term turbulent fluxes from the indirect method as shown in current global sea surface turbulent flux datasets. This study uses synergized data of satellite global radiation, precipitation, and other meteorological variable observations to estimate sea surface turbulent heat fluxes. The radiation observations are made by the satellite Clouds and the Earth’s Radiant Energy System (CERES) sensors, while the global precipitation data is from NASA’s satellite Global Precipitation Climatology Project (GPCP). Other data includes satellite sea surface water temperature and wind speed observations. These datasets are obtained from a wide range of space sensors from passive to active instruments and from visible and near infrared to thermal infrared and microwave spectral sounders. One significant feature of these datasets is that they have multi-decades long climate records. Top-of-atmosphere (TOA) radiation and its anomaly represents the net heat energy input to the climate system, and the oceanic precipitation and ocean-land moisture transport can be used to quantify sea surface latent heat release. Based on the synergized datasets and the principle of global water and energy balances, global mean turbulent heat fluxes are estimated. The turbulent heat anomalies for the first two decades of the 21st century are, then, obtained mainly from global CERES radiation and GPCP precipitation anomalies, along with CERES derived ocean-land heat transports. Analysis indicates that the uncertainties in the estimated turbulent flux anomalies may be reduced considerably. The results suggest a strong needs in synchronized and synergized observations of atmospheric radiation, precipitation, and oceanic meteorological variables for long-term climate studies.

Bing Lin↗

Advanced Electric Propulsion for RLV Launched Geosynchronous Spacecraft

Solar Electric Propulsion (SEP) when used for station keeping and final orbit insertion has been shown to increase a geostationary satellite's payload when launched by existing expendable launch vehicles. In the case of reusable launch vehicles or expendable launch vehicles where an upper stage is an expensive option, this methodology can be modified by using the existing on-board apogee chemical system to perform a perigee burn and then letting the electric propulsion system complete the transfer to geostationary orbit. The elimination of upper stages using on-board chemical and electric propulsion systems was thus examined for GEO spacecraft. Launch vehicle step-down from an Atlas IIAR to a Delta 7920 (no upper stage) was achieved using expanded on-board chemical tanks, 40 kW payload power for electric propulsion, and a 60 day elliptical to GEO SEP orbit insertion. Optimal combined chemical and electric trajectories were found using SEPSPOT. While Hall and ion thrusters provided launch vehicle step-down and even more payload for longer insertion times, NH3 arcjets had insufficient performance to allow launch vehicle step-down. Degradation levels were only 5% to 7% for launch step-down cases using advanced solar arrays. Results were parameterized to allow comparisons for future reusable launch vehicles. Results showed that for an 8 W/kg initial power/launch mass power density spacecraft, 50% to 100% more payload can be launched using this method.

Oleson, Steven↗

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗

Finding the global minimum: a fuzzy end elimination implementation

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

NASA Program Exobiology↗

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Wheaton, Wm. A.↗

A ray tracing model of gravity wave propagation and breakdown in the middle atmosphere

Gravity wave ray tracing and wave packet theory is used to parameterize wave breaking in the mesosphere. Rays are tracked by solving the group velocity equations, and the interaction with the basic state is determined by considering the evolution of the packet wave action density. The ray tracing approach has a number of advantages over the steady state parameterization as the effects of gravity wave focussing and refraction, local dissipation, and wave response to rapid changes in the mean flow are more realistically considered; however, if steady state conditions prevail, the method gives identical results. The ray tracing algorithm is tested using both interactive and noninteractive models of the basic state. In the interactive model, gravity wave interaction with the polar night jet on a beta-plane is considered. The algorithm produces realistic polar night jet closure for weak topographic forcing of gravity waves. Planetary scale waves forced by local transfer of wave action into the basic flow in turn transfer their wave action into the zonal mean flow. Highly refracted rays are also found not to contribute greatly to the climatology of the mesosphere, as their wave action is severely reduced by dissipation during their lateral travel.

Schoeberl, M. R.↗