Search NASA⌕ Search

SEARCH · Search NASA

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,279 records · Page 71

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Bui, Trong T.↗

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

An approach to verification and validation of a reliable multicasting protocol

This paper describes the process of implementing a complex communications protocol that provides reliable delivery of data in multicast-capable, packet-switching telecommunication networks. The protocol, called the Reliable Multicasting Protocol (RMP), was developed incrementally using a combination of formal and informal techniques in an attempt to ensure the correctness of its implementation. Our development process involved three concurrent activities: (1) the initial construction and incremental enhancement of a formal state model of the protocol machine; (2) the initial coding and incremental enhancement of the implementation; and (3) model-based testing of iterative implementations of the protocol. These activities were carried out by two separate teams: a design team and a V&V team. The design team built the first version of RMP with limited functionality to handle only nominal requirements of data delivery. In a series of iterative steps, the design team added new functionality to the implementation while the V&V team kept the state model in fidelity with the implementation. This was done by generating test cases based on suspected errant or off-nominal behaviors predicted by the current model. If the execution of a test was different between the model and implementation, then the differences helped identify inconsistencies between the model and implementation. The dialogue between both teams drove the co-evolution of the model and implementation. Testing served as the vehicle for keeping the model and implementation in fidelity with each other. This paper describes (1) our experiences in developing our process model; and (2) three example problems found during the development of RMP.

Callahan, John R.↗

An Approach to Verification and Validation of a Reliable Multicasting Protocol

This paper describes the process of implementing a complex communications protocol that provides reliable delivery of data in multicast-capable, packet-switching telecommunication networks. The protocol, called the Reliable Multicasting Protocol (RMP), was developed incrementally using a combination of formal and informal techniques in an attempt to ensure the correctness of its implementation. Our development process involved three concurrent activities: (1) the initial construction and incremental enhancement of a formal state model of the protocol machine; (2) the initial coding and incremental enhancement of the implementation; and (3) model-based testing of iterative implementations of the protocol. These activities were carried out by two separate teams: a design team and a V&V team. The design team built the first version of RMP with limited functionality to handle only nominal requirements of data delivery. In a series of iterative steps, the design team added new functionality to the implementation while the V&V team kept the state model in fidelity with the implementation. This was done by generating test cases based on suspected errant or offnominal behaviors predicted by the current model. If the execution of a test was different between the model and implementation, then the differences helped identify inconsistencies between the model and implementation. The dialogue between both teams drove the co-evolution of the model and implementation. Testing served as the vehicle for keeping the model and implementation in fidelity with each other. This paper describes (1) our experiences in developing our process model; and (2) three example problems found during the development of RMP.

Callahan, John R.↗

An Assessment of the Potential Predictability of Drought Over the United States Based on Climate Model Simulations with Specified SST

The USCLIV AR working group on drought recently initiated a series of global climate model simulations forced with idealized SST anomaly patterns, designed to address a number of uncertainties regarding the impact of SST forcing and the role of land-atmosphere feedbacks on regional drought. The runs were done with several global atmospheric models including NASA/NSIPP-l, NCEP/GFS, GFDLlAM2, and NCAR CCM3 and CAM3.5. Specific questions that the runs are designed to address include: What are mechanisms that maintain drought across the seasonal cycle and from one year to the next. To what extent can droughts develop independently of ocean variability due to year-to-year memory that may be inherent to the land. What is the role of the different ocean basins? Here we focus on the potential predictability of drought conditions over the United States. Specific issues addressed include the seasonality and regionality of the signal-to-noise ratios associated with Pacific and Atlantic SST forcing, and the sensitivity of the results to the climatological stationary waves simulated by the different AGCMs.

Schubert, Siegfried↗

Analysis of stratospheric ozone, temperature, and minor constituent data

The objective of this research is to use available satellite measurements of temperature and constituent concentrations to test the conceptual picture of stratospheric chemistry and transport. This was originally broken down into two sub-goals: first, to use the constituent data to search for critical tests of our understanding of stratospheric chemistry and second, to examine constituent transport processes emphasizing interactions with chemistry on various time scales. A third important goal which has evolved is to use the available solar backscattered ultraviolet (SBUV) and Total Ozone Mapping Spectrometer (TOMS) data from Nimbus 7 to describe the morphology of recent changes in Antarctic and global ozone with emphasis on searching for constraints to theories. The major effort now being pursued relative to the two original goals is our effort as a theoretical team for the Arctic Airborne Stratospheric Expedition (AASE). Our effort for the AASE is based on the 3D transport and chemistry model at Goddard. Our goal is to use this model to place the results from the mission data in a regional and global context. Specifically, we set out to make model runs starting in late December and running through March of 1989, both with and without heterogeneous chemistry. The transport is to be carried out using dynamical fields from a 4D data assimilation model being developed under separate funding from this task. We have successfully carried out a series of single constituent transport experiments. One of the things demonstrated by these runs was the difficulty in obtaining observed low N2O abundances in the vortex without simultaneously obtaining very high ozone values. Because the runs start in late December, this difficulty arises in the attempt to define consistent initial conditions for the 3D model. To accomplish a consistent set of initial conditions, we are using the 2D photochemistry-transport model of Jackman and Douglass and mapping in potential temperature, potential vorticity space as developed by Schoeberl and coworkers.

Stolarski, Richard S.↗

Documentation of the Goddard Laboratory for atmospheres fourth-order two-layer shallow water model

The theory and numerical treatment used in the 2-level GLA fourth-order shallow water model are described. This model was designed to emulate the horizontal finite differences used by the GLA Fourth-Order General Circulation Model (Kalnay et al., 1983) in addition to its grid structure, form of high-latitude and global filtering, and time-integration schemes. A user's guide is also provided instructing the user on how to create initial conditions, execute the model, and post-process the data history.

Takacs, L. L.↗

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Kuang, Wei-Jia↗

Heating and melting of small icy satellites by the decay of 26Al

We study the effect of radiogenic heating due to 26Al on the thermal evolution of small icy satellites. Our object is to find the extent of internal melting as a function of the satellite radius and of the initial 26Al abundance. The implicit assumption, based on observations of young stars, is that planet and satellite accretion occurred on a time scale of approximately 10(6) yr (comparable with the lifetime of 26Al). The icy satellites are modeled as spheres of initially amorphous ice, with chondritic abundances of 40K, 232Th, 235U, 238U, corresponding to an ice/dust mass ratio of 1. Evolutionary calculations are carried out, spanning 4.5 x 10(9) yr, for different combinations of the two free parameters. Heat transfer by subsolidus convection is neglected for these small satellites. Our main conclusion is that the initial 26Al abundance capable of melting icy bodies of satellite size to a significant extent is more than 10 times lower than that prevailing in the interstellar medium (or that inferred from the Ca-Al rich inclusions of the Allende meteorite, approximately 7 x 10(-7) by mass). We find, for example, that an initial 26Al mass fraction of approximately 4 x 10(-8) is sufficient for melting almost completely icy spheres with radii of 800 km, typical of the larger icy planetary satellites. We also find that for any given 26Al abundance, there is a narrow range of radii below which only marginal melting occurs and above which most of the ice melts (and refreezes later). Since extensive melting may have important consequences, such as differentiation, gas release, and volcanic activity, the effect of 26Al should be included in future studies of satellite interiors.

Non-NASA Center↗

Determining the Shallow Surface Velocity at the Apollo 17 Landing Site

Many studies have been performed to determine the shallow surface velocity model at the Apollo 17 landing site. The Lunar Seismic Profiling Experiment (LSPE) had both an active component with eight explosive packages (EPs) and a passive experiment collecting data at various time intervals. Using the eight EPs, the initial shallow surface velocity model was determined to be 250 m/s in the first layer of depth 248 m, 1200 m/s with a depth of 927 m in the second layer, and 4000 m/s down to a depth of 2 km in the third layer. Have performed variations on this study to produce new velocity models shown. Recent studies have also been reanalyzing the passive LSPE data and have found three different thermal moonquake event types occurring at different times within the lunar day. The current goal of the project is to collocate the thermal moonquakes to physical surface features to determine the breakdown of lunar rocks. However, to locate shallow surface events, an accurate velocity model is needed. Presented a thermal moonquake location algorithm using first order approximation, including surface events only. To improve these approximations, a shallow surface velocity is needed.

Phillips, D.↗

Impact of the initial specification of moisture and vertical motion on precipitation forecasts with a mesoscale model Implications for a satellite mesoscale data base

Using a numerical mesoscale model, four simulations were performed to determine the effects of suppressing the initial mesoscale information in the moisture and wind fields on the precipitation forecasts. The simulations included a control forecast 12-h simulation that began at 1200 GMT March 1982 and three experiment simulations with modifications to the moisture and vertical motion fields incorporated at 1800 GMT. The forecasts from 1800 GMT were compared to the second half of the control forecast. It was found that, compared to the control forecast, suppression of the moisture and/or wind initial field(s) produces a drier forecast. However, the characteristics of the precipitation forecasts of the experiments were not different enough to conclude that either mesoscale moisture or mesoscale vertical velocity at the initial time are more important for producing a forecast closer to that of the control.

Mlynczak, Pamela E.↗

Ultraviolet imaging telescope and optical emission-line observations of H II regions in M81

Images of the type Sab spiral galaxy M81 were obtained in far-UV and near-UV bands by the Ultraviolet Imaging Telescope (UIT) during the Astro-1 Spacelab mission of 1990 December. Magnitudes in the two UV bands are determined for 52 H II regions from the catalog of Petit, Sivan, & Karachentsev (1988). Fluxes of the H-alpha and H-beta emission lines are determined from CCD images. Extinctions for the brightest H II regions are determined from observed Balmer decrements. Fainter H II regions are assigned the average of published radio-H-alpha extinctions for several bright H II regions. The radiative transfer models of Witt, Thronson, & Capuano (1992) are shown to predict a relationship between Balmer Decrement and H-alpha extinction consistent with observed line and radio fluxes for the brightest 7 H II regions and are used to estimate the UV extinction. Ratios of Lyman continuum with ratios predicted by model spectra computed for initial mass function (IMF) slope equal to -1.0 and stellar masses ranging from 5 to 120 solar mass. Ages and masses are estimated by comparing the H-alpha and far-UV fluxes and their ratio with the models. The total of the estimated stellar masses for the 52 H II regions is 1.4 x 10(exp 5) solar mass. The star-formation rate inferred for M81 from the observed UV and H-alpha fluxes is low for a spiral galaxy at approximately 0.13 solar mass/yr, but consistent with the low star-formation rates obtained by Kennicutt (1983) and Caldwell et al. (1991) for early-type spirals.

Hill, Jesse K.↗

Influence of radiative heating and cumulus convection on development of mean monsoon circulation in July

Numerical simulations of July mean monsoon circulation in the tropics are described. The model used in the simulations was based on a series of primitive equations for the combined effects of variations of solar radiation, radiative diurnal warming, and large-scale and deep cumulus condensation, and the kinematic effects of topography. The initial states of the model were derived from the observed mean distributions of pressure and humidity. Analysis of the numerical results showed that the large-scale features of the mean July monsoon circulation in the tropics are created mainly by differential diabatic heating under the influence of the specific topography. The time necessary to establish the large scale features was only about 5 days when the diurnal variation of solar radiation was taken into account. Graphic illustrations of the simulated mean July flow conditions are provided.

Kuo, H. L.↗

Antarctic springtime ozone depletion computed from temperature observations

An observationally based, mechanistic dynamical model is used to simulate the decline of total ozone during September and October for the years 1979 through 1986. Vertical velocities derived from observed stratospheric temperature changes and computed radiative heating rates are used to advect an ozone mixing ratio profile during the Antarctic spring period. An early August 1982 Syowa balloonsonde ozone profile is used to initialize the computations. The model reasonably simulates the September and October changes in total ozone, considering the uncertainties in the observed data and the radiative heating. The simulated decline is found to be very sensitive to the choice of initial ozone profile and to small changes in the radiative heating. The results of this study suggest that the dynamical hypothesis of the Antarctic ozone depletion is both quantitatively credible and consistent with the observed temperature changes.

Rosenfield, Joan E.↗

Heating and melting of small icy satellites by the decay of Al-26

The effect of radiogenic heating due to Al-26 on the thermal evolution of small icy satellites is studied. The object is to find the extent of internal melting as a function of the satellite radius and of the initial Al-26 abundance. The implicit assumption, based on observations of young stars, is that planet and satellite accretion occurred on a time scale of about 10 to the 6th yr (comparable with the lifetime of Al-26. The icy satellites are modeled as spheres of initially amorphous ice, with chondritic abundances of K-40, Th-232, U-235, and U-238, corresponding to an ice/dust mass ratio of 1. Evolutionary calculations are carried out, spanning 4.5 x 10 to the 9th yr, for different combinations of the two free parameters. Heat transfer by subsolidus convection is neglected for these small satellites. The main conclusion is that the initial Al-26 abundance capable of melting icy bodies of satellite size to a significant extent is more than 10 times lower than that prevailing in the interstellar medium (or that inferred from the Ca-Al rich inclusions of the Allende meteorite, about 7 x 10 to the -7th by mass).

Prialnik, Dina↗

The utilization of Seasat-A data in numerical weather prediction

The present paper provides the results from several numerical experiments which have been conducted in cases in which SASS (Seasat-A Satellite Scatterometer) winds were either included or excluded. Global scale experiments are discussed along with a regional area model. The regional model is essentially that described by Duffy (1981). A description of the initial conditions used in the regional model experiment is also provided. It is found that Seasat-A data can be important in defining the initial conditions which results in rapid maritime cyclogenesis. Currently conducted work is concerned with the development of a higher horizontal and vertical resolution global model, with an improved planetary boundary layer parameterization.

Duffy, D. G.↗

The homunculus of Eta Carinae: An interacting stellar winds paradigm

We simulate the origin and evolution of the bipolar nebula surrounding Eta Car using numerical two-dimensional gasdynamic models. The generalized interacting stellar winds scenario, wherein a stellar wind interacts with an aspherical circumstellar environment, is adopted. The eruption wind of 1840-1860, which is taken to be spherically symmetric, interacts with a preeruption toroidal density environment. Using reasonable assumptions of initial conditions and eruption parameters based on archival data, we have performed over 30 simulations in an effort to bracket the initial parameters which produce models that best match observations. We find that models with high pole-to-equator density contrasts (greater than 100) and toroidal density configurations nicely account for the observed morphology and kinematics of the homunculus.

Frank, Adam↗

A software development and evolution model based on decision-making

Design is a complex activity whose purpose is to construct an artifact which satisfies a set of constraints and requirements. However the design process is not well understood. The software design and evolution process is the focus of interest, and a three dimensional software development space organized around a decision-making paradigm is presented. An initial instantiation of this model called 3DPM(sub p) which was partly implemented, is presented. Discussion of the use of this model in software reuse and process management is given.

Wild, J. Christian↗