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 361 records · Page 20

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a nondimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗

Multiple large-scale coherent mode interactions in a developing round jet

The integral energy method has been used to study the nonlinear interactions of the large-scale coherent structure in a spatially developing round jet. The streamwise development of a jet is obtained in terms of the mean flow shear-layer momentum thickness, the wave-mode kinetic energy and the wave-mode phase angle. With the energy method, a system of partial differential equations is reduced to a system of ordinary differential equations. The nonlinear differential equations are solved with initial conditions which are given at the nozzle exit. It is shown that the initial wave-mode energy densities as well as the initial phase angles play a significant role in the streamwise evolution of the large-scale coherent wave modes and the mean flow.

Lee, Sang S.↗

Three-dimensional Structure of Straight and Curved Plane Wakes

Although the plane wake is marked by the formation of strong spanwise vortices, the initially two-dimensional Karman-like vortices soon develop a three-dimensional structure in the form of secondary streamwise vortices. So far, this streamwise vortex structure has been studied mostly through flow visualization and at relatively low Reynolds numbers. The primary objective of the present program was to investigate the origin and evolution of the three-dimensional structure of straight and curved plane wakes at relatively high Reynolds numbers (Re(sub b) = 28,000) through detailed measurements of the mean and turbulent properties at several streamwise locations. The experiments were conducted in three phases. In the first phase, the development of a straight plane wake was investigated. In the second phase, the effects of imposed streamwise curvature on the wake development were examined. The streamwise curvature was of constant radius and very mild in terms of the curvature ratio (b/(square root of R) is less than 2 percent). In both the first and second phases, the role of initial conditions was examined in wakes generated from both untripped (laminar) and tripped (turbulent) initial boundany layers. In the third phase, the effects of injecting streamwise vorticity and the effects of increased Reynolds number on the tripped wake structure and development were investigated.

Weygandt, James H.↗

Stability of the Dual-Frequency Radar Equations and a New Method Applied to the GPM’s Dual-Frequency Precipitation Radar (DPR) Data

A new algorithm is proposed that estimates two parameters of the particle size distribution (PSD) at each range bin from the Global Precipitation Measurement’s (GPM’s) Dual-frequency Precipitation Radar data. The equation that expresses the relationship of the PSD parameters between adjacent range bins is derived. By including the attenuation effect within the bin in the discretized equation, the new algorithm alleviates the double-solution problem when attenuation within the bin is sufficiently large. The stability of the solutions to the equation depends on the value of the mean diameter D(sub m) and its gradient with respect to range in the case of liquid precipitation. If the critical diameter above which the dual-frequency ratio (DFR) of radar reflectivity factors becomes a monotonically increasing function of D(sub m), the backward processing of the equation provides a stable or moderately diverging solution, unlike the forward processing that often gives unstable solutions. To provide a set of initial conditions without using the surface reference technique (SRT) in the backward processing, an initialization method using the Hitschfeld-Bordan attenuation correction method is proposed and tested. The proposed algorithm may provide a tool for investigating the assumptions used in various algorithms.

GPM↗

CRRES: The combined release and radiation effects satellite program directory

As a result of natural processes, plasma clouds are often injected into the magnetosphere. These chemical releases can be used to study many aspects of such injections. When a dense plasma is injected into the inner magnetosphere, it is expected to take up the motion of the ambient plasma. However, it has been observed in previous releases at moderate altitudes that the cloud preserved its momentum for some time following the release and that parts of the cloud peeled off from the main cloud presumable due to the action of an instability. As one moves outward into the magnetosphere, the mirror force becomes less dominant and the initial conditions following a release are dominated by the formation of a diamagnetic cavity since the initial plasma pressure from the injected Ba ions is greater than the magnetic field energy density. A previous high-altitude release (31,300 km) showed this to be the case initially, but at later times there was evidence for acceleration of the Ba plasma to velocities corresponding to 60,000 K. This effect is not explained. This series of experiments is therefore designed to inject plasma clouds into the magnetosphere under widely varying conditions of magnetic field strength and ambient plasma density. In this way the coupling of injected clouds to the ambient plasma and magnetic field, the formation of striations due to instabilities, and possible heating and acceleration of the injected Ba plasma can be studied over a wide range of magnetosphere parameters. Adding to the scientific yield will be the availability of measurements for the DOD/SPACERAD instruments which can monitor plasma parameters, electric and magnetic fields, and waves before, during and after the releases.

Layman, Laura D.↗

An initial state perturbation experiment with the GISS model

Monthly mean global forecasts for January 1975 have been computed with the Goddard Institute for Space Studies model from four slightly different sets of initial conditions - a 'control' state and three random perturbations thereof - to simulate the effects of initial state uncertainty on forecast quality. Differences among the forecasts are examined in terms of energetics, synoptic patterns and forecast statistics. The 'noise level' of the model predictions is depicted on global maps of standard deviations of sea level pressures, 500 mb heights and 850 mb temperatures for the set of four forecasts. Initial small-scale random errors do not appear to result in any major degradation of the large-scale monthly mean forecast beyond that generated by the model itself, nor do they appear to represent the major source of large-scale forecast error.

Spar, J.↗

Numerical prediction of the large scale tropical flow

The mean error characteristics of two series of forecasts of the tropical flow were studied. By using as initial conditions analyses made with and without the FGGE special observing system, the impact that initial data has on the accuracy of the forecasts is estimated.

Kalnay, E.↗

The influence of land surface properties on Sahel climate. Part 1: Desertification

This is a general circulation model sensitivity study of the physical mechanisms of the effects of desertification on the Sahel drought. The model vegetation types were changed in the prescribed desertification area, which led to changes in the surface characteristics. The model was integrated for three months (June, July, August) with climatological surface conditions (control) and desertification conditions (anomaly) to examine the summer season response to the changed surface conditions. The control and anomaly experiments consisted of five pairs of integrations with different initial conditions and/or sea surface temperature boundary conditions. In the desertification experiment, the moisture flux convergence and rainfall were reduced in the test area and increased to the immediate south of this area. The simulated anomaly dipole pattern was similar to the observed African drought patterns in which the axis of the maximum rainfall shifts to the south. The circulation changes in the desertification experiment were consistent with those observed during sub-Saharan dry years. The tropical easterly jet was weaker and the African easterly jet was stronger than normal. Further, in agreement with the observations, the easterly wave disturbances were reduced in intensity but not in number. Descending motion dominated the desertification area. The surface energy budget and hydrological cycle were also changed substantially in the anomaly experiment.

Xue, Yongkang↗

Spurious fields in time domain computations of scattering problems

In this paper two-dimensional electromagnetic scattering problems with a time-periodic incident field are considered. The scatterer is a perfect conductor, and an artificial boundary condition is used. The large time behavior of solutions, depending on (divergence-free) initial conditions, is characterized. It turns out that in addition to the expected time-periodic solution the limiting solution may also contain a spurious stationary field. The source of the stationary field is explained and equations describing it are obtained. Several avoidance strategies are discussed, and numerical comparisons of these techniques are given.

Kangro, Urve↗

On hierarchical solutions to the BBGKY hierarchy

It is thought that the gravitational clustering of galaxies in the universe may approach a scale-invariant, hierarchical form in the small separation, large-clustering regime. Past attempts to solve the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy in this regime have assumed a certain separable hierarchical form for the higher order correlation functions of galaxies in phase space. It is shown here that such separable solutions to the BBGKY equations must satisfy the condition that the clustered component of the solution has cluster-cluster correlations equal to galaxy-galaxy correlations to all orders. The solutions also admit the presence of an arbitrary unclustered component, which plays no dyamical role in the large-clustering regime. These results are a particular property of the specific separable model assumed for the correlation functions in phase space, not an intrinsic property of spatially hierarchical solutions to the BBGKY hierarchy. The observed distribution of galaxies does not satisfy the required conditions. The disagreement between theory and observation may be traced, at least in part, to initial conditions which, if Gaussian, already have cluster correlations greater than galaxy correlations.

Hamilton, A. J. S.↗

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber and Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Rogers, Michael M.↗

Magellan observations of extended impact crater related features on the surface of Venus

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber & Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Campbell, D. B.↗

Next generation initiation techniques

Four-dimensional data assimilation strategies can generally be classified as either current or next generation, depending upon whether they are used operationally or not. Current-generation data-assimilation techniques are those that are presently used routinely in operational-forecasting or research applications. They can be classified into the following categories: intermittent assimilation, Newtonian relaxation, and physical initialization. It should be noted that these techniques are the subject of continued research, and their improvement will parallel the development of next generation techniques described by the other speakers. Next generation assimilation techniques are those that are under development but are not yet used operationally. Most of these procedures are derived from control theory or variational methods and primarily represent continuous assimilation approaches, in which the data and model dynamics are 'fitted' to each other in an optimal way. Another 'next generation' category is the initialization of convective-scale models. Intermittent assimilation systems use an objective analysis to combine all observations within a time window that is centered on the analysis time. Continuous first-generation assimilation systems are usually based on the Newtonian-relaxation or 'nudging' techniques. Physical initialization procedures generally involve the use of standard or nonstandard data to force some physical process in the model during an assimilation period. Under the topic of next-generation assimilation techniques, variational approaches are currently being actively developed. Variational approaches seek to minimize a cost or penalty function which measures a model's fit to observations, background fields and other imposed constraints. Alternatively, the Kalman filter technique, which is also under investigation as a data assimilation procedure for numerical weather prediction, can yield acceptable initial conditions for mesoscale models. The third kind of next-generation technique involves strategies to initialize convective scale (non-hydrostatic) models.

Warner, Tom↗

The Delta(dot) B = O Constraint vs. Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. On the other hand, the MHD equations can be derived from basic principles in either conservative or non-conservative form. The non-conservative system has a better conditioned eigensystem. The Delta(dot)B = 0 constraint of the A4HD equations is only an initial condition constraint. One does not need the Delta(dot)B condition to close the MHD system. We formulate our new low dissipative high order scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the Delta(dot)B condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗

Simulation of Multiple Shocks in the Solar Atmosphere using a High-Order Shock-Capturing Scheme

We present a simulation of shock phenomena in the solar atmosphere using high-order shock capturing central-upwind schemes, based on the Euler equation in gravity. These simulations use observation-based initial conditions, with results matching observation. Various difficulties are highlighted, including very high gradients in the initial data, multiple scales, and initialization in hydrostatic equilibrium. We focus on the design of algorithms to address these problems, as well as the interaction with observational data.

Bryson, Steve↗

The DeepMIP Contribution to PMIP4: Experimental Design for Model Simulations of the EECO, PETM, and pre-PETM (version 1.0)

Past warm periods provide an opportunity to evaluate climate models under extreme forcing scenarios, in particular high ( greater than 800 ppmv) atmospheric CO2 concentrations. Although a post hoc intercomparison of Eocene (approximately 50 Ma) climate model simulations and geological data has been carried out previously, models of past high-CO2 periods have never been evaluated in a consistent framework. Here, we present an experimental design for climate model simulations of three warm periods within the early Eocene and the latest Paleocene (the EECO, PETM, and pre-PETM). Together with the CMIP6 pre-industrial control and abrupt 4(times) CO2 simulations, and additional sensitivity studies, these form the first phase of DeepMIP - the Deep-time Model Intercomparison Project, itself a group within the wider Paleoclimate Modeling Intercomparison Project (PMIP). The experimental design specifies and provides guidance on boundary conditions associated with palaeogeography, greenhouse gases, astronomical configuration, solar constant, land surface processes, and aerosols. Initial conditions, simulation length, and output variables are also specified. Finally, we explain how the geological data sets, which will be used to evaluate the simulations, will be developed.

climate models↗

An Exploratory Study of the Butterfly Effect Using Agent-Based Modeling

This paper provides insights about the behavior of chaotic complex systems, and the sensitive dependence of the system on the initial starting conditions. How much does a small change in the initial conditions of a complex system affect it in the long term? Do complex systems exhibit what is called the "Butterfly Effect"? This paper uses an agent-based modeling approach to address these questions. An existing model from NetLogo library was extended in order to compare chaotic complex systems with near-identical initial conditions. Results show that small changes in initial starting conditions can have a huge impact on the behavior of chaotic complex systems. The term the "butterfly effect" is attributed to the work of Edward Lorenz [1]. It is used to describe the sensitive dependence of the behavior of chaotic complex systems on the initial conditions of these systems. The metaphor refers to the notion that a butterfly flapping its wings somewhere may cause extreme changes in the ecological system's behavior in the future, such as a hurricane.

Khasawneh, Mahmoud T.↗

Excitation of nearly steady finite-amplitude barotropic waves

An exact nonlinear stationary solution is obtained for barotropic waves in a beta-plane channel and it is shown that it can be excited under a range of initial conditions. Results show that a finite-amplitude wave in a constant shear flow, given an initial phase tilt against the shear and a sufficient initial amplitude, interacts with the mean flow to produce a nearly steady state close to the exact stationary solution. This equilibration process involves nonlinear transients; in particular, as the flow equilibrates, the emergence of critical levels is accompanied by the neutralization of local mean vorticity gradients at these levels, thus allowing the solution to attain a nonsingular modal structure.

Hou, A. Y.↗