Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 1,585 records · Page 88

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗

Asteroseismic modelling of solar-type stars: a deeper look at the treatment of initial helium abundance

ABSTRACT Detailed understanding of stellar physics is essential towards a robust determination of stellar properties (e.g. radius, mass, and age). Among the vital input physics used in the modelling of solar-type stars which remain poorly constrained, is the initial helium abundance. To this end, when constructing stellar model grids, the initial helium abundance is estimated either (i) by using the semi-empirical helium-to-heavy element enrichment ratio, (ΔY/ΔZ), anchored to the standard big bang nucleosynthesis value, or (ii) by setting the initial helium abundance as a free variable. Adopting 35 low-mass, solar-type stars with multiyear Kepler photometry from the asteroseismic ‘LEGACY’ sample, we explore the systematic uncertainties on the inferred stellar parameters (i.e. radius, mass, and age) arising from the treatment of the initial helium abundance in stellar model grids. The stellar masses and radii derived from grids with free initial helium abundance are lower compared to those from grids based on a fixed ΔY/ΔZ ratio. We find the systematic uncertainties on mean density, radius, mass, and age arising from grids which employ a fixed value of ΔY/ΔZ and those with free initial helium abundance to be ∼ 0.9 per cent, ∼ 2 per cent, ∼ 5 per cent, and ∼ 29 per cent, respectively. We report that the systematic uncertainties on the inferred masses and radii arising from the treatment of initial helium abundance in stellar grids lie within the expected accuracy limits of ESA’s PLATO, although this is not the case for the age.

Nsamba, Benard↗

Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains

This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.

54 ENVIRONMENTAL SCIENCES↗

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Proton Exchange Membrane Fuel Cell Engineering Model Powerplant. Test Report: Benchmark Tests in Three Spatial Orientations

Proton exchange membrane (PEM) fuel cell technology is the leading candidate to replace the aging alkaline fuel cell technology, currently used on the Shuttle, for future space missions. This test effort marks the final phase of a 5-yr development program that began under the Second Generation Reusable Launch Vehicle (RLV) Program, transitioned into the Next Generation Launch Technologies (NGLT) Program, and continued under Constellation Systems in the Exploration Technology Development Program. Initially, the engineering model (EM) powerplant was evaluated with respect to its performance as compared to acceptance tests carried out at the manufacturer. This was to determine the sensitivity of the powerplant performance to changes in test environment. In addition, a series of tests were performed with the powerplant in the original standard orientation. This report details the continuing EM benchmark test results in three spatial orientations as well as extended duration testing in the mission profile test. The results from these tests verify the applicability of PEM fuel cells for future NASA missions. The specifics of these different tests are described in the following sections.

Loyselle, Patricia↗

Milestone 1.2.10: Steady-state H 2 “roll over” point data for aluminum alloys 1100 and 6061

Extended (> 50 years) dry storage is being evaluated by the U.S. Department of Energy (DOE) for the disposition of ~ 18 metric tons of aluminum-clad spent nuclear fuel (ASNF). Transition of the current ASNF inventory into dry storage—using the standard DOE canister—necessitates a rigorous, predictive understanding of the long-term physical and chemical factors that may influence the integrity of the proposed storage canister, including radiolytic molecular hydrogen (H 2 ) generation. Current model predictions employ initial radiolytic yields of H 2 , the values of which change as the cladding’s H 2 -precursor inventory is depleted and H 2 itself becomes progressively more involved in radiolytic and surface dissociation processes. Consequently, the absorbed radiation dose that this steady-state H 2 yield corresponds to is essential for the evaluation and improvement of model predictions. Here, we report our findings on the long-term generation of H 2 from the gamma irradiation (≤ 36 MGy) of corroded AA1100 and AA6061 coupons in helium environments at ambient temperature and ~ 50% RH. Our findings show that AA1100 systems reached steady-state by ~ 36 MGy, while higher doses were necessary for AA6061 systems. This discrepancy was attributed to the AA6061 coupons developing a thicker corrosion layer that led to the trapping of H 2 and its precursors, and potentially additional chemistries, ultimately delaying the depletion of H 2 precursors and the system’s “roll over” point. Further, current model predictions—based on previous AA1100 data—do not show steady-state attainment until above 120 MGy, which is not the case for the AA1100 data collected here. Consequently, the new alloy dependent data presented here are important for the continued improvement of predictive computer models for evaluating the feasibility of extended storage of ASNF in helium backfilled canisters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A numerical investigation of severe thunderstorm gust fronts

A numerical model was developed to simulate the evolution and structure of severe thunderstorm gust fronts. The model is a non-hydrostatic, fine resolution, cross-sectional primitive equation model. Two-dimensional horizontal and vertical equations of motion, the continuity equation, and the thermodynamic energy equation were utilized. It was shown that two dominant factors influencing gust front configuration are surface friction and the solenoidal field coincident with the front. It is suggested that solenoidal accelerations oppose the deceleration of surface friction. After a downdraft is initiated in the model, these opposing tendencies soon reach a balance and the gust front achieves a quasi-steady configuration. Thus, the experiments indicate that surface friction does not induce a cycle of front formation and collapse. In addition, the effect of evaporative cooling in producing a vigorous downdraft was parameterized by a local cooling function. Greater cooling in the downdraft results in a more intense gust front that exhibits stronger wind maximums and greater shears. The ambient air stability was shown to be an important factor influencing the depth of the cold outflow.

Mitchell, K. E.↗

Massive superclusters as a probe of the nature and amplitude of primordial density fluctuations

It is pointed out that correlation studies of galaxy positions have been widely used in the search for information about the large-scale matter distribution. The study of rare condensations on large scales provides an approach to extend the existing knowledge of large-scale structure into the weakly clustered regime. Shane (1975) provides a description of several apparent massive condensations within the Shane-Wirtanen catalog, taking into account the Serpens-Virgo cloud and the Corona cloud. In the present study, a description is given of a model for estimating the frequency of condensations which evolve from initially Gaussian fluctuations. This model is applied to the Corona cloud to estimate its 'rareness' and thereby estimate the rms density contrast on this mass scale. An attempt is made to find a conflict between the density fluctuations derived from the Corona cloud and independent constraints. A comparison is conducted of the estimate and the density fluctuations predicted to arise in a universe dominated by cold dark matter.

Kaiser, N.↗

The GOES-R Lightning Mapper Sensor

The Lightning Mapper Sensor on GOES-R builds on previous measurements of lightning from low earth orbit by the OTD (Optical Transient Detector) and LIS (Lightning Imaging Sensor) sensors. Unlike observations from low earth orbit, the GOES-R platform will allow continuous monitoring of lightning activity over the Continental United States and southern Canada, Central and South America, and portions of the Atlantic and Pacific Oceans. The LMS will detect total (cloud-to-ground and intracloud) lightning at storm scale resolution (approx. 8 km) using a highly sensitive Charge Coupled Device (CCD) detector array. Discrimination between lightning optical transients and a bright sunlit background scene is accomplished by employing spectral, spatial, and temporal filtering along with a background subtraction technique. The result is 24 hour detection capability of total lightning. These total lightning observations can be made available to users within about 20 seconds. Research indicates a number of ways that total lightning observations from LMS could benefit operational activities, including 1) potential increases in lead times and reduced false alarms for severe thunderstorm and tornado Warnings, 2) improved routing of &rail around thunderstorms, 3) support for spacecraft launches and landings, 4) improved ability to monitor tropical cyclone intensity, 5) ability to monitor thunderstorm intensification/weakening during radar outages or where radar coverage is poor, 6) better identification of deep convection for the initialization of numerical prediction models, 7) improved forest fire forecasts, 8) identification of convective initiation, 9) identification of heavy convective snowfall, and 10) enhanced temporal resolution of storm evolution (1 minute) than is available from radar observations. Total lightning data has been used in an operational environment since July 2003 at the Huntsville, Alabama National Weather Service office. Total lightning measurements are obtained by the North Alabama Lightning Mapping Array (LMA) and have successfully been used in warning decisions. Every 2 minutes, total lightning counts in 2 km by 2 km horizontal, 1 km vertical grids are available to forecasters on an AWIPS (Advanced Weather Interactive Processing System) workstation. Experience with the LMA total lightning data is used to illustrate the potential use of LMS data that would be available to forecasters across the US. This abstract is for submission as a presentation to the National Weather Association Annual Meeting to be held 16-21 October 2004 in Portland, OR. This abstract will be published in the conference proceedings.

Buechler, Dennis↗

Improving the Representation of Land Surface Processes using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here we apply an Ensemble Kalman Filter DA system as part of the Data Assimilation Research Testbed to a variety of land surface simulations. First, we describe the use of remotely sensed biomass observations to provide improved simulations of plant phenology, carbon and water cycling for regions highly sensitive to climate change (Western US, China, and Arctic). We discuss approaches to account for systemic biases between models and observations, including the use of spatially-varying adaptive ensemble inflation as an alternative approach to re-scaling soil moisture observations. Finally, we discuss a strategy to incorporate complementary observations (snow water equivalent, solar-induced fluorescence) to better constrain the representation of carbon and water cycling across complex terrain.

DART↗

Improving the Representation of Land Surface Processes Using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here, we emphasize the implementation of tools and approaches to overcome challenges related to land DA to constrain carbon and water cycling. In particular, we discuss the implementation of adaptive inflation to modify ensemble spread in response to time-varying networks of gridded observations. We also discuss methods to generate ensemble spread through boundary condition (meteorology) forcing that can be applied to site-level applications. Next, we describe the application of vertical localization upon surface soil moisture observations, and forward operators specifically designed for the assimilation of snow and solar-induced fluorescence observations. Finally, we discuss the potential benefit of a quantile conserving filter used to update bounded quantities (state or parameter values).

Brett Raczka↗

THERMAL STRESS ANALYSIS OF A SPENT NUCLEAR FUEL CANISTER

The potential for chloride induced stress corrosion cracking (CISCC) in spent nuclear fuel dry storage canisters is a current topic of research by the US Department of Energy Spent Fuel and Waste Science and Technology program. One of the important prerequisites for CISCC is a tensile residual stress state. This study utilizes computational models to provide an initial analysis of thermal stresses in a generic spent nuclear fuel canister. A STAR-CCM+ thermal fluid model provides a temperature profile analysis of the canister based on four different ambient temperature conditions. An ANSYS APDL finite element model incorporates the temperature profiles to analyze the thermal stresses in the canister. The calculated thermal stress magnitudes are in the range of 10 MPa to 80 MPa, which could be significant for crack propagation through the canister wall. The maximum thermal stress is not a concern for the structural failure of the canister, but it is high enough that it is a significant contribution to the total stress state of the canister wall, which also includes residual stress from fabrication, internal pressure from helium cover gas, and potential transient mechanical loads, such as earthquakes. This paper describes the initial finite element thermal stress analysis that was completed in 2020, reports the initial findings, and identifies areas where model refinement is still needed to complete this work in the future.

Jensen, Ben J.↗

Chromium Plume at LANL - From Conceptual Site Model to Development of Complex Remediation Strategies - 18682

A hexavalent chromium plume is present in deep groundwater beneath Los Alamos National Laboratory (LANL or the Laboratory). The plume originated from releases of potassium dichromate from cooling towers from 1956 - 1972. Estimates are that up to 72,000 kg of chromium was released during that period. Data from surface water-balance studies, core-hole and groundwater characterization, geomorphic investigations and geology and geophysics provide the basis for an evolved conceptual model that chromium initially traveled rapidly via an effluent-supported stream for several kilometers, then infiltrated through a stratigraphically complex 275-300 meters thick vadose zone and eventually into a deep (regional) aquifer that is the water supply for the Laboratory and surrounding municipality. The travel time from the release point to the regional aquifer is estimated at 5-10 years based on modeling and contaminant tracers. Monitoring data from two-screen wells and piezometers within the regional aquifer, and data derived from stratified sampling conducted during sonic coring and from analysis of sonic core indicate that the chromium contamination occurs within the upper 23 meters of the highly stratified regional aquifer. Different geologic units at the water table and structural dip of depositional bedding appear to have little effect on chromium migration and plume shape and thickness which instead appear to be a manifestation of multiple breakthrough locations, interconnectedness of preferential hydraulic strata, and a negligible vertical gradient. Recent monitoring data indicate that the down-gradient plume edge (as defined by the New Mexico groundwater standard of 50 μg/L) may be advancing near the Laboratory's boundary with a neighboring Indian reservation. The Laboratory is addressing the increasing concentrations with a regulatory-driven interim measure that uses hydraulic control to manage plume migration. The interim measure involves limited-duration (2-3 years) pumping, treatment, and injection to rapidly establish and maintain the 50 μg/L plume edge within the Laboratory boundary while a final remediation strategy is developed and approved by the State regulator. A final remediation strategy will depend heavily on the ability to cost-effectively access the deep aquifer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Spatially Explicit Individual-Based Modeling Approach to Evaluate the Cumulative Effects of Wind Energy Development on the Greater Sage-Grouse: Sensitivity Analysis and Validation of Key Parameters

To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.

17 WIND ENERGY↗

Evaluating skill in predicting the Interdecadal Pacific Oscillation in initialized decadal climate prediction hindcasts in E3SMv1 and CESM1 using two different initialization methods and a small set of start years

Abstract It is a daunting challenge to conduct initialized hindcasts with enough ensemble members and associated start years to form a drifted climatology from which to compute the anomalies necessary to quantify the skill of the hindcasts when compared to observations. This limits the ability to experiment with case studies and other applications where only a few initial years are needed. Here we run a set of hindcasts with CESM1 and E3SMv1 using two different initialization methods for a limited set of start years and use the respective uninitialized free-running historical simulations to form the model climatologies. Since the drifts from the observed initial states in the hindcasts toward the uninitialized model state are large and rapid, after a few years the drifted initialized models approach the uninitialized model climatological errors. Therefore, hindcasts from the limited start years can use the uninitialized climatology to represent the drifted model states after about lead year 3, providing a means to compute forecast anomalies in the absence of a large hindcast sample. There is comparable skill for predicting spatial patterns of multi-year Pacific sea surface temperature anomalies in the domain of the Interdecadal Pacific Oscillation using this method compared to the conventional methodology with a large hindcast data set, though there is a model dependence to the drifts in the two initialization methods.

54 ENVIRONMENTAL SCIENCES↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Automated Model Updating Effect on the Linear and Nonlinear Dynamic Responses of Historical Masonry Structures

Finite element models supply the simulation of the problems to obtain the structural behaviors under different loading conditions. The updating processes based on changing the uncertain parameters such as material properties, boundary conditions, meshing size etc. are made to create actual finite element models in order to represent the real conditions of the problems. This procedure can be made using two methods, namely manual model updating and global/local automatic model updating. The global/local automatic model updating procedures are very popular owing to the applicability to all kind of engineering structures, minimizing the differences nearly to zero, availability of damage localization and structural health monitoring. In this paper a historical masonry armory building selected as an application and its finite element analyses, experimental measurements and automated model updating processes are carried out. The finite element model of the structures was developed with ANSYS software and first three natural frequencies are attained. The experimental dynamic characteristics such as natural frequencies, damping ratios and mode shapes are obtained using ambient vibration method. The differences between experimental and numerical frequencies are minimized up to the range of 0–1%, by using local updating procedure. The linear and nonlinear seismic analyses of the building are conducted on initial and updated models. At the end of the analyses, maximum displacements, principal stresses and strains are given with detail. In addition, crack distributions diagrams are attained to validate the stresses accumulation points. It can be seen that there is an increasing trend in the displacement, stress and strain values after finite element model updating both linear and nonlinear analysis. In addition, linear and nonlinear analysis results are compared with each other to evaluate the nonlinear analysis effect. It is concluded that the displacement and strain values are increased after nonlinear analysis. But, there is not any agreement in stress values.

42 ENGINEERING↗

Short-range precipitation forecasts using assimilation of simulated satellite water vapor profiles and column cloud liquid water amounts

These observing system simulation experiments investigate the assimilation of satellite-observed water vapor and cloud liquid water data in the initialization of a limited-area primitive equations model with the goal of improving short-range precipitation forecasts. The assimilation procedure presented includes two aspects: specification of an initial cloud liquid water vertical distribution and diabatic initialization. The satellite data is simulated for the next generation of polar-orbiting satellite instruments, the Advanced Microwave Sounding Unit (AMSU) and the High-Resolution Infrared Sounder (HIRS), which are scheduled to be launched on the NOAA-K satellite in the mid-1990s. Based on cloud-top height and total column cloud liquid water amounts simulated for satellite data a diagnostic method is used to specify an initial cloud water vertical distribution and to modify the initial moisture distribution in cloudy areas. Using a diabatic initialization procedure, the associated latent heating profiles are directly assimilated into the numerical model. The initial heating is estimated by time averaging the latent heat release from convective and large-scale condensation during the early forecast stage after insertion of satellite-observed temperature, water vapor, and cloud water formation. The assimilation of satellite-observed moisture and cloud water, together withy three-mode diabatic initialization, significantly alleviates the model precipitation spinup problem, especially in the first 3 h of the forecast. Experimental forecasts indicate that the impact of satellite-observed temperature and water vapor profiles and cloud water alone in the initialization procedure shortens the spinup time for precipitation rates by 1-2 h and for regeneration of the areal coverage by 3 h. The diabatic initialization further reduces the precipitation spinup time (compared to adiabatic initialization) by 1 h.

Wu, Xiaohua↗