Search NASA⌕ Search

SEARCH · Search NASA

Results for “initial conditions”

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 127 records · Page 7

Computer program for calculating laminar and turbulent boundary layer development in compressible flow

A computer program is described which performs a numerical integration of the equations of motion for a compressible two-dimensional boundary layer. Boundary layer calculations may be carried out for both laminar and turbulent flow for arbitrary Reynolds number and free stream Mach number distribution on planar or axisymmetric bodies with wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and a rough or a smooth wall. A variety of options are available as initial conditions. The program can generate laminar initial conditions such as Falkner-Skan similarity solutions (so that initial wedge flows can be simulated including Blasius or stagnation point flow) or approximate equilibrium turbulent profiles. Alternatively, initial profile input data can be utilized.

Herring, H. J.↗

Numerical prediction experiments related to the summer 1980 U.S. heat wave

The GLA (Goddard Laboratory for Atmosphere) fourth-order model predictions of the summer 1980 heat wave are examined to assess the skill of the model in numerically forecasting the initiation of the heat wave and the sensitivity of the model predictions to initial conditions and boundary forcing. Results of experiments indicate that the model predictions are sensitive to the initial conditions. The sensitivity of the model to anomalous boundary forcing is found to be generally larger than its sensitivity to initial conditions at extended ranges. The results of the study suggest that models may be useful in forecasting similar events if adequate determinations of the evolving boundary forcing can be made in real time.

Atlas, R.↗

Photoemission of Single Dust Grains for Heliospheric Conditions

Initial results of an experiment to measure the photoemission of single dust grains as a function of far ultraviolet wavelengths are presented. Coulombic forces dominate the interaction of the dust grains in the heliosphere. Knowledge of the charge state of dust grains, whether in a dusty plasma (Debye length < intergrain distance) or in the diffuse interplanetary region, is key to understanding their interaction with the solar wind and other solar system constituents. The charge state of heliospheric grains is primarily determined by primary electron and ion collisions, secondary electron emission and photoemission due to ultraviolet sunlight. We have established a unique experimental technique to measure the photoemission of individual micron-sized dust grains in vacuum. This technique resolves difficulties associated with statistical measurements of dust grain ensembles and non-static dust beams. The photoemission yield of Aluminum Oxide 3-micron grains For wavelengths from 120-300 nm with a spectral resolution of 1 nm FWHM is reported. Results are compared to interplanetary conditions.

Spann, James F., Jr.↗

Response of Winter Forecasts Made with the GLA Fourth Order GCM to Changes in the Horizontal Grid Resolution

The GLA Fourth Order General Circulation Model was modified in 1984 to allow it to be run at a finer horizontal grid resolution. Previously the standard model versions integrated the primitive equations on a 4 degree latitude by 5 degree longitude grid. Forecasts made at this resolution are compared to others which were made using the finer 2 degree latitude by 2 1/2 degree longitude grid resolution. Three winter cases are included with initial conditions from 15 December 1978, 5 January 1979 and 21 January 1979. The 21 January 1979 initial conditions were taken from a GLAS analysis (experiment 2728) which incorporated the GLAS temperature retrievals. The other two cases used European Center Analysis for initial conditions. The forecasts were all verified against European Center Analysis.

Pfaendtner, J.↗

Stability and bifurcation of quasiparallel Alfven solitons

The inverse scattering transformation (IST) is used to study the one-parameter and two-parameter soliton families of the derivative nonlinear Schroedinger (DNLS) equation. The two-parameter soliton family is determined by the discrete complex eigenvalue spectrum of the Kaup-Newell scattering problem and the one-parameter soliton family corresponds to the discrete real eigenvalue spectrum. The structure of the IST is exploited to discuss the existence of discrete real eigenvalues and to prove their structural stability to perturbations of the initial conditions. Also, though the two-parameter soliton is structurally stable in general, it is shown that a perturbation of the initial conditions may change the two-parameter soliton into a degenerate soliton which, in turn, is structurally unstable. This degenerate, or double pole, soliton may bifurcate due to a perturbation of the initial conditions into a pair of one-parameter solitons. If the initial profile is on compact support, then this pair of one-parameter solitons must be compressive and rarefactive respectively. Finally, the Gelfand-Levitan equations appropriate for the double pole soliton are solved.

Hamilton, R. L.↗

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗

The Impacts of Microphysics and Planetary Boundary Layer Physics on Model Simulations of U.S. Deep South Summer Convection

Convection-allowing numerical weather simula- tions have often been shown to produce convective storms that have significant sensitivity to choices of model physical parameterizations. Among the most important of these sensitivities are those related to cloud microphysics, but planetary boundary layer parameterizations also have a significant impact on the evolution of the convection. Aspects of the simulated convection that display sensitivity to these physics schemes include updraft size and intensity, simulated radar reflectivity, timing and placement of storm initi- ation and decay, total storm rainfall, and other storm features derived from storm structure and hydrometeor fields, such as predicted lightning flash rates. In addition to the basic parameters listed above, the simulated storms may also exhibit sensitivity to im- posed initial conditions, such as the fields of soil temper- ature and moisture, vegetation cover and health, and sea and lake water surface temperatures. Some of these sensitivities may rival those of the basic physics sensi- tivities mentioned earlier. These sensitivities have the potential to disrupt the accuracy of short-term forecast simulations of convective storms, and thereby pose sig- nificant difficulties for weather forecasters. To make a systematic study of the quantitative impacts of each of these sensitivities, a matrix of simulations has been performed using all combinations of eight separate microphysics schemes, three boundary layer schemes, and two sets of initial conditions. The first version of initial conditions consists of the default data from large-scale operational model fields, while the second features specialized higher- resolution soil conditions, vegetation conditions and water surface temperatures derived from datasets created at NASA's Short-term Prediction and Operational Research Tran- sition (SPoRT) Center at the National Space Science and Technology Center (NSSTC) in Huntsville, AL. Simulations as outlined above, each 48 in number, were conducted for five midsummer weakly sheared coastal convective events each at two sites, Mobile, AL (MOB) and Houston, TX (HGX). Of special interest to operational forecasters at MOB and HGX were accuracy of timing and placement of convective storm initiation, reflectivity magnitudes and coverage, rainfall and inferred lightning threat.

McCaul, Eugene W., Jr.↗

Mesoscale Simulation Data for Initializing Fast-Time Wake Transport and Decay Models

The fast-time wake transport and decay models require vertical profiles of crosswinds, potential temperature and the eddy dissipation rate as initial conditions. These inputs are normally obtained from various field sensors. In case of data-denied scenarios or operational use, these initial conditions can be provided by mesoscale model simulations. In this study, the vertical profiles of potential temperature from a mesoscale model were used as initial conditions for the fast-time wake models. The mesoscale model simulations were compared against available observations and the wake model predictions were compared with the Lidar measurements from three wake vortex field experiments.

Ahmad, Nashat N.↗

GCM sensitivity to 1982 - 1983 equatorial Pacific sea surface temperature anomalies

The response of the GLAS climate model to the much larger 1982-83 sea surface temperatures (SST) anomalies is discussed. Two separate 75 day experiments (control and anomaly simulation pairs) were started from observed initial conditions on 16 Dec. 1982 and initial conditions on 16 Dec. 1979 taken from a 2 year model control run after one year of simulation, respectively. The January control and anomaly SST fields used in both experiments are given. Notable is the greatly extended region of very warm (approximately equal to 29 C) SST water in the anomaly simulation. The January SST anomaly field is representative of the other months of the experiments, all of which had a much larger region of very warm SST in the anomaly simulation than in the previously noted general circulation model studies. The model used is an improved version of the GLAS B-grid GCM used by Shukla and Wallace (1983). The most important physical change in the model is the inclusion of the surface flux parameterization of Deardorff (1972) as modified by Randall (1976). An important improvement in the model simulations is the removal of the climate drift towards unrealistically high temperatures in the tropics, which was moted by Shukla and Wallace (1983).

Fennessy, M. J.↗

Full-Envelope Launch Abort System Performance Analysis Methodology

The implementation of a new dispersion methodology is described, which dis-perses abort initiation altitude or time along with all other Launch Abort System (LAS) parameters during Monte Carlo simulations. In contrast, the standard methodology assumes that an abort initiation condition is held constant (e.g., aborts initiated at altitude for Mach 1, altitude for maximum dynamic pressure, etc.) while dispersing other LAS parameters. The standard method results in large gaps in performance information due to the discrete nature of initiation conditions, while the full-envelope dispersion method provides a significantly more comprehensive assessment of LAS abort performance for the full launch vehicle ascent flight envelope and identifies performance "pinch-points" that may occur at flight conditions outside of those contained in the discrete set. The new method has significantly increased the fidelity of LAS abort simulations and confidence in the results.

Aubuchon, Vanessa V.↗

Robust Landing Using Time-to-Collision Measurement with Actuator Saturation

This paper considers a landing problem for an MAV that uses only a monocular camera for guidance. Although this sensor cannot measure the absolute distance to the target, by using optical flow algorithms, time-to-collision to the target is obtained. Existing work has applied a simple proportional feedback control to simple dynamics and demonstrated its potential. However, due to the singularity in the time-to-collision measurement around the target, this feedback could require an infinite control action. This paper extends the approach into nonlinear dynamics. In particular, we explicitly consider the saturation of the actuator and include the effect of the aerial drag. It is shown that the convergence to the target is guaranteed from a set of initial conditions, and the boundaries of such initial conditions in the state space are numerically obtained. The paper then introduces parametric uncertainties in the vehicle model and in the time-to-collision measurements. Using an argument similar to the nominal case, the robust convergence to the target is proven, but the region of attraction is shown to shrink due to the existence of uncertainties. The numerical simulation validates these theoretical results.

parametric uncertainty↗

Nonlinear structure formation in flat cosmological models

This paper describes the formation of nonlinear structure in flat (zero curvature) Friedmann cosmological models. We consider models with two components: the usual nonrelativistic component that evolves under gravity and eventually forms the large-scale structure of the universe, and a uniform dark matter component that does not clump under gravity, and whose energy density varies with the scale factor a(t) like a(t)(sup -n), where n is a free parameter. Each model is characterized by two parameters: the exponent n and the present density parameter Omega(sub 0) of the nonrelativistic component. The linear perturbation equations are derived and solved for these models, for the three different cases n = 3, n is greater than 3, and n is less than 3. The case n = 3 is relevant to model with massive neutrinos. The presence of the uniform component strongly reduces the growth of the perturbation compared with the Einstein-de Sitter model. We show that the Meszaros effect (suppression of growth at high redshift) holds not only for n = 4, radiation-dominated models, but for all models with n is greater than 3. This essentially rules out any such model. For the case n is less than 3, we numerically integrate the perturbation equations from the big bang to the present, for 620 different models with various values of Omega(sub 0) and n. Using these solutions, we show that the function f(Omega(sub 0), n) = (a/delta(sub +))d(delta)(sub +)/da, which enters in the relationship between the present density contrast delta(sub 0) and peculiar velocity field u(sub 0) is essentially independent of n. We derive approximate solutions for the second-order perturbation equations. These second-order solutions are tested against the exact solutions and the Zel'dovich approximation for spherically symmetric perturbations in the marginally nonlinear regime (the absolute value of delta is less than or approximately 1). The second-order and Zel'dovich solutions have comparable accuracy, significantly higher than the accuracy of the linear solutions. We then investigate the dependence of the delta(sub 0) - u(sub 0) relationship upon the value of n in the nonlinear regime using the second-order solutions for marginally nonlinear, general perturbations, and the exact solutions for strongly nonlinear, spherically symmetric perturbations. In both cases, we find that the delta(sub 0) - u(sub 0) relationship remains independent of n. We speculate that this result extends to strongly nonlinear, general perturbations as well. This eliminates any hope to determine the presence of the uniform component or the value of n using dynamical methods. Finally, we compute the nonlinear evolution of the skewness of the distribution of values of delta, assuming Gaussian initial conditions. We find that the skewness is not only independent of n, but also of Omega(sub 0). Thus the skewness cannot be used to discriminate among various models with Gaussian initial conditions. However, it can be used for testing the Gaussianity of the initial conditions themselves. We conclude that the uniform component leaves no observable signature in the present large-scale structure of the universe. To determine its presence and nature, we must investigate the relationship between the past and present universe, using redshift-dependent tests.

Martel, Hugo↗

The motion of throw away detectors relative to the space shuttle

The motions of throw away detectors (TAD's) are analyzed using the linearized relative motion equations. The TAD's are to be used in the amps program as diagnostic instruments for making various measurements near the shuttle. The TAD's are ejected from the shuttle in arbitrary directions with small relative velocities (0.1 to 1.0 m/s) their subsequent trajectories relative to the shuttle are analyzed. Initial conditions that are likely to result in recontact between the TAD and the shuttle are identified. The sensitivity of the motion to variations in the initial conditions, possibly resulting from inaccuracy in the ejection mechanism, are analyzed as are effects of atmospheric drag. A targeting method, a method of giving the TAD correct initial conditions such that it will pass through a given point relative to the shuttle at a given time, is developed. The results of many specific cases are presented in graphical form.

Mullins, L. D.↗

Sensitivities of Modeled Tropical Cyclones to Surface Friction and the Coriolis Parameter

In this investigation the sensitivities of a 2-D tropical cyclone (TC) model to surface frictional coefficient and the Coriolis parameter are studied and their implication is discussed. The model used is an axisymmetric version of the latest version of the Goddard cloud ensemble model. The model has stretched vertical grids with 33 levels varying from 30 m near the bottom to 1140 m near the top. The vertical domain is about 21 km. The horizontal domain covers a radius of 962 km (770 grids) with a grid size of 1.25 km. The time step is 10 seconds. An open lateral boundary condition is used. The sea surface temperature is specified at 29C. Unless specified otherwise, the Coriolis parameter is set at its value at 15 deg N. The Newtonian cooling is used with a time scale of 12 hours. The reference vertical temperature profile used in the Newtonian cooling is that of Jordan. The Newtonian cooling models not only the effect of radiative processes but also the effect of processes with scale larger than that of TC. Our experiments showed that if the Newtonian cooling is replaced by a radiation package, the simulated TC is much weaker. The initial condition has a temperature uniform in the radial direction and its vertical profile is that of Jordan. The initial winds are a weak Rankin vortex in the tangential winds superimposed on a resting atmosphere. The initial sea level pressure is set at 1015 hPa everywhere. Since there is no surface pressure perturbation, the initial condition is not in gradient balance. This initial condition is enough to lead to cyclogenesis, but the initial stage (say, the first 24 hrs) is not considered to resemble anything observed. The control experiment reaches quasi-equilibration after about 10 days with an eye wall extending from 15 to 25 km radius, reasonable comparing with the observations. The maximum surface wind of more than 70 m/s is located at about 18 km radius. The minimum sea level pressure on day 10 is about 886 hPa. Thus the overall simulation is considered successful and the model is considered adequate for our investigation.

Chao, Winston C.↗

Improved Modeling of Land-Atmosphere Interactions using a Coupled Version of WRF with the Land Information System

The exchange of energy and moisture between the Earth's surface and the atmospheric boundary layer plays a critical role in many hydrometeorological processes. Accurate and high-resolution representations of surface properties such as sea-surface temperature (SST), vegetation, soil temperature and moisture content, and ground fluxes are necessary to better understand the Earth-atmosphere interactions and improve numerical predictions of weather and climate phenomena. The NASA/NWS Short-term Prediction Research and Transition (SPORT) Center is currently investigating the potential benefits of assimilating high-resolution datasets derived from the NASA moderate resolution imaging spectroradiometer (MODIS) instruments using the Weather Research and Forecasting (WRF) model and the Goddard Space Flight Center Land Information System (LIS). The LIS is a software framework that integrates satellite and ground-based observational and modeled data along with multiple land surface models (LSMs) and advanced computing tools to accurately characterize land surface states and fluxes. The LIS can be run uncoupled to provide a high-resolution land surface initial condition, and can also be run in a coupled mode with WRF to integrate surface and soil quantities using any of the LSMs available in LIS. The LIS also includes the ability to optimize the initialization of surface and soil variables by tuning the spin-up time period and atmospheric forcing parameters, which cannot be done in the standard WRF. Among the datasets available from MODIS, a leaf-area index field and composite SST analysis are used to improve the lower boundary and initial conditions to the LIS/WRF coupled model over both land and water. Experiments will be conducted to measure the potential benefits from using the coupled LIS/WRF model over the Florida peninsula during May 2004. This month experienced relatively benign weather conditions, which will allow the experiments to focus on the local and mesoscale impacts of the high-resolution MODIS datasets and optimized soil and surface initial conditions. Follow-on experiments will examine the utility of such an optimized WRF configuration for more complex weather scenarios such as convective initiation. This paper will provide an overview of the experiment design and present preliminary results from selected cases in May 2004.

Case, Jonathan L.↗

Dynamic Modeling of Off-Nominal Operation in Advanced Life Support Systems

System failures, off-nominal operation, or unexpected interruptions in processing capability can cause unanticipated instabilities in Advanced Life Support (ALS) systems, even long after they are repaired. Much current modeling assumes ALS systems are static and linear, but ALS systems are actually dynamic and nonlinear, especially when failures and off nominal operation are considered. Modeling and simulation provide a way to study the stability and time behavior of nonlinear dynamic ALS systems under changed system configurations or operational scenarios. The dynamic behavior of a nonlinear system can be fully explored only by computer simulation over the full range of inputs and initial conditions. Previous simulations of BIO-Plex in SIMULINK, a toolbox of Matlab, were extended to model the off-nominal operation and long-term dynamics of partially closed physical/chemical and bioregenerative life support systems. System nonlinearity has many interesting potential consequences. Different equilibrium points may be reached for different initial conditions. The system stability can depend on the exact system inputs and initial conditions. The system may oscillate or even in rare cases behave chaotically. Temporary internal hardware failures or external perturbations in ALS systems can lead to dynamic instability and total ALS system failure. Appropriate control techniques can restore reliable operation and minimize the effects of dynamic instabilities due to anomalies or perturbations in a life support system.

Jones, Harry↗

Controlled experiments in cosmological gravitational clustering

A systematic study is conducted of gravitational instability in 3D on the basis of power-law initial spectra with and without spectral cutoff, emphasizing nonlinear effects and measures of nonlinearity; effects due to short and long waves in the initial conditions are separated. The existence of second-general pancakes is confirmed, and it is noted that while these are inhomogeneous, they generate a visually strong signal of filamentarity. An explicit comparison of smoothed initial conditions with smoothed envelope models also reconfirms the need to smooth over a scale larger than any nonlinearity, in order to extrapolate directly by linear theory from Gaussian initial conditions.

Melott, Adrian L.↗

Effects of initial spanwise perturbations and three-dimensionality in a plane mixing layer

An experimental study concerning the influence of small changaes in initial conditions on the near- and far-field evolution of the th ree-dimensional structure in a plane mixing layer has been conducted. A two-stream mixing layer with a velocity ratio of 0.6 was generated with the initial boundary layers laminar and nominally two-dimensional. The initial conditions were changed slightly by interchanging the high- and low-speed sides of the wind tunnel, while maintaining the same velocity ratio. The results indicate that, even with this minor change in initial conditions, the near-field regions of the two cases differ significantly, with the peak Reynolds stress levels varying by up to 100 percent. In addition, the distribution of the individual streamwise vortical structures differs in intensity and position for the two cases, although the overall qualitative description is comparable.

Plesniak, M. W.↗