Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time Evolution”

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 541 records · Page 30

Dissociative Recombination (DR) and Associative Ionization (AI) Cross Section Calculations for the NO+ + e yields N + O (2D+3P, 2P+3P, and 2D+1D) Reaction for Atmospheric Entry Modeling

During entry of a space craft into Earth atmosphere, the flow surrounding the vehicle becomes partially ionized leading to significant cation and free-electron production. Subsequently, electron impact excitation and dissociation form radiating excited state species that contribute to the heat load on the vehicle. While experimental data on selected total cross sections are available, few experiments address the need for accurate AI cross sections for metastable atomic states at the high temperatures realized in atmospheric re-entry. In order to maintain desired safety margins during atmospheric entry of a space vehicle, chemical reaction models need to accurately account for this process. For the present study, we have computed vibrationally resolved cross sections for the DR of NO+ for electron energies between 0.01 to 10 eV and apply microscopic reversibility to obtain the AI cross sections and rate coefficients. For the DR cross section calculations, we use state-of-the-art MRCI potential energy curves. The resulting adiabatic potential curves are transformed to a diabatic representation, which is used in time-dependent wave packet calculations to describe the nuclear motion of the dissociating cation upon collision with an electron. Based on the wave function evolution in time, these calculations provide T-matrix and cross sections for the DR including recombination into the low energy metastable atomic states. The DR and AI rate coefficients are compared with the available experimental data.

Ewa Papajak↗

NICER Observations Reveal that the X-ray Transient MAXI J1348-630 is a Black Hole X-ray Binary

We studied the outburst evolution and timing properties of the recently discovered X-ray transient MAXI J1348−630 as observed with NICER. We produced the fundamental diagrams commonly used to trace the spectral evolution, and power density spectra to study the fast X-ray variability. The main outburst evolution of MAXI J1348−630 is similar to that commonly observed in blackhole transients. The source evolved from the hard state (HS), through hard- and soft-intermediate states, into the soft state in the outburst rise, and back to the HS in reverse during the outburst decay. At the end of the outburst, MAXI J1348−630 underwent two reflares with peak fluxes approximately one and two orders of magnitude fainter than the main outburst, respectively. During the reflares, the source remained in the HS only, without undergoing any state transitions, which is similar to the so-called ‘failed outbursts’. Different types of quasi-periodic oscillations (QPOs) are observed at different phases of the outburst. Based on our spectral-timing results, we conclude that MAXI J1348−630 is a black hole candidate.

L Zhang↗

Observing Supraglacial Lakes Using Deep Learning and PlanetScope Imagery

Supraglacial lakes (SGL)s result from melt water accumulation in topographic depressions on the surface of glaciers. SGLs primarily affect glacial dynamics through a positive feedback loop in which the albedo-lowering effect of SGLs can escalate surface melt leading to increases in lake extent and depth, amplifying the afore mentioned albedo-lowering effect. The implications of accelerated glacial melt include increased sea level rise and modifications to ocean primary productivity. SGLs are critical indicators of surface melt and its downstream impacts and should be monitored efficiently. In situ observations and measurements of SGLs are time consuming, cost-prohibitive and difficult to scale. Earth observation data and machine learning enable scalable monitoring of SGLs through pattern detection and quantification of lake evolution over time [1]. This work presents a model developed by training a convolutional neural network with imagery and labels from NASA Operation IceBridge and predicting SGLs in high temporal and spatial resolution PlanetScope imagery.

Supraglacial lake↗

Comparing Global Atmospheric CO2 Flux and Transport Models with Remote Sensing (and Other) Observations

We report recent progress derived from comparison of global CO2 flux and transport models with new remote sensing and other sources of CO2 data including those from satellite. The overall objective of this activity is to improve the process models that represent our understanding of the workings of the atmospheric carbon cycle. Model estimates of CO2 surface flux and atmospheric transport processes are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, to provide the basic framework for carbon data assimilation, and ultimately for future projections of carbon-climate interactions. Models can also be used to test consistency within and between CO2 data sets under varying geophysical states. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 2000 through 2009. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at 1x1 degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-3), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to remote sensing observations from TCCON, GOSAT, and AIRS as well as relevant in situ observations. Examples of the influence of key process representations are shown from both forward and inverse model comparisons. We find that the model can resolve much of the synoptic, seasonal, and interannual variability in the observations, although reasons for persistent discrepancies in northern hemisphere vegetation uptake are examined. At this time, we do not find any serious shortcomings in the model transport representation, but this is still the subject of close scrutiny. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in CO2 fluxes and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗

The effect of thicker oceanic crust in the Archaean on the growth of continental crust through time

Present crustal evolution models fail to account for the generation of the large volume of continental crust in the required time intervals. All Archaean plate tectonic models, whether invoking faster spreading rates, similar to today's spreading rates, or longer ridge lengths, essentially propose that continental crust has grown by island arc accretion due to the subduction of oceanic crust. The petrological differences that characterize the Archaean from later terrains result from the subduction of hotter oceanic crust into a hotter mantle. If the oceanic crust was appreciably thicker in the Archaean, as geothermal models would indicate, this thicker crust is surely going to have an effect on tectonic processes. A more valid approach is to compare the possible styles of convergence of thick oceanic crust with modern convergence zones. The best modern analog occurs where thick continental crust is colliding with thick continental crust. Oceanic crustal collision on the scale of the present-day Himalayan continental collision zone may have been a frequent occurrence in the Archaean, resulting in extensive partial melting of the hydrous underthrust oceanic crust to produce voluminous tonalite melts, leaving a depleted stabilized basic residuum. Present-day island arc accretion may not have been the dominant mechanism for the growth of the early Archaean crust.

Wilks, M. E.↗

Evaluating the Capacity of Global CO2 Flux and Atmospheric Transport Models to Incorporate New Satellite Observations

As we enter the new era of satellite remote sensing for CO2 and other carbon cyclerelated quantities, advanced modeling and analysis capabilities are required to fully capitalize on the new observations. Model estimates of CO2 surface flux and atmospheric transport are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, and ultimately for future projections of carbon-climate interactions. For application to current, planned, and future remotely sensed CO2 data, it is desirable that these models are accurate and unbiased at time scales from less than daily to multi-annual and at spatial scales from several kilometers or finer to global. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 1998 through 2006. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at lxi degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-2), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to in situ observations at sites in Northern mid latitudes and the continental tropics. The influence of key process representations is inferred. We find that the model can resolve much of the hourly to synoptic variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The seasonal cycle and its interannual variations generally respond adequately, but discrepancies in the tropics suggest the need for a refinement of the soil moisture dependence of the respiration flux in CASA. Examples and inferences for interpretation of satellite data will be discussed. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in the terrestrial CO2 sink and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗

Assimilation of DAWN Doppler Wind Lidar Data During the 2017 Convective Processes Experiment (CPEX): Impact on Precipitation and Flow Structure

An improved representation of 3-D air motion and precipitation structure through forecast models and assimilation of observations is vital for improvements in weather forecasting capabilities. However, there are few independent data to properly validate a model forecast of precipitation structure when the underlying dynamics are evolving on short convective timescales. Using data from the JPL Ku/Ka-band Airborne Precipitation Radar (APR-2) and the 2 μmDoppler Aerosol Wind (DAWN) lidar collected during the2017 Convective Processes Experiment (CPEX), the NASA Unified Weather Research and Forecasting (WRF) Ensemble Data Assimilation System (EDAS) modeling system was used to quantify the impact of high-resolution sparsely sampled DAWN measurements on the analyzed variables and on the forecast when the DAWN winds were assimilated. Over-all, the assimilation of the DAWN wind profiles had a discernible impact on the wind field as well as the evolution and timing of the 3-D precipitation structure. Analysis of individual variables revealed that the assimilation of the DAWN winds resulted in important and coherent modifications of the environment. It led to an increase in the near-surface convergence, temperature, and water vapor, creating more favorable conditions for the development of convection exactly where it was observed (but not present in the control run). Comparison to APR-2 and observations by the Global Precipitation Measurement (GPM) satellite shows a much-improved forecast after the assimilation of the DAWN winds – development of precipitation where there was none, more organized precipitation where there was some, and a much more intense and organized cold pool, similar to the analysis of the dropsonde data. The onset of the vertical evolution of the precipitation showed similar radar-derived cloud-top heights, but delayed in time. While this investigation was limited to a single CPEX flight date, the investigation design is appropriate for further investigation of the impact of airborne Doppler wind lidar observations upon short-term convective precipitation forecasts

DAWN↗

Gasdynamic models of the solar wind/interstellar medium interaction

The interaction between the solar wind and the interstellar medium is modeled self-consistently using numerical solutions of the time-dependent gasdynamic equations in spherical and cylindrical coordinates. For the results presented here it is assumed that the solar system moves through the surrounding medium with a supersonic velocity. After an initial (nonequilibrium) state has been specified, the numerical solution follows the evolution in time until the interaction relaxes to a dynamic equilibrium. As would be expected, the solutions show the formation of a bow shock upstream of the traveling solar system to deflect the interstellar plasma around the cavity created by the solar wind. A terminiation shock also forms to slow and compress the solar wind plasma. For the simulation in spherical coordinates, the downstream portion of the termination shock reaches equilibrium more than three times further from the Sun than the equilibrium distance to the termination shock on the upstream side.

Steinolfson, R. S.↗

Dynamical behavior of surface tension on rotating fluids in low and microgravity environments

Consideration is given to the time-dependent evolutions of the free surface profile (bubble shapes) of a cylindrical container, partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry in low and microgravity environments. The dynamics of the bubble shapes are calculated for four cases: linear time-dependent functions of spin-up and spin-down in low and microgravity, linear time-dependent functions of increasing and decreasing gravity at high and low rotating cylinder speeds, time-dependent step functions of spin-up and spin-down in low gravity, and sinusoidal function oscillation of the gravity environment in high and low rotating cylinder speeds. It is shown that the computer algorithms developed by Hung et al. (1988) may be used to simulate the profile of time-dependent bubble shapes under variations of centrifugal, capillary, and gravity forces.

Hung, R. J.↗

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗

A Precessing Be Disc as a Possible Model for Occultation Events in GX 304-1

We report on the RXTE detection of a sudden increase in the absorption column density, NH, during the 2011 May outburst of GX 304−1. The NH increased up to ∼16 × 10(exp 22) atoms cm(exp −2), which is a factor of 3–4 larger than what is usually measured during the outbursts of GX 304−1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-Proportional Counter Array light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disc around the Be-type companion, we are able to qualitatively explain the NH evolution over time. We are able to constrain the Be disc density to be of the order of 10(exp −11)g cm(exp −3). Our model strengthens the idea of inclined Be discs as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

Kühnel, M.↗

Compilation and Analysis of 20 and 30 GHz Rain Fade Events at the ACTS NASA Ground Station: Statistics and Model Assessment

The purpose of the propagation studies within the ACTS Project Office is to acquire 20 and 30 GHz rain fade statistics using the ACTS beacon links received at the NGS (NASA Ground Station) in Cleveland. Other than the raw, statistically unprocessed rain fade events that occur in real time, relevant rain fade statistics derived from such events are the cumulative rain fade statistics as well as fade duration statistics (beyond given fade thresholds) over monthly and yearly time intervals. Concurrent with the data logging exercise, monthly maximum rainfall levels recorded at the US Weather Service at Hopkins Airport are appended to the database to facilitate comparison of observed fade statistics with those predicted by the ACTS Rain Attenuation Model. Also, the raw fade data will be in a format, complete with documentation, for use by other investigators who require realistic fade event evolution in time for simulation purposes or further analysis for comparisons with other rain fade prediction models, etc. The raw time series data from the 20 and 30 GHz beacon signals is purged of non relevant data intervals where no rain fading has occurred. All other data intervals which contain rain fade events are archived with the accompanying time stamps. The definition of just what constitutes a rain fade event will be discussed later. The archived data serves two purposes. First, all rain fade event data is recombined into a contiguous data series every month and every year; this will represent an uninterrupted record of the actual (i.e., not statistically processed) temporal evolution of rain fade at 20 and 30 GHz at the location of the NGS. The second purpose of the data in such a format is to enable a statistical analysis of prevailing propagation parameters such as cumulative distributions of attenuation on a monthly and yearly basis as well as fade duration probabilities below given fade thresholds, also on a monthly and yearly basis. In addition, various subsidiary statistics such as attenuation rate probabilities are derived. The purged raw rain fade data as well as the results of the analyzed data will be made available for use by parties in the private sector upon their request. The process which will be followed in this dissemination is outlined in this paper.

Manning, Robert M.↗

The Effect on the Lunar Exosphere of a Coroual Mass Ejection Passage

Solar wind bombardment onto exposed surfaces in the solar system produces an energetic component to the exospheres about those bodies. The solar wind energy and composition are highly dependent on the origin of the plasma. Using the measured composition of the slow wind, fast wind, solar energetic particle (SEP) population, and coronal mass ejection (CME), broken down into their various components, we have estimated the total sputter yield for each type of solar wind. We show that the heavy ion component, especially the He++ and 0+7 can greatly enhance the total sputter yield during times when the heavy ion population is enhanced. Folding in the flux, we compute the source rate for several species during different types of solar wind. Finally, we use a Monte Carlo model developed to simulate the time-dependent evolution of the lunar exosphere to study the sputtering component of the exosphere under the influence of a CME passage. We simulate the background exosphere of Na, K, Ca, and Mg. Simulations indicate that sputtering increases the mass of those constituents in the exosphere a few to a few tens times the background values. The escalation of atmospheric density occurs within an hour of onset The decrease in atmospheric density after the CME passage is also rapid, although takes longer than the increase, Sputtered neutral particles have a high probability of escaping the moon,by both Jeans escape and photo ionization. Density and spatial distribution of the exosphere can be tested with the LADEE mission.

Killen, R. M.↗

Dynamics of surface tension in microgravity environment

Time-dependent evolutions of the profile of free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, in low and microgravity environments, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out (1) linear time-dependent functions of spin-up and spin-down in low and microgravity environments, (2) linear time-dependent functions of increasing and decreasing gravity environment at high and low rotating cylinder speeds, (3) time-dependent step functions of spin-up and spin-down in a low-gravity environment, and (4) sinusoidal function oscillation of the gravity environment at high and low rotating cylinder speeds.

Hung, R. J.↗

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗

Tangle-Free Finite Element Mesh Motion for Ablation Problems

Mesh motion is the process by which a computational domain is updated in time to reflect physical changes in the material the domain represents. Such a technique is needed in the study of the thermal response of ablative materials, which erode when strong heating is applied to the boundary. Traditionally, the thermal solver is coupled with a linear elastic or biharmonic system whose sole purpose is to update mesh node locations in response to altering boundary heating. Simple mesh motion algorithms rely on boundary surface normals. In such schemes, evolution in time will eventually cause the mesh to intersect and "tangle" with itself, causing failure. Furthermore, such schemes are greatly limited in the problems geometries on which they will be successful. This paper presents a comprehensive and sophisticated scheme that tailors the directions of motion based on context. By choosing directions for each node smartly, the inevitable tangle can be completely avoided and mesh motion on complex geometries can be modeled accurately.

Droba, Justin↗

Transition and mixing in axisymmetric jets and vortex rings

A class of impulsively started, axisymmetric, laminar jets produced by a time dependent joint source of momentum are considered. These jets are different flows, each initially at rest in an unbounded fluid. The study is conducted at three levels of detail. First, a generalized set of analytic creeping flow solutions are derived with a method of flow classification. Second, from this set, three specific creeping flow solutions are studied in detail: the vortex ring, the round jet, and the ramp jet. This study involves derivation of vorticity, stream function, entrainment diagrams, and evolution of time lines through computer animation. From entrainment diagrams, critical points are derived and analyzed. The flow geometry is dictated by the properties and location of critical points which undergo bifurcation and topological transformation (a form of transition) with changing Reynolds number. Transition Reynolds numbers were calculated. A state space trajectory was derived describing the topological behavior of these critical points. This state space derivation yielded three states of motion which are universal for all axisymmetric jets. Third, the axisymmetric round jet is solved numerically using the unsteady laminar Navier Stokes equations. These equations were shown to be self similar for the round jet. Numerical calculations were performed up to a Reynolds number of 30 for a 60x60 point mesh. Animations generated from numerical solution showed each of the three states of motion for the round jet, including the Re = 30 case.

Allen, G. A., Jr.↗