Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid modeling”

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 919 records · Page 51

Computational Analysis of Cryogenic Flow Through a Control Valve

The initial efforts to develop the capability to model valves used in rocket engine component testing at Stennis Space Center are documented. An axisymmetric model of a control valve with LN2 as the working fluid was developed. The goal was to predict the effect of change in the plug/sear region of the valve prior to testing. The valve flow coefficient was predicted for a range of plug positions. Verification of the calculations was carried out to quantify the uncertainty in the numerical answer. The modeled results compared well qualitatively to experimental trends. Additionally, insights into the flow processes in the valve were obtained. Benefits from the verification process included the ability to use coarser grids and insight into ways to reduce computational time by using double precision accuracy and non-integer grid ratios. Future valve modeling activities will include shape optimization of the valve/seat region and dynamic grid modeling.

Danes, Russell↗

A Model Calibration Method for Grid-Forming Inverters Using Iterative Bayesian Optimization

As inverter-based resources (IBRs) are rapidly deployed, especially at the distribution and microgrid levels, the need to include their accurate representations in power systems models increases. With well-calibrated models, IBR-interactions can be examined, preventing any stability and operational issues, especially in islanded systems. This paper proposes a method to calibrate generic inverter models using Bayesian optimization, but with a parameter-grouping approach to improve speed and accuracy. This approach is illustrated for a grid-forming inverter using synthetic data and field measurements from a simple IBRbased microgrid. Tests show that the calibration algorithm has modest computation cost and is robust to noisy measurements..

Biswas, Shuchismita↗

A Crop Yield Change Emulator for Use in GCAM and Similar Models: Persephone v1.0

Future changes in Earth system state will impact agricultural yields and, through these changed yields, can have profound impacts on the global economy. Global gridded crop models estimate the influence of these Earth system changes on future crop yields but are often too computationally intensive to dynamically couple into global multisector economic models, such as the Global Change Assessment Model (GCAM) and other similar-in-scale models. Yet, generalizing a faster site-specific crop model’s results to be used globally will introduce inaccuracies, and the question of which model to use is unclear given the wide variation in yield response across crop models. To examine the feedback loop among socioeconomics, Earth system changes, and crop yield changes, rapidly generated yield responses with some quantification of crop response uncertainty are desirable. The Persephone v1.0 response functions presented in this work are based on the Agricultural Model Intercomparison and Improvement Project (AgMIP) Coordinated Climate-Crop Modeling Project (C3MP) sensitivity test data set and are focused on providing GCAM and similar models with a tractable number of rapid to evaluate dynamic yield response functions corresponding to a range of the yield response sensitivities seen in the C3MP data set. With the Persephone response functions, a new variety of agricultural impact experiments will be open to GCAM and other economic models: for example, examining the economic impacts of a multi-year drought in a key agricultural region and how economic changes in response to the drought can, in turn, impact the drought.

Snyder, Abigail↗

Grid artifact reduction for direct digital radiography detectors based on rotated stationary grids with homomorphic filtering

Purpose: Grid artifacts are caused when using the antiscatter grid in obtaining digital x-ray images. In this paper, research on grid artifact reduction techniques is conducted especially for the direct detectors, which are based on amorphous selenium. Methods: In order to analyze and reduce the grid artifacts, the authors consider a multiplicative grid image model and propose a homomorphic filtering technique. For minimal damage due to filters, which are used to suppress the grid artifacts, rotated grids with respect to the sampling direction are employed, and min-max optimization problems for searching optimal grid frequencies and angles for given sampling frequencies are established. The authors then propose algorithms for the grid artifact reduction based on the band-stop filters as well as low-pass filters. Results: The proposed algorithms are experimentally tested for digital x-ray images, which are obtained from direct detectors with the rotated grids, and are compared with other algorithms. It is shown that the proposed algorithms can successfully reduce the grid artifacts for direct detectors. Conclusions: By employing the homomorphic filtering technique, the authors can considerably suppress the strong grid artifacts with relatively narrow-bandwidth filters compared to the normal filtering case. Using rotated grids also significantly reduces the ringing artifact. Furthermore, for specific grid frequencies and angles, the authors can use simple homomorphic low-pass filters in the spatial domain, and thus alleviate the grid artifacts with very low implementation complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗

Multi-Sensor Cloud and Aerosol Retrieval Simulator and Remote Sensing from Model Parameters : Aerosols - Part 2

The Multi-sensor Cloud Retrieval Simulator (MCRS) produces a simulated radiance product from any high-resolution general circulation model with interactive aerosol as if a specific sensor such as the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of the atmospheric column and land ocean surface at a specific location. Previously the MCRS code only included contributions from atmosphere and clouds in its radiance calculations and did not incorporate properties of aerosols. In this paper we added a new aerosol properties module to the MCRS code that allows users to insert a mixture of up to 15 different aerosol species in any of 36 vertical layers. This new MCRS code is now known as MCARS (Multi-sensor Cloud and Aerosol Retrieval Simulator). Inclusion of an aerosol module into MCARS not only allows for extensive, tightly controlled testing of various aspects of satellite operational cloud and aerosol properties retrieval algorithms, but also provides a platform for comparing cloud and aerosol models against satellite measurements. This kind of two-way platform can improve the efficacy of model parameterizations of measured satellite radiances, allowing the assessment of model skill consistently with the retrieval algorithm. The MCARS code provides dynamic controls for appearance of cloud and aerosol layers. Thereby detailed quantitative studies of the impacts of various atmospheric components can be controlled. In this paper we illustrate the operation of MCARS by deriving simulated radiances from various data field output by the Goddard Earth Observing System version 5 (GEOS-5) model. The model aerosol fields are prepared for translation to simulated radiance using the same model sub grid variability parameterizations as are used for cloud and atmospheric properties profiles, namely the ICA technique. After MCARS computes modeled sensor radiances equivalent to their observed counterparts, these radiances are presented as input to operational remote-sensing algorithms. Specifically, the MCARS-computed radiances are input into the processing chain used to produce the MODIS Data Collection 6 aerosol product (MOYD04). TheMOYD04 product is of course normally produced from MOYD021KM MODIS Level-1B radiance product directly acquired by the MODIS instrument. MCARS matches the format and metadata of a MOYD021KM product. The resulting MCARS output can be directly provided to MODAPS (MODIS Adaptive Processing System) as input to various operational atmospheric retrieval algorithms. Thus the operational algorithms can be tested directly without needing to make any software changes to accommodate an alternative input source. We show direct application of this synthetic product in analysis of the performance of the MOD04 operational algorithm. We use biomass-burning case studies over Amazonia employed in a recent Working Group on Numerical Experimentation (WGNE)-sponsored study of aerosol impacts on numerical weather prediction (Freitas et al., 2015). We demonstrate that a known low bias in retrieved MODIS aerosol optical depth appears to be due to a disconnect between actual column relative humidity and the value assumed by the MODIS aerosol product.

aerosol retrieval↗

Self-Secure Inverters Against Malicious Setpoints

The next generation of grid-interactive inverters brings a communication feature that allows data sharing from utility supervisory controllers and smart devices that are connected to the same network. This feature enhances the control capabilities of grid-interactive inverters to provide services beyond active power injection. However, communication networks entail more vulnerable surfaces to malicious attacks that may result in modifying active and reactive power setpoints and causing weak-grid conditions or abnormal inverter operation. In this paper, steady-state and the dynamic behavior of the inverter for the incoming setpoints are analyzed to detect false data injection attacks and provide device-level security. The steady-state behavior of the inverter in the operating region is determined from the grid parameters such as the grid voltage and the grid impedance. These estimations are accomplished by the proposed self-security technique through a low-frequency signal injection-based approach combined with the recursive least square method. Moreover, a reduced fourth-order inverter model is used as the dynamic reference model, and grid parameters as well as the incoming setpoints are implemented to the reference model to verify whether the dynamic behavior of the inverter is inside the permissible region of operation. The validity and performance of the proposed method are verified experimentally through Allen-Bradley Powerflex 755 three-phase inverter and a 12 kW NHR 9410 regenerative power grid emulator. The results show that the self-secure smart-inverter is able to accept or reject the incoming commands and thus is protected from malicious cyber-physical attacks.

Hossen, Tareq↗

A crop yield change emulator for use in GCAM and similar models: Persephone v1.0

Abstract. Future changes in Earth system state will impact agricultural yields and, through these changed yields, can have profound impacts on the global economy. Global gridded crop models estimate the influence of these Earth system changes on future crop yields but are often too computationally intensive to dynamically couple into global multi-sector economic models, such as the Global Change Assessment Model (GCAM) and other similar-in-scale models. Yet, generalizing a faster site-specific crop model's results to be used globally will introduce inaccuracies, and the question of which model to use is unclear given the wide variation in yield response across crop models. To examine the feedback loop among socioeconomics, Earth system changes, and crop yield changes, rapidly generated yield responses with some quantification of crop response uncertainty are desirable. The Persephone v1.0 response functions presented in this work are based on the Agricultural Model Intercomparison and Improvement Project (AgMIP) Coordinated Climate-Crop Modeling Project (C3MP) sensitivity test data set and are focused on providing GCAM and similar models with a tractable number of rapid to evaluate dynamic yield response functions corresponding to a range of the yield response sensitivities seen in the C3MP data set. With the Persephone response functions, a new variety of agricultural impact experiments will be open to GCAM and other economic models: for example, examining the economic impacts of a multi-year drought in a key agricultural region and how economic changes in response to the drought can, in turn, impact the drought.

58 GEOSCIENCES↗

Seismic Tomographic Modeling of the Crust and Upper Mantle beneath Israel and the Middle East: Improved Resolution through Optimized Model Parameterization

Accurate regional seismic travel-time (RSTT) predictions rely on regional phases (e.g., Pg, Lg, Pn, Sn ) to account for 3D effects in the crust and upper mantle that are not captured by 1D models traditionally used for real-time location. The RSTT prediction model accounts for regional-scale crust and upper mantle structure globally by incorporating regional seismic phases into its travel-time calculations. Previous versions of the RSTT model have used a constant grid cell size of 1°. To improve the tomographic accuracy of recovering velocity structure at regional scales, we, in this study, perform data-driven grid refinement on the RSTT model down to a 0.125° grid (~14 km) in pursuit of two main goals: (1) to test the limits of RSTT capability and accuracy of determined velocity structure through variable grid refinement and (2) to image smaller structures in Israel and the Middle East and illuminate upper mantle dynamics operating in this complex tectonic area. We investigate the effects of model parameterization as grid cell size decreases and the trade-offs between recovered velocity structures. Our final dataset includes 4751 events and 499 stations that recorded 79,344 Pn and 7489 Pg . The variable grid refinement method allows recovery of finer-scale velocity structures and reduces travel-time residuals in areas with the highest data coverage. At smaller grid cell sizes, longer paths need to be upweighted to stabilize the inversion. Results illuminate tectonic features undefined in coarser grid-size models; in particular, we observe mantle perturbations related to the subduction zone around the Cyprian arc and crustal anomalies near the Dead Sea fault and throughout the Anatolian plate.

58 GEOSCIENCES↗

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the atmosphere–land–ocean–sea ice interface processes in the Regional Arctic System Model version 1 (RASM1) using local and globally gridded observations

The Regional Arctic System Model version 1 (RASM1) has been developed to provide high-resolution simulations of the Arctic atmosphere–ocean–sea ice–land system. Here, we provide a baseline for the capability of RASM to simulate interface processes by comparing retrospective simulations from RASM1 for 1990–2014 with the Community Earth System Model version 1 (CESM1) and the spread across three recent reanalyses. Evaluations of surface and 2 m air temperature, surface radiative and turbulent fluxes, precipitation, and snow depth in the various models and reanalyses are performed using global and regional datasets and a variety of in situ datasets, including flux towers over land, ship cruises over oceans, and a field experiment over sea ice. These evaluations reveal that RASM1 simulates precipitation that is similar to CESM1, reanalyses, and satellite gauge combined precipitation datasets over all river basins within the RASM domain. Snow depth in RASM is closer to upscaled surface observations over a flatter region than in more mountainous terrain in Alaska. The sea ice–atmosphere interface is well simulated in regards to radiation fluxes, which generally fall within observational uncertainty. RASM1 monthly mean surface temperature and radiation biases are shown to be due to biases in the simulated mean diurnal cycle. At some locations, a minimal monthly mean bias is shown to be due to the compensation of roughly equal but opposite biases between daytime and nighttime, whereas this is not the case at locations where the monthly mean bias is higher in magnitude. These biases are derived from errors in the diurnal cycle of the energy balance (radiative and turbulent flux) components. Therefore, the key to advancing the simulation of SAT and the surface energy budget would be to improve the representation of the diurnal cycle of radiative and turbulent fluxes. The development of RASM2 aims to address these biases. Still, an advantage of RASM1 is that it captures the interannual and interdecadal variability in the climate of the Arctic region, which global models like CESM cannot do.

58 GEOSCIENCES↗

A High-resolution Regional Wave Resource Characterization For The U.S. West Coast

Objectives/Scope: Wave resource characterization is a critical step for wave energy converter deployment in the coastal ocean and relies on long-term, high-resolution wave datasets. This study presents a detailed modeling study of the wave resource along the U.S. West Coast (Washington, Oregon, and California), a coastal region that was identified with high wave energy potential in earlier studies. Methods, Procedures, Process: The wave hindcast covers a 32-year period from 1979 to 2010 and is based on a multi-resolution, unstructured-grid SWAN model framework. Model configuration closely follows and meets the requirements recommended by the International Electrotechnical Commission Technical Specification (IEC TS) for wave energy resource assessment and characterization (Class 2 - feasibility study). The model domain covers the entire U.S. Exclusive Economic Zone (EEZ) in the West Coast and has a spatial resolution varying from ~300 m in the nearshore region (20 km from the shoreline) to ~2500 m within the EEZ and ~5000 m at the open boundary, which extends beyond the EEZ. The model was forced by hourly 2-D wave spectra produced by a two-way nested WaveWatch III model, which covers the global ocean domain and the broader U.S. West Coast region domain with spatial resolutions of 0.5 degree and 10 arc-minutes, respectively. Both wave models are forced by hourly, 0.5-degree wind forcing obtained from NCEP’s Climate Forecast System Reanalysis (CFSR) product. Results, Observations, Conclusions: The standard model output for the SWAN model includes 3-hourly output for the six IEC wave resource parameters (e.g., omnidirectional wave power) at each grid point and hourly 2-D spectra at more than 50 NDBC buoys. Extensive model validation was achieved by comparing the six model-predicted IEC parameters with those derived from field observations at representative NDBC buoys. The error statistics indicated the model’s satisfactory performance. Further analyses were conducted to systematically evaluate the temporal and spatial distributions of wave energy potential and wave climate along the U.S. West Coast. Results suggest that Washington and Oregon coasts have similar nearshore wave resource, which is significantly higher than resources in Southern California. Strong seasonal variations are also observed, e.g., high wave energy tends to occur in the winter months. In summary, this study produced the first high-resolution, comprehensive dataset on wave energy distribution along the U.S. West Coast. Novel/Additive Information: The results are being used by the National Renewable Energy Laboratory to update the MHK Atlas, which was originally derived from NOAA’s 4-arc-minute WaveWatch III model output. In addition, the monthly averaged wave energy climatology dataset can be readily shared to support a variety of research and application efforts within the EEZ of the U.S. West Coast.

Wang, Taiping↗

Spectral Irradiance Calibration in the Infrared: New Composite Spectra, Comparison with Model Atmospheres, and Far-Infrared Extrapolations - Part 7

We present five new absolutely calibrated continuous stellar spectra constructed as far as possible from spectral fragments observed from the ground, the Kuiper Airborne Observatory (KAO), and the IRAS Low Resolution Spectrometer. These stars-alpha Boo, gamma Dra, alpha Cet, gamma Cru, and mu UMa-augment our six, published, absolutely calibrated spectra of K and early-M giants. All spectra have a common calibration pedigree. A revised composite for alpha Boo has been constructed from higher quality spectral fragments than our previously published one. The spectrum of gamma Dra was created in direct response to the needs of instruments aboard the Infrared Space Observatory (ISO); this star's location near the north ecliptic pole renders it highly visible throughout the mission. We compare all our low-resolution composite spectra with Kurucz model atmospheres and find good agreement in shape, with the obvious exception of the SiO fundamental, still lacking in current grids of model atmospheres. The CO fundamental seems slightly too deep in these models, but this could reflect our use of generic models with solar metal abundances rather than models specific to the metallicities of the individual stars. Angular diameters derived from these spectra and models are in excellent agreement with the best observed diameters. The ratio of our adopted Sirius and Vega models is vindicated by spectral observations. We compare IRAS fluxes predicted from our cool stellar spectra with those observed and conclude that, at 12 and 25 microns, flux densities measured by IRAS should be revised downwards by about 4.1% and 5.7%, respectively, for consistency with our absolute calibration. We have provided extrapolated continuum versions of these spectra to 300 microns, in direct support of ISO (PHT and LWS instruments). These spectra are consistent with IRAS flux densities at 60 and 100 microns.

Cohen, Martin↗

Model Stellar Atmospheres and Real Stellar Atmospheres and Status of the ATLAS12 Opacity Sampling Program and of New Programs for Rosseland and for Distribution Function Opacity

I discuss errors in theory and in interpreting observations that are produced by the failure to consider resolution in space, time, and energy. I discuss convection in stellar model atmospheres and in stars. Large errors in abundances are possible such as the factor of ten error in the Li abundance for extreme Population II stars. Finally I discuss the variation of microturbulent velocity with depth, effective temperature, gravity, and abundance. These variations must be dealt with in computing models and grids and in any type of photometric calibration. I have also developed a new opacity-sampling version of my model atmosphere program called ATLAS12. It recognizes more than 1000 atomic and molecular species, each in up to 10 isotopic forms. It can treat all ions of the elements up through Zn and the first 5 ions of heavier elements up through Es. The elemental and isotopic abundances are treated as variables with depth. The fluxes predicted by ATLAS12 are not accurate in intermediate or narrow bandpass intervals because the sample size is too small. A special stripped version of the spectrum synthesis program SYNTHE is used to generate the surface flux for the converged model using the line data on CD-ROMs 1 and 15. ATLAS12 can be used to produce improved models for Am and Ap stars. It should be very useful for investigating diffusion effects in atmospheres. It can be used to model exciting stars for H II regions with abundances consistent with those of the H II region. These programs and line files will be distributed on CD-ROMs.

Kurucz, Robert L.↗

Electric Propulsion Interactions Code (EPIC): Recent Enhancements and Goals for Future Capabilities

The Electric Propulsion Interactions Code (EPIC) is the leading interactive computer tool for assessing the effects of electric thruster plumes on spacecraft subsystems. EPIC, developed by SAIC under the sponsorship of the Space Environments and Effects (SEE) Program at the NASA Marshall Space Flight Center, has three primary modules. One is PlumeTool, which calculates plumes of electrostatic thrusters and Hall-effect thrusters by modeling the primary ion beam as well as elastic scattering and charge-exchange of beam ions with thruster-generated neutrals. ObjectToolkit is a 3-D object definition and spacecraft surface modeling tool developed for use with several SEE Program codes. The main EPIC interface integrates the thruster plume into the 3-D geometry of the spacecraft and calculates interactions and effects of the plume with the spacecraft. Effects modeled include erosion of surfaces due to sputtering, re-deposition of sputtered materials, surface heating, torque on the spacecraft, and changes in surface properties due to erosion and deposition. In support of Prometheus I (JIMO), a number of new capabilities and enhancements were made to existing EPIC models. Enhancements to EPIC include adding the ability to scale and view individual plume components, to import a neutral plume associated with a thruster (to model a grid erosion plume, for example), and to calculate the plume from new initial beam conditions. Unfortunately, changes in program direction have left a number of desired enhancements undone. Variable gridding over a surface and resputtering of deposited materials, including multiple bounces and sticking coefficients, would significantly enhance the erosion/deposition model. Other modifications such as improving the heating model and the PlumeTool neutral plume model, enabling time dependent surface interactions, and including EM1 and optical effects would enable EPIC to better serve the aerospace engineer and electric propulsion systems integrator. We review EPIC S overall capabilities and recent modifications, and discuss directions for future enhancements.

Gardner, Barbara M.↗

Criticality Accident Alarm System Modeling for the Uranium Processing Facility

The Criticality Accident Alarm System (CAAS) for the Uranium Processing Facility (UPF) provides coverage for two primary processing facilities. For UPF, detailed documentation and drawings of the design are readily available. While this design information allows for some CAAS modeling features to match the actual design with a high level of accuracy, assumptions and simplifications are necessary as it is unrealistic and unnecessary to replicate many aspects of the design. A set of modeling guidelines and techniques were used to develop consistency between the MCNP6 models. While ensuring conservatism, overly conservative choices must be avoided to demonstrate the minimum accident of concern can be detected. To ensure proper modeling choices were made for the final models, a group of parametric studies were performed to evaluate the impact on the dose rates detected at the individual detectors. Some of these modeling choices and/or associated parametric studies include: - Modeling the grid steel area above the detectors as homogeneous stainless steel and determining the impact on detectability with an increasing density. - Modeling the presence of equipment in close proximity to the CAAS detectors. - Determining an appropriate material density to represent miscellaneous items not specifically modeled. - Evaluating the impact of sprinkler activation. The results of UPF CAAS analyses show that the modeling assumptions made based on the detailed design and parametric studies are appropriate and demonstrate accurate yet conservative models have been developed for detecting minimum accident of concern within the UPF.

61 RADIATION PROTECTION AND DOSIMETRY↗

Spectral properties of cool stars: Extended abundance analysis of 1,617 planet-search stars

We present a catalog of uniformly determined stellar properties and abundances for 1,617 F, G, and K stars using an automated spectral synthesis modeling procedure. All stars were observed using the HIRES spectrograph at Keck Observatory. Our procedure used a single line list to fit model spectra to observations of all stars to determine effective temperature, surface gravity, metallicity, projected rotational velocity, and the abundances of 15 elements (C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, and Y). Sixty percent of the sample had Hipparcos parallaxes and V-band photometry, which we combined with the spectroscopic results to obtain mass, radius, and luminosity. Additionally, we used the luminosity, effective temperature, metallicity and α-element enhancement to interpolate in the Yonsei–Yale isochrones to derive mass, radius, gravity, and age ranges for those stars. Finally, we determined new relations between effective temperature and macroturbulence for dwarfs and subgiants. Our analysis achieved precisions of 25 K in T{sub eff}, 0.01 dex in [M/H], 0.028 dex for logg, and 0.5 km s{sup −1} in vsin i based on multiple observations of the same stars. The abundance results were similarly precise, between ∼0.01 and ∼0.04 dex, though trends with respect to T{sub eff} remained for which we derived empirical corrections. The trends, though small, were much larger than our uncertainties and are shared with published abundances. We show that changing our model atmosphere grid accounts for most of the trend in [M/H] between 5000 and 5500 K, indicating a possible problem with the atmosphere models or opacities.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of Near-Surface Models in Large-Eddy Simulations of a Neutrally Stratified Atmospheric Boundary Layer

Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Senocak, I.↗