SEARCH · Engineering Papers

Results for “generalized additive models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗

NASA NTRS1971

Vibration characteristics of Z-ring-stiffened 60 deg conical shell models of a planetary entry spacecraft

An experimental investigation of the vibration characteristics of a 60 deg conical shell model of a planetary entry vehicle is described and the results presented. Model configurations include the shell with or without one or two Z-ring stiffeners and with or without a simulated payload. Tests were conducted with the model clamped at the small diameter and with the model suspended at the simulated payload. Additionally, calculated results obtained from application of several analytical procedures reported in the literature are presented together with comparisons between experimental and calculated frequencies and meridional mode shapes. Generally, very good frequency agreement between experimental and calculated results was obtained for all model configurations. For small values of circumferential mode number, however, the frequency agreement decreased as the number of ring stiffeners increased. Overall agreement between experimental and calculated mode shapes was generally good. The calculated modes usually showed much larger curvatures in the vicinity of the rings than were observed in the experimentally measured mode shapes. Dual resonances associated with modal preference were noted for the shell without Z-ring stiffeners, whereas the addition of stiffeners produced resonances for which the model responded in two or more modes over different sections of the shell length.

Naumann, E. C.↗

DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗

NASA NTRS2017

A Time-Regularized Multiple Gravity-Assist Low-Thrust Bounded-Impulse Model for Trajectory Optimization

The multiple gravity assist low-thrust (MGALT) trajectory model combines the medium-fidelity Sims-Flanagan bounded-impulse transcription with a patched-conics flyby model and is an important tool for preliminary trajectory design. While this model features fast state propagation via Kepler's equation and provides a pleasingly accurate estimation of the total mass budget for the eventual flight-suitable integrated trajectory it does suffer from one major drawback, namely its temporal spacing of the control nodes. We introduce a variant of the MGALT transcription that utilizes the generalized anomaly from the universal formulation of Kepler's equation as a decision variable in addition to the trajectory phase propagation time. This results in two improvements over the traditional model. The first is that the maneuver locations are equally spaced in generalized anomaly about the orbit rather than time. The second is that the Kepler propagator now has the generalized anomaly as its independent variable instead of time and thus becomes an iteration-free propagation method. The new algorithm is outlined, including the impact that this has on the computation of Jacobian entries for numerical optimization, and a motivating application problem is presented that illustrates the improvements that this model has over the traditional MGALT transcription.

Optimization↗

NASA NTRS1997

Diurnal Forcing of Planetary Atmospheres

Much progress has been made on calculations of the Martian seasonal water cycle using the Mars Climate Model developed for this purpose. Two papers, documenting the model and the water transport results obtained with it have been published in the Journal of Geophysical Research - Planets. An additional paper describing results related to the evolution of the seasonal water cycle as a result of orbital changes was published in Advances in Space Research. Since that time, further studies have concentrated on the consequences of the soil adsorption required to match the observed water cycle and its relation to the stability of ground ice and other potential water reservoirs. Earth-related studies have concentrated on incorporating an efficient and realistic microphysical model into the Ames Stratospheric General Circulation Model used to simulate the spread of the ML Pinatubo and other volcanic clouds in the stratosphere. In addition, visualizations of the simulations are being incorporated into a video describing the UARS mission. A paper describing the new stratospheric aerosol microphysics package (and its consequences for volcanic cloud evolution) will be submitted in the near future. The paper will discuss the relative importance of condensation and coagulation to early particle growth and the separation of the cloud by sedimentation of the larger particles. A more general paper which highlights the observation that particle number densities did not increase dramatically after the ML Pinatubo eruption is planned. Simulations of atmospheric transport will be extended to include studies of terrestrial tropospheric tracers using the Fifth-Generation Penn State/NCAR Mesoscale Model.

Houben, Howard C.↗

DOE OSTI2022

Search for a W' boson decaying to a vector-like quark and a top or bottom quark in the all-jets final state at $\sqrt{\mathrm{s}}$ = 13 TeV

A search is presented for a heavy W' boson resonance decaying to a B or T vector-like quark and a t or a b quark, respectively. The analysis is performed using proton-proton collisions collected with the CMS detector at the LHC. The data correspond to an integrated luminosity of 138 fb -1 at a center-of-mass energy of 13 TeV. Both decay channels result in a signature with a t quark, a Higgs or Z boson, and a b quark, each produced with a significant Lorentz boost. The all-hadronic decays of the Higgs or Z boson and of the t quark are selected using jet substructure techniques to reduce standard model backgrounds, resulting in a distinct three-jet W' boson decay signature. No significant deviation in data with respect to the standard model background prediction is observed. Upper limits are set at 95% confidence level on the product of the W' boson cross section and the final state branching fraction. A W' boson with a mass below 3.1 TeV is excluded, given the benchmark model assumption of democratic branching fractions. In addition, limits are set based on generalizations of these assumptions. These are the most sensitive limits to date for this final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Modeling PAH Mixture Interactions in a Human In Vitro Organotypic Respiratory Model

One of the most significant challenges in human health risk assessment is to evaluate hazards from exposure to environmental chemical mixtures. Polycyclic aromatic hydrocarbons (PAHs) are a class of ubiquitous contaminants typically found as mixtures in gaseous and particulate phases in ambient air pollution associated with petrochemicals from Superfund sites and the burning of fossil fuels. However, little is understood about how PAHs in mixtures contribute to toxicity in lung cells. To investigate mixture interactions and component additivity from environmentally relevant PAHs, two synthetic mixtures were created from PAHs identified in passive air samplers at a legacy creosote site impacted by wildfires. The primary human bronchial epithelial cells differentiated at the air–liquid interface were treated with PAH mixtures at environmentally relevant proportions and evaluated for the differential expression of transcriptional biomarkers related to xenobiotic metabolism, oxidative stress response, barrier integrity, and DNA damage response. Component additivity was evaluated across all endpoints using two independent action (IA) models with and without the scaling of components by toxic equivalence factors. Both IA models exhibited trends that were unlike the observed mixture response and generally underestimated the toxicity across dose suggesting the potential for non-additive interactions of components. Overall, this study provides an example of the usefulness of mixture toxicity assessment with the currently available methods while demonstrating the need for more complex yet interpretable mixture response evaluation methods for environmental samples.

3D in vitro models↗

DOE OSTI2026

Bayesian Optimization of Catalysis with In-Context Learning

Large language models (LLMs) can perform accurate classification with zero or few examples through in-context learning (ICL), allowing the model to observe query-relevant examples at inference time and eliminating the need for additional weight updates to generalize beyond its original training data. We extend this capability to regression with uncertainty estimation using frozen LLMs (e.g., GPT-4o, Gemini), enabling Bayesian optimization (BO) in natural language without explicit model training or feature engineering. We apply this to materials discovery by representing materials as synthesis and testing procedures for use in natural language prompts. This Bayesian, design-first approach prioritizes optimization toward target material properties before detailed characterization, in contrast to conventional experimental workflows that often emphasize characterization of suboptimal materials. On benchmarks like aqueous solubility and oxidative coupling of methane (OCM), BO-ICL matches or outperforms Gaussian processes. In live experiments on the reverse water–gas shift (RWGS) reaction, BO-ICL identifies multimetallic catalysts that approach equilibrium CO yield within 6 and 10 iterations from a pool of 3,700 and 360,000 candidates, respectively. Our method redefines materials representation and accelerates discovery, with broad applications across catalysis, materials science, and AI.

Calibration↗

DOE OSTI2022

Theory, phenomenology, and experimental avenues for dark showers: a Snowmass 2021 report

In this work, we consider the case of a strongly coupled dark/hidden sector, which extends the Standard Model (SM) by adding an additional non-Abelian gauge group. These extensions generally contain matter fields, much like the SM quarks, and gauge fields similar to the SM gluons. We focus on the exploration of such sectors where the dark particles are produced at the LHC through a portal and undergo rapid hadronization within the dark sector before decaying back, at least in part and potentially with sizeable lifetimes, to SM particles, giving a range of possibly spectacular signatures such as emerging or semi-visible jets. Other, non-QCD-like scenarios leading to soft unclustered energy patterns or glueballs are also discussed. After a review of the theory, existing benchmarks and constraints, this work addresses how to build consistent benchmarks from the underlying physical parameters and present new developments for the PYTHIA Hidden Valley module, along with jet substructure studies. Finally, a series of improved search strategies is presented in order to pave the way for a better exploration of the dark showers at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1994

The extratropical 40-day oscillation in the UCLA general circulation model. Part 1: Atmospheric angular momentum

Variations in atmospheric angular momentum (AAM) are examined in a three-year simulation of the large-scale atmosphere with perpetual January forcing. The simulation is performed with a version of the University of California at Los Angeles (UCLA) general circulation model that contains no tropical Madden-Julian Oscillation (MJO). In addition, the results of three shorter experiments with no topography are analyzed. The three-year standard topography run contains no significant intraseasonal AAM periodicity in the tropics, consistent with the lack of the MJO, but produces a robust, 42-day AAM oscillation in the Northern Hemisphere (NH) extratropics. The model tropics undergoes a barotropic, zonally symmetric oscillation, driven by an exchange of mass with the NH extratropics. No intraseasonal periodicity is found in the average tropical latent heating field, indicating that the model oscillation is dynamically rather than thermodynamically driven. The no-mountain runs fail to produce an intraseasonal AAM oscillation, consistent with a topographic origin for the NH extratropical oscillation in the standard model. The spatial patterns of the oscillation in the 500-mb height field, and the relationship of the extratropical oscillation to intraseasonal variations in the tropics, will be discussed in Part 2 of this study.

Marcus, S. L.↗

NASA NTRS1997

The Effect of Enhanced Diabatic Heating on Stratospheric Circulation

The objective of this research focuses on the stratospheric dynamical response to the increase in aerosol loading and subsequent enhanced diabatic heating resulting from the eruption of Mt. Pinatubo. The Langley research Center three dimensional general circulation model and modifications made to that model for this study are described (addition of hydrogen fluoride tracer and diabatic heating enhancement). Unperturbed hydrogen fluoride distribution is compared to the hydrogen fluoride distribution measured by HALOE. A comparison of control and perturbed model runs is presented.

Kleb, Mary M.↗

NASA NTRS2004

Seasonal-to-Interannual Variability and Land Surface Processes

Atmospheric chaos severely limits the predictability of precipitation on subseasonal to interannual timescales. Hope for accurate long-term precipitation forecasts lies with simulating atmospheric response to components of the Earth system, such as the ocean, that can be predicted beyond a couple of weeks. Indeed, seasonal forecasts centers now rely heavily on forecasts of ocean circulation. Soil moisture, another slow component of the Earth system, is relatively ignored by the operational seasonal forecasting community. It is starting, however, to garner more attention. Soil moisture anomalies can persist for months. Because these anomalies can have a strong impact on evaporation and other surface energy fluxes, and because the atmosphere may respond consistently to anomalies in the surface fluxes, an accurate soil moisture initialization in a forecast system has the potential to provide additional forecast skill. This potential has motivated a number of atmospheric general circulation model (AGCM) studies of soil moisture and its contribution to variability in the climate system. Some of these studies even suggest that in continental midlatitudes during summer, oceanic impacts on precipitation are quite small relative to soil moisture impacts. The model results, though, are strongly model-dependent, with some models showing large impacts and others showing almost none at all. A validation of the model results with observations thus naturally suggests itself, but this is exceedingly difficult. The necessary contemporaneous soil moisture, evaporation, and precipitation measurements at the large scale are virtually non-existent, and even if they did exist, showing statistically that soil moisture affects rainfall would be difficult because the other direction of causality - wherein rainfall affects soil moisture - is unquestionably active and is almost certainly dominant. Nevertheless, joint analyses of observations and AGCM results do reveal some suggestions of land-atmosphere feedback in the observational record, suggestions that soil moisture can affect precipitation over seasonal timescales and across certain large continental areas. The strength of this observed feedback in nature is not large but is still significant enough to be potentially useful, e.g., for forecasts. This talk will address all of these issues. It will begin with a brief overview of land surface modeling in atmospheric models but will then focus on recent research - using both observations and models - into the impact of land surface processes on variability in the climate system.

Koster, Randal↗

DOE OSTI2024

Revealing the unusual rate-dependent mechanical behaviors of nematic liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit unique mechanical properties of soft elasticity and enhanced energy dissipation with rate dependency. They are potentially transformative materials for applications in mechanical impact mitigation and vibration isolation. However, previous studies have primarily focused on the mechanics of LCEs under equilibrium and quasistatic loading conditions. Critical knowledge gaps exist in understanding their rate-dependent behaviors, which are a complex mixture of traditional network viscoelasticity and the soft elastic behaviors with changes in the mesogen orientation and order parameter. Together, these inelastic mechanisms lead to unusual rate-dependent energy absorption responses of LCEs. In this work, we developed a viscoelastic constitutive theory for monodomain nematic LCEs to investigate how multiple underlying sources of inelasticity manifest in the rate-dependent and dissipative behaviors of monodomain LCEs. Here, the theoretical modeling framework combines the neo-classical network theory with evolution rules for the mesogen orientation and order parameter with conventional viscoelasticity. The model is calibrated with uniaxial tension and compression data spanning six decades of strain rates. The established 3D constitutive model enables general loading predictions taking the initial mesogen orientation and order parameter as inputs. Additionally, parametric studies were performed to further understand the rate dependence of monodomain LCEs in relation to their energy absorption characteristics. Based on the parametric studies, particularly loading scenarios are identified as conditions where LCEs outperform conventional elastomers regarding energy absorption.

42 ENGINEERING↗

NASA NTRS2014

Selection of Thermal Worst-Case Orbits via Modified Efficient Global Optimization

Efficient Global Optimization (EGO) was used to select orbits with worst-case hot and cold thermal environments for the Stratospheric Aerosol and Gas Experiment (SAGE) III. The SAGE III system thermal model changed substantially since the previous selection of worst-case orbits (which did not use the EGO method), so the selections were revised to ensure the worst cases are being captured. The EGO method consists of first conducting an initial set of parametric runs, generated with a space-filling Design of Experiments (DoE) method, then fitting a surrogate model to the data and searching for points of maximum Expected Improvement (EI) to conduct additional runs. The general EGO method was modified by using a multi-start optimizer to identify multiple new test points at each iteration. This modification facilitates parallel computing and decreases the burden of user interaction when the optimizer code is not integrated with the model. Thermal worst-case orbits for SAGE III were successfully identified and shown by direct comparison to be more severe than those identified in the previous selection. The EGO method is a useful tool for this application and can result in computational savings if the initial Design of Experiments (DoE) is selected appropriately.

Moeller, Timothy M.↗

DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗

NASA NTRS2022

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

CFD modeling↗

NASA NTRS1985

Some Recent Developments in Numerical Modelling at ECMWF

A new atmospheric model was introduced into operational forecasting at the European Center for Medium Range Weather Forecasts, on 21 April 1983. The principal differences between this model and the Center's first operational model were in the adiabatic formulation, which, in the new model, includes use of a spectral representation in the horizontal, a more general vertical coordinate, and a modified, more-efficient, time-stepping scheme. In addition, new programming techniques and standards were adopted to facilitate both the model's use as a research tool and its adaptation to make full use of the features of the recently-acquired CRAY X-MP computer. A number of revisions were also made in detailed aspects of the formulation of the parameterization schemes. The operational change to this new model was accompanied by a second important change, namely the use of a higher envelope orography in the lower boundary conditions of the model. Problems in the operational performance of the model and development preparatory to the future operational implementation of a high resolution version of the model on the CRAY X-MP are discussed.

Simmons, A. J.↗