Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 235 records · Page 13

The development of a population of 4D pediatric XCAT phantoms for imaging research and optimization

Purpose: We previously developed a set of highly detailed 4D reference pediatric extended cardiac-torso (XCAT) phantoms at ages of newborn, 1, 5, 10, and 15 yr with organ and tissue masses matched to ICRP Publication 89 values. In this work, we extended this reference set to a series of 64 pediatric phantoms of varying age and height and body mass percentiles representative of the public at large. The models will provide a library of pediatric phantoms for optimizing pediatric imaging protocols. Methods: High resolution positron emission tomography-computed tomography data obtained from the Duke University database were reviewed by a practicing experienced radiologist for anatomic regularity. The CT portion of the data was then segmented with manual and semiautomatic methods to form a target model defined using nonuniform rational B-spline surfaces. A multichannel large deformation diffeomorphic metric mapping algorithm was used to calculate the transform from the best age matching pediatric XCAT reference phantom to the patient target. The transform was used to complete the target, filling in the nonsegmented structures and defining models for the cardiac and respiratory motions. The complete phantoms, consisting of thousands of structures, were then manually inspected for anatomical accuracy. The mass for each major tissue was calculated and compared to linearly interpolated ICRP values for different ages. Results: Sixty four new pediatric phantoms were created in this manner. Each model contains the same level of detail as the original XCAT reference phantoms and also includes parameterized models for the cardiac and respiratory motions. For the phantoms that were 10 yr old and younger, we included both sets of reproductive organs. This gave them the capability to simulate both male and female anatomy. With this, the population can be expanded to 92. Wide anatomical variation was clearly seen amongst the phantom models, both in organ shape and size, even for models of the same age and sex. The phantoms can be combined with existing simulation packages to generate realistic pediatric imaging data from different modalities. Conclusions: This work provides a large cohort of highly detailed pediatric phantoms with 4D capabilities of varying age, height, and body mass. The population of phantoms will provide a vital tool with which to optimize 3D and 4D pediatric imaging devices and techniques in terms of image quality and radiation-absorbed dose.

60 APPLIED LIFE SCIENCES↗

Extrapolation of Type Ia Supernova Spectra into the Near-infrared Using Principal Component Analysis

Abstract We present a method of extrapolating the spectroscopic behavior of Type Ia supernovae (SNe Ia) in the near-infrared (NIR) wavelength regime up to 2.30 μ m using optical spectroscopy. Such a process is useful for accurately estimating K-corrections and other photometric quantities of SNe Ia in the NIR. A principal component analysis is performed on data consisting of Carnegie Supernova Project I & II optical and NIR FIRE spectra to produce models capable of making these extrapolations. This method differs from previous spectral template methods by not parameterizing models strictly by photometric light-curve properties of SNe Ia, allowing for more flexibility of the resulting extrapolated NIR flux. A difference of around −3.1% to −2.7% in the total integrated NIR flux between these extrapolations and the observations is seen here for most test cases including Branch core-normal and shallow-silicon subtypes. However, larger deviations from the observation are found for other tests, likely due to the limited high-velocity and broad-line SNe Ia in the training sample. Maximum-light principal components are shown to allow for spectroscopic predictions of the color-stretch light-curve parameter, s BV , within approximately ±0.1 units of the value measured with photometry. We also show these results compare well with NIR templates, although in most cases the templates are marginally more fitting to observations, illustrating a need for more concurrent optical+NIR spectroscopic observations to truly understand the diversity of SNe Ia in the NIR.

Astronomy & Astrophysics↗

Interactions Between Clouds and Wind-Driven Surface Heat Exchanges over Land

Earth system model experiments show that increasing horizontal resolution fundamentally alters the simulated soil-moisture-precipitation feedback. Kilometer-scale simulations often produce weaker or even negative feedback compared to coarse-resolution models. A key difference of kilometer-scale models is that they resolve mesoscale secondary circulations, including boundary layer horizontal rolls and cellular structures, in addition to cold pools and downdrafts associated with convective precipitation. However, because the relevant processes occur on yet-smaller scales, these circulations are often poorly resolved. This project demonstrated that boundary layer secondary circulations significantly affect surface heat exchanges and wind gusts, and that current model parameterizations can misrepresent these processes at kilometer-scale resolution. Using DOE Atmospheric Radiation Measurement (ARM) observations and targeted experiments with the DOE Energy Exascale Earth System Model (E3SM), we identified physically unrealistic wind gust and surface flux responses to secondary circulations, diagnosed a systematic overestimation of wind shear in convective cold pools, and uncovered a multivariate relationship between land surface fluxes and the scales of updrafts that form shallow cumulus clouds. These findings provide observation-based recommendations for improving parameterizations of surface fluxes and wind gusts in high-resolution Earth system models, thereby reducing uncertainty in convective storm prediction and land-atmosphere feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Adaptive Grid Lens Modeling of the Cosmic Horseshoe Using Hubble Space Telescope Imaging

The Cosmic Horseshoe gravitational lensing system is an extraordinary example of strong gravitational lensing both due to the nearly complete Einstein ring formed in this system and due to the star-forming nature of the lensed $z$ = 2.38 Lyman-break galaxy. In this paper, we describe the development of a new lens modeling package and the lens models produced using the Hubble Space Telescope imaging. Our new lens modeling package uses adaptive grid methods, based especially on the algorithm described by Vegetti & Koopmans. The new lens modeling package introduces the K-means method to deal with multiple background sources. We utilize two parameterized models, the singular isothermal ellipsoid and the Navarro–Frenk–White model in order to optimize the Bayesian penalty function. Color–color diagrams show two distinct colors exist in the Einstein ring, which suggests that the Cosmic Horseshoe is formed from two background galaxies or from a highly irregular galaxy. We find that the best lens model includes two components and the total mass within the Einstein ring is $(5.43\pm 0.08)\times {10}^{12}\,{M}_{\odot }$. Furthermore, the background source reconstruction shows two peaks in the source plane and the most probable lens profile is consistent with the light profile of the foreground galaxy when comparing axis ratios and position angles.

79 ASTRONOMY AND ASTROPHYSICS↗

Global electric field determination in the Earth's outer magnetosphere using energetic charged particles

Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental measurement (for zeroth order MHD equilibrium) has been made infrequently and with poor spatial coverage - the global electric field. This oversight is due in part to the neglect of theorists. However, there is renewed interest in the convection electric field because it is now realized to be central to many magnetospheric processes, including the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere ionosphere coupling, ring current and radiation belt transport, substorm injections, and several acceleration mechanisms. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models) and we are left with an overly simplistic theoretical field, the Volland-Stern electric field model. Single point measurements of the plasmapause were used to infer the appropriate amplitudes of this model, parameterized by K(sub p). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 20 years. The analysis of current data sets requires a great deal more accuracy than can be provided by the Volland-Stern model. The variability of electric field shielding has not been properly addressed although effects of penetrating magnetospheric electric fields has been seen in mid-and low-latitude ionospheric data sets. The growing interest in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus we proposed and developed algorithms for extracting electric fields from particle data taken in the Earth's magnetosphere. As a test of the effectiveness of these new techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit from 1984 to 1989.

Eastman, Timothy E.↗

Implementation of a generalized actuator line model for wind turbine parameterization in the Weather Research and Forecasting model

A generalized actuator line (GAL) wind turbine parameterization is implemented within the Weather Research and Forecasting model to enable high-fidelity large-eddy simulations of wind turbine interactions with boundary layer flows under realistic atmospheric forcing conditions. Numerical simulations using the GAL parameterization are evaluated against both an already implemented generalized actuator disk (GAD) wind turbine parameterization and two field campaigns that measured the inflow and near-wake regions of a single turbine. The representation of wake wind speed, variance, and vorticity distributions is examined by comparing fine-resolution GAL and GAD simulations and GAD simulations at both fine and coarse-resolutions. The higher-resolution simulations show slightly larger and more persistent velocity deficits in the wake and substantially increased variance and vorticity when compared to the coarse-resolution GAD. The GAL generates distinct tip and root vortices that maintain coherence as helical tubes for approximately one rotor diameter downstream. Coarse-resolution simulations using the GAD produce similar aggregated wake characteristics to both fine-scale GAD and GAL simulations at a fraction of the computational cost. Furthermore, the GAL parameterization provides the capability to resolve near wake physics, including vorticity shedding and wake expansion.

17 WIND ENERGY↗

A spectrum model for weakly anisotropic turbulence

A simple model, parameterized by the Reynolds stress anisotropy, is proposed for the spectrum of weakly anisotropic turbulence. It contains a model constant that affects its region of realizability. This spectrum model is used to derive a one-point closure to the rapid pressure-strain term. The derived pressure-strain model is linear in the Reynolds stress anisotropy and is of the same form as the closure model of Launder et al. (1975). The spectrum model becomes unrealizable in some regions of wave space for sufficiently large anisotropy of the Reynolds stress, and this is used to infer the region of validity of the linear closure model. It is found that the extent of the valid region is very small when the model constant is set to match rapid distortion theory, and largest for a model-constant set close to the value suggested by LRR. However, even the largest valid domain does not extend very far from isotropy, suggesting inherent weakness in the linear pressure-strain models.

Shih, Tsan-Hsing↗

Quantifying the impacts of land surface schemes and dynamic vegetation on the model dependency of projected changes in surface energy and water budgets

Assessing and quantifying the uncertainties in projected future changes of energy and water budgets over land surface are important steps toward improving our confidence in climate change projections. In our study, the contribution of land surface models to the inter-GCM variation of projected future changes in land surface energy and water fluxes are assessed based on output from 19 global climate models (GCMs) and offline Community Land Model version 4 (CLM4) simulations driven by meteorological forcing from the 19 GCMs. Similar offline simulations using CLM4 with its dynamic vegetation submodel are also conducted to investigate how dynamic vegetation feedback, a process that is being added to more earth system models, may amplify or moderate the intermodel variations of projected future changes. Projected changes are quantified as the difference between the 2081–2100 period from the Representative Concentration Pathway 8.5 (RCP8.5) future experiment and the 1981–2000 period from the historical simulation. Under RCP8.5, projected changes in surface water and heat fluxes show a high degree of model dependency across the globe. Although precipitation is very likely to increase in the high latitudes of the Northern Hemisphere, a high degree of model-related uncertainty exists for evapotranspiration, soil water content, and surface runoff, suggesting discrepancy among land surface models (LSMs) in simulating the surface hydrological processes and snow-related processes. Large model-related uncertainties for the surface water budget also exist in the Tropics including southeastern South America and Central Africa. Moreover, these uncertainties would be reduced in the hypothetical scenario of a single near-perfect land surface model being used across all GCMs, suggesting the potential to reduce uncertainties through the use of more consistent approaches toward land surface model development. Under such a scenario, the most significant reduction is likely to be seen in the Northern Hemisphere high latitudes. Including representation of vegetation dynamics is expected to further amplify the model-related uncertainties in projected future changes in surface water and heat fluxes as well as soil moisture content. This is especially the case in the high latitudes of the Northern Hemisphere (e.g., northwestern North America and central North Asia) where the projected vegetation changes are uncertain and in the Tropics (e.g., the Amazon and Congo Basins) where dense vegetation exists. Finally, findings from this study highlight the importance of improving land surface model parameterizations related to soil and snow processes, as well as the importance of improving the accuracy of dynamic vegetation models.

58 GEOSCIENCES↗

Global electric field determination in the Earth's outer magnetosphere using charged particles

Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental (for a zeroeth order magnetohydrodynamic (MHD) equilibrium) measurement was made infrequently and with poor spatial coverage: the global electric field. This oversight is in part due to the difficulty of measuring a plasma electric field, and in part due to the difficulty of measuring a plasma electric field, and in part due to the neglect of theorists. However, there is renewed interest in the convection electric field, since it has been realized that it is vital for understanding many aspects of the magnetosphere: the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere-ionosphere coupling, ring current and radiation belt transport, substorm injections, acceleration mechanisms, etc. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models), and we are left with an overly simplistic theoretical field, the Volland-Stern electric field mode. Again, single point measurements of the plasma pause were used to infer the appropriate amplitudes of the model, parameterized by Kp (Maynard & Chen, JGR 1975). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 15 years. However, the data sets being taken today require a great deal more accuracy than can be provided by the Volland-Stern model. Nor has the variability of the electric field shielding been properly addressed, although effects of penetrating magnetospheric electric fields has been seen in mid- and low-latitude ionospheric data sets. The growing interests in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus, we proposed and developed algorithms for extracting electric fields from particle data taken in the earth's magnetosphere. As a test of the effectiveness of these techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit between 1984-1988. Some analytic tools had to be developed before construction of computer algorithms, and they are discussed.

Eastman, T.↗

Comparison of Daytime Low-Level Cloud Properties Derived From GOES and ARM SGP Measurements

Large-scale satellite data are critical for both verifying and improving general circulation model parameterizations of clouds and radiation for climate prediction. For reliable application of satellite data sets in cloud processes and climate models, it is important to have a reasonable estimate of the errors in the derived cloud properties. The daytime single-layered low-level cloud properties retrieved by the Geostationary Operational Environmental Satellite system (GOES) are compared with ground-based observations and retrievals over the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Central Facility from June 1998 to December 2006. The GOES retrievals are made via the Visible-Infrared Solar-infrared Split-window Technique. They are spatially averaged within a 0.15° × 0.15° box centered on the ARM SGP site, and the ARM surface observations are temporally averaged ±15 min around the GOES scans to produce collocated pairs. Comparisons are made for monthly means, diurnal means, and one-to-one GOES and ARM collocated pairs. GOES T eff is highly correlated with ARM T top cloud temperature, having an R 2 value of 0.75, though GOES exhibits a cold bias. GOES-retrieved τ and liquid water path have very good agreement with ARM retrievals with R 2 s of 0.45 and 0.47, while r e (GOES), on average, is about 2 μm greater than ARM r e . An examination of solar and viewing geometry has shown that GOES-retrieved mean r e and τ values are impacted by solar zenith angle and especially scattering angle, which is not unexpected and needs to be accounted for by users.

54 ENVIRONMENTAL SCIENCES↗

The seasonal CO2 cycle on Mars - An application of an energy balance climate model

Energy balance climate models of the Budyko-Sellers variety are applied to the carbon-dioxide cycle on Mars. Recent data available from the Viking mission, in particular the seasonal pressure variations measured by Viking landers, are used to constrain the models. No set of parameters was found for which a one-dimensional model parameterized in terms of ground temperature gave an adequate fit to the observed pressure variations. A modified, two-dimensional model including the effects of dust storms and the polar hood reasonably reproduces the pressure curve, however. The implications of these results for Martian climate changes are discussed.

James, P. B.↗

Satellite Derived Land Surface Temperature for Model Assimilation

Studies have shown that land surface temperature (LST) tendencies are sensitive to the surface moisture availability which is a function of soil moisture and vegetation. The assimilation of satellite derived LST tendencies into the surface energy budget of mesoscale models has shown promise in improving the representation of the complex effects of both soil moisture and vegetation within the models for short term simulations. LST derived from geostationary satellites has the potential of providing the temporal and spatial resolution needed for an LST assimilation process. This paper presents an analysis comparing the LST derived from GOES-8 infrared measurements with LST calculated by the MM5 numerical model. The satellite derived LSTs are calculated using a physical split window approach using channels 4 and 5 of GOES-8. The differences in the LST data sets, especially the tendencies, are presented and examined. Quantifying the differences between the data sets provide insight of possible weaknesses in the model parameterizations affecting the surface energy budget calculations and an indication of the potential effectiveness o f assimilating LST into the models.

Suggs, Ronnie J.↗

CMIP5 Climate Models Overestimate Cooling by Volcanic Aerosols

We compare projections of the observed hemispherical mean surface temperature (HadCRUT4.6.0.0) and the ensemble mean of CMIP5 climate models' simulations on a set of standard regression model forcing variables. We find that the volcanic aerosol regression coefficients of the CMIP5 simulations are consistently significantly larger (by 40–49%) than the volcanic aerosol coefficients of the observed temperature. The probability that the observed differences are caused just by chance is much less than 0.01. The overestimate is due to the climate models' response to volcanic aerosol radiative forcing. The largest overestimate occurs in the winter season of each hemisphere. We hypothesize that the models' parameterization of aerosol-cloud interactions within ice and mixed phase clouds is a likely source of this discrepancy. Furthermore, the models significantly underestimate the effect of solar variability on temperature for both hemispheres.

54 ENVIRONMENTAL SCIENCES↗

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗

GeoTess: A generalized Earth model software utility

GeoTess is a model parameterization and software support library that manages the construction, population, storage, and interrogation of data stored in 2D and 3D Earth models. Here, the software is available in Java and C++, with a C interface to the C++ library.

58 GEOSCIENCES↗

Simulation Studies of the Effect of Forest Spatial Structure on InSAR Signature

The height of scattering phase retrieved from InSAR data is considered being correlated with the tree height and the spatial structure of the forest stand. Though some researchers have used simple backscattering models to estimate tree height from the height of scattering center, the effect of forest spatial structure on InSAR data is not well understood yet. A three-dimensional coherent radar backscattering model for forest canopies based on realistic three-dimensional scene was used to investigate the effect in this paper. The realistic spatial structure of forest canopies was established either by field measurements (stem map) or through use of forest growth model. Field measurements or a forest growth model parameterized using local environmental parameters provides information of forest species composition and tree sizes in certain growth phases. A fractal tree model (L-system) was used to simulate individual 3- D tree structure of different ages or heights. Trees were positioned in a stand in certain patterns resulting in a 3-D medium of discrete scatterers. The radar coherent backscatter model took the 3-D forest scene as input and simulates the coherent radar backscattering signature. Interferometric SAR images of 3D scenes were simulated and heights of scattering phase centers were estimated from the simulated InSAR data. The effects of tree height, crown cover, crown depth, and the spatial distribution patterns of trees on the scattering phase center were analyzed. The results will be presented in the paper.

Sun, Guoqing↗

The Effect of Satellite Observing System Changes on MERRA Water and Energy Fluxes

Because reanalysis data sets offer state variables and fluxes at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Although one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics, another issue, perhaps less appreciated, is the fact that continual but discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (approximately the last 30 years). In particular we highlight the effect on the reanalysis of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We first examine MERRA fluxes from the perspective of how physical modes of variability (e.g. ENSO events, Pacific Decadal Variability) are contaminated by artificial step-like trends induced by the onset of new moisture data these two satellite observing systems. Secondly, we show how Redundancy Analysis, a statistical regression methodology, is effective in relating these artifact signals in the moisture and temperature analysis increments to their presence in the physical flux terms (e.g. precipitation, radiation). This procedure is shown to be effective greatly reducing the artificial trends in the flux quantities.

Robertson, Franklin R.↗

The Effect of Satellite Observing System Changes on MERRA Water and Energy Fluxes

Because reanalysis data sets offer state variables and fluxes at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Although one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics, another issue, perhaps less appreciated, is the fact that continual but discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (approx. the last 30 years). In particular we highlight the effect on the reanalysis of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 and, more prominently, with improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We first examine MERRA fluxes from the perspective of how physical modes of variability (e.g. ENSO events, Pacific Decadal Variability) are contained by artificial step-like trends induced by the onset of new moisture data these two satellite observing systems. Secondly, we show how Redundancy Analysis, a statistical regression methodology, is effective in relating these artifact signals in the moisture and temperature analysis increments to their presence in the physical flux terms (e.g. precipitation, radiation). This procedure is shown to be effective greatly reducing the artificial trends in the flux quantities.

Robertson, Franklin R.↗