Search NASA⌕ Search

SEARCH · Search NASA

Results for “Continuity equation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Rocket Investigation of Mesospheric Eddy Diffusion Effects on Airglow and Oxygen Chemistry

A Terrier Orion rocket was launched at 0750 Z on 02/25/98 about seven minutes after the Clemson University chemical release rocket. Measurements made of the electron density by a dc probe calibrated by a capacitance probe showed several layers of electron density on a rocket ascent in the altitude range from 90 to 110 km. Rocket descent results showed several but not all of the ascent structure. From power spectral analysis of the measured electron densities, turbulent parameters are derived Measurements were made on rocket ascent and descent by an infrared radiometer of the OH Meinel (3-1) band and O2 singlet delta emissions. Profiles of the emissions are presented and discussed on both rocket ascent and descent an enhancement of the OH emission monitored by the OH radiometer was observed above 90 km. The glow was not defected by the O2 radiometer and was significantly reduced on rocket descent. Using these data and a mechanistic analysis, a profile proportional to atomic oxygen is obtained. This profile is compared to one from the ATOX probe on the rocket. A one-dimensional (1-D) photochemical model that solves the time-dependent continuity equations is used with the rocket data to investigate the odd-oxygen concentration in the near equatorial mesosphere.

Ulwick, James C.↗
NASA NTRS2001

Design of Fiber Optic Sensors for Measuring Hydrodynamic Parameters

The science of optical hydrodynamics involves relating the optical properties to the fluid dynamic properties of a hydrodynamic system. Fiber-optic sensors are being designed for measuring the hydrodynamic parameters of various systems. As a flowing fluid makes an encounter with a flat surface, it forms a boundary layer near this surface. The region between the boundary layer and the flat plate contains information about parameters such as viscosity, compressibility, pressure, density, and velocity. An analytical model has been developed for examining the hydrodynamic parameters near the surface of a fiber-optic sensor. An analysis of the conservation of momentum, the continuity equation and the Navier-Stokes equation for compressible flow were used to develop expressions for the velocity and the density as a function of the distance along the flow and above the surface. When examining the flow near the surface, these expressions are used to estimate the sensitivity required to perform direct optical measurements and to derive the shear force for indirect optical measurements. The derivation of this result permits the incorporation of better design parameters for other fiber-based sensors. Future work includes analyzing the optical parametric designs of fiber-optic sensors, modeling sensors to utilize the parameters for hydrodynamics and applying different mixtures of hydrodynamic flow. Finally, the fabrication of fiber-optic sensors for hydrodynamic flow applications of the type described in this presentation could enhance aerospace, submarine, and medical technology.

Lyons, Donald R.↗
NASA NTRS2002

The Fractional Step Method Applied to Simulations of Natural Convective Flows

This paper describes research done to apply the Fractional Step Method to finite-element simulations of natural convective flows in pure liquids, permeable media, and in a directionally solidified metal alloy casting. The Fractional Step Method has been applied commonly to high Reynold's number flow simulations, but is less common for low Reynold's number flows, such as natural convection in liquids and in permeable media. The Fractional Step Method offers increased speed and reduced memory requirements by allowing non-coupled solution of the pressure and the velocity components. The Fractional Step Method has particular benefits for predicting flows in a directionally solidified alloy, since other methods presently employed are not very efficient. Previously, the most suitable method for predicting flows in a directionally solidified binary alloy was the penalty method. The penalty method requires direct matrix solvers, due to the penalty term. The Fractional Step Method allows iterative solution of the finite element stiffness matrices, thereby allowing more efficient solution of the matrices. The Fractional Step Method also lends itself to parallel processing, since the velocity component stiffness matrices can be built and solved independently of each other. The finite-element simulations of a directionally solidified casting are used to predict macrosegregation in directionally solidified castings. In particular, the finite-element simulations predict the existence of 'channels' within the processing mushy zone and subsequently 'freckles' within the fully processed solid, which are known to result from macrosegregation, or what is often referred to as thermo-solutal convection. These freckles cause material property non-uniformities in directionally solidified castings; therefore many of these castings are scrapped. The phenomenon of natural convection in an alloy under-going directional solidification, or thermo-solutal convection, will be explained. The development of the momentum and continuity equations for natural convection in a fluid, a permeable medium, and in a binary alloy undergoing directional solidification will be presented. Finally, results for natural convection in a pure liquid, natural convection in a medium with a constant permeability, and for directional solidification will be presented.

Westra, Douglas G.↗
NASA NTRS2004

Development of a Computationally Efficient, High Fidelity, Finite Element Based Hall Thruster Model

This report documents the development of a two dimensional finite element based numerical model for efficient characterization of the Hall thruster plasma dynamics in the framework of multi-fluid model. Effect of the ionization and the recombination has been included in the present model. Based on the experimental data, a third order polynomial in electron temperature is used to calculate the ionization rate. The neutral dynamics is included only through the neutral continuity equation in the presence of a uniform neutral flow. The electrons are modeled as magnetized and hot, whereas ions are assumed magnetized and cold. The dynamics of Hall thruster is also investigated in the presence of plasma-wall interaction. The plasma-wall interaction is a function of wall potential, which in turn is determined by the secondary electron emission and sputtering yield. The effect of secondary electron emission and sputter yield has been considered simultaneously, Simulation results are interpreted in the light of experimental observations and available numerical solutions in the literature.

Jacobson, David↗
NASA NTRS2008

A Semi-Empirical Model for Forecasting Relativistic Electrons at Geostationary Orbit

We developed a new prediction model for forecasting relativistic (>2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/Interplanetary Magnetic Field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is about 0.9. The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible. The correlation coefficient between predicted and actual electron fluxes is stable and incredibly high.

Lyatsky, Wladislaw↗
NASA NTRS2007

Theory of the Lattice Boltzmann Equation: Symmetry properties of Discrete Velocity Sets

In the lattice Boltzmann equation, continuous particle velocity space is replaced by a finite dimensional discrete set. The number of linearly independent velocity moments in a lattice Boltzmann model cannot exceed the number of discrete velocities. Thus, finite dimensionality introduces linear dependencies among the moments that do not exist in the exact continuous theory. Given a discrete velocity set, it is important to know to exactly what order moments are free of these dependencies. Elementary group theory is applied to the solution of this problem. It is found that by decomposing the velocity set into subsets that transform among themselves under an appropriate symmetry group, it becomes relatively straightforward to assess the behavior of moments in the theory. The construction of some standard two- and three-dimensional models is reviewed from this viewpoint, and procedures for constructing some new higher dimensional models are suggested.

Rubinstein, Robert↗
NASA NTRS2008

Relativistic Electrons at Geostationary Orbit: Modeling Results

We developed a new prediction model for forecasting relativistic (greater than 2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/interplanetary magnetic field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is stable and incredibly high (about 0.9). The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible.

Khazanov, George V.↗
NASA NTRS2008

Prediction Model for Relativistic Electrons at Geostationary Orbit

We developed a new prediction model for forecasting relativistic (greater than 2MeV) electrons, which provides a VERY HIGH correlation between predicted and actually measured electron fluxes at geostationary orbit. This model implies the multi-step particle acceleration and is based on numerical integrating two linked continuity equations for primarily accelerated particles and relativistic electrons. The model includes a source and losses, and used solar wind data as only input parameters. We used the coupling function which is a best-fit combination of solar wind/interplanetary magnetic field parameters, responsible for the generation of geomagnetic activity, as a source. The loss function was derived from experimental data. We tested the model for four year period 2004-2007. The correlation coefficient between predicted and actual values of the electron fluxes for whole four year period as well as for each of these years is stable and incredibly high (about 0.9). The high and stable correlation between the computed and actual electron fluxes shows that the reliable forecasting these electrons at geostationary orbit is possible.

Khazanov, George V.↗
NASA NTRS2015

Hybrid Glacier Inventory, Gravimetry and Altimetry (HIGA) Mass Balance Product for Greenland and the Canadian Arctic

We present a novel inversion algorithm that generates a mass balance field that is simultaneously consistent with independent observations of glacier inventory derived from optical imagery, cryosphere-attributed mass trends derived from satellite gravimetry, and ice surface elevation trends derived from airborne and satellite altimetry. We use this algorithm to assess mass balance across Greenland and the Canadian Arctic over the Sep-2003 to Oct- 2009 period at 26 kilometers resolution. We evaluate local algorithm-inferred mass balance against forty in situ point observations. This evaluation yields a root mean squared error (RMSE) of 0.15 mWE/a ( 0.15 meters (water equivalent) per annum), and highlights a paucity of in situ observations from regions of high dynamic mass loss and peripheral glaciers. We assess mass losses of 212 plus or minus 67 Gigatons per annum to the Greenland ice sheet proper, 38 Gigatons per annum to peripheral glaciers in Greenland, and 42 Gigatons per annum to glaciers in the Canadian Arctic. These magnitudes of mass loss are dependent on the gravimetry-derived spherical harmonic mass trend we invert. We spatially partition the transient glacier continuity equation by differencing algorithm inferred mass balance from modeled surface mass balance, in order to solve the horizontal divergence of ice flux as a residual. This residual ice dynamic field infers flux divergence (or submergent flow) in the ice sheet accumulation area and at tidewater margins, and flux convergence (or emergent flow) in land-terminating ablation areas, which is consistent with continuum mechanics theory.

Altimetry↗
NASA NTRS2015

Electron-Impact Excitation Cross Sections for Modeling Non-Equilibrium Gas

In order to provide a database for modeling hypersonic entry in a partially ionized gas under non-equilibrium, the electron-impact excitation cross sections of atoms have been calculated using perturbation theory. The energy levels covered in the calculation are retrieved from the level list in the HyperRad code. The downstream flow-field is determined by solving a set of continuity equations for each component. The individual structure of each energy level is included. These equations are then complemented by the Euler system of equations. Finally, the radiation field is modeled by solving the radiative transfer equation.

Electron↗
NASA NTRS2016

Advancements in Afterbody Radiative Heating Simulations for Earth Entry

Four advancements to the simulation of backshell radiative heating for Earth entry are presented. The first of these is the development of a flow field model that treats electronic levels of the dominant backshell radiator, N, as individual species. This is shown to allow improvements in the modeling of electron-ion recombination and two-temperature modeling, which are shown to increase backshell radiative heating by 10 to 40%. By computing the electronic state populations of N within the flow field solver, instead of through the quasi-steady state approximation in the radiation code, the coupling of radiative transition rates to the species continuity equations for the levels of N, including the impact of non-local absorption, becomes feasible. Implementation of this additional level of coupling between the flow field and radiation codes represents the second advancement presented in this work, which is shown to increase the backshell radiation by another 10 to 50%. The impact of radiative transition rates due to non-local absorption indicates the importance of accurate radiation transport in the relatively complex flow geometry of the backshell. This motivates the third advancement, which is the development of a ray-tracing radiation transport approach to compute the radiative transition rates and divergence of the radiative flux at every point for coupling to the flow field, therefore allowing the accuracy of the commonly applied tangent-slab approximation to be assessed for radiative source terms. For the sphere considered at lunar-return conditions, the tangent-slab approximation is shown to provide a sufficient level of accuracy for the radiative source terms, even for backshell cases. This is in contrast to the agreement between the two approaches for computing the radiative flux to the surface, which differ by up to 40%. The final advancement presented is the development of a nonequilibrium model for NO radiation, which provides significant backshell radiation at velocities below 10 km/s. The developed model reduces the nonequilibrium NO radiation by 50% relative to the previous model.

Johnston, Christopher O.↗
NASA NTRS2016

Evidence for the Kinematic Sunyaev-Zel'dovich Effect with the Atacama Cosmology Telescope and Velocity Reconstruction from the Baryon Oscillation Spectroscopic Survey

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope at 146 GHz, together with the "Constant Mass" CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Zel'dovich (kSZ) effect over the redshift range z1/4 0.4-0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the microwave temperature at the location of each halo, weighted by the corresponding reconstructed velocity. We vary the size of the aperture photometry filter used, thus probing the free electron profile of these halos from within the virial radius out to three virial radii, on the scales relevant for investigating the missing baryons problem. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3 and 2.9 sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The data suggest that the baryon profile is shallower than the dark matter in the inner regions of the halos probed here, potentially due to energy injection from active galactic nucleus or supernovae. Thus, by constraining the gas profile on a wide range of scales, this technique will be useful for understanding the role of feedback in galaxy groups and clusters. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Zel'dovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

Schaan, Emmanuel S.↗
NASA NTRS2014

Seasonal Evolution of Saturn's Polar Temperatures and Composition

The seasonal evolution of Saturn's polar atmospheric temperatures and hydrocarbon composition is derived from a decade of Cassini Composite Infrared Spectrometer (CIRS) 7-16 micrometers thermal infrared spectroscopy. We construct a near-continuous record of atmospheric variability poleward of 60 deg from northern winter/southern summer (2004, Ls = 293 deg) through the equinox (2009, Ls= 0 deg) to northern spring/southern autumn (2014, Ls = 56 deg). The hot tropospheric polar cyclones that are entrained by pro-grade jets within 2-3 deg of each pole, and the hexagonal shape of the north polar belt, are both persistent features throughout the decade of observations. The hexagon vertices rotated westward by approx. equal to 30 deg longitude between March 2007 and April 2013, confirming that they are not stationary in the Voyager-defined System III longitude system as previously thought. Tropospheric temperature contrasts between the cool polar zones (near 80-85 deg) and warm polar belts (near 75-80 deg) have varied in both hemispheres, resulting in changes to the vertical wind shear on the zonal jets in the upper troposphere and lower stratosphere. The extended region of south polar stratospheric emission has cooled dramatically poleward of the sharp temperature gradient near 75 deg S (by approximately -5 K/yr), coinciding with a depletion in the abundances of acetylene (0030 +/- 0.005 ppm/yr) and ethane (0.35 +/- 0.1 ppm/yr), and suggestive of stratospheric upwelling with vertical wind speeds of w approx. equal to +0.1 mm/s. The upwelling appears most intense within 5 deg latitude of the south pole. This is mirrored by a general warming of the northern polar stratosphere (+5 K/yr) and an enhancement in acetylene (0.030 +/- 0.003 ppm/yr) and ethane (0.45 +/- 0.1 ppm/yr) abundances that appears to be most intense poleward of 75 deg N, suggesting subsidence at w approx. equal to -0.15 mm/ s. However, the sharp gradient in stratospheric emission expected to form near 75 deg N by northern summer solstice (2017, Ls = 90 deg) has not yet been observed, so we continue to await the development of a northern summer stratospheric vortex. The peak stratospheric warming in the north occurs at lower pressure levels (p less than 1 mbar) than the peak stratospheric cooling in the south (p greater than 1 mbar). Vertical motions are derived from both the temperature field (using the measured rates of temperature change and the deviations from the expectations of radiative equilibrium models) and hydrocarbon distributions (solving the continuity equation). Vertical velocities tend towards zero in the upper troposphere where seasonal temperature contrasts are smaller, except within the tropospheric polar cyclones where w approx. equal to +0.02 mm/s. North polar minima in tropospheric and stratospheric temperatures were detected in 2008-2010 (lagging one season, or 6-8 years, behind winter solstice); south polar maxima appear to have occurred before the start of the Cassini observations (1-2 years after summer solstice), consistent with the expectations of radiative climate models. The influence of dynamics implies that the coldest winter temperatures occur in the 75-80 deg region in the stratosphere, and in the cool polar zones in the troposphere, rather than at the poles themselves. In addition to vertical motions, we propose that the UV-absorbent polar stratospheric aerosols entrained within Saturn's vortices contribute significantly to the radiative budget at the poles, adding to the localized enhancement in the south polar cooling and north polar warming poleward of +/-75 deg.

polar temperatures↗
NASA NTRSDate not supplied

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

We demonstrate for state dynamics governed by the continuity equation and related hyperbolic partial differential equations that significant, spurious variance loss occurs during covariance propagation by traditional methods used in data assimilation, even at full rank. This inaccurate variance evolution is caused not by numerical dissipation, but rather by a discontinuous change in the continuum covariance dynamics as correlation lengths tend to zero.

Shay Gilpin↗
NASA NTRSDate not supplied

Prediction of High-Latitude Ionospheric Electrodynamics Using the Machine Learning Based Auroral Ionospheric Electrodynamics Model

We introduce a new framework for Machine-Learning (ML) based Auroral Ionosphere Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents (FACs) of Kunduri et al., 2020 (https://doi.org/10.1029/2020JA027908), the FAC-derived aurora conductance model of Robinson et al., 2020 (https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen & Brekke (1993). The ML-AIM inputs are 60min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pederson/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity on 14 May 2013 and a geomagnetic storm on 7-8 September 2017. ML-AIM produces reasonable ionospheric potential patterns such as two cell convection patterns and the enhancement of electric potentials during active times. The cross polar cap potential drop from ML-AIM is also comparable to the ones from the Weimer 2005 model, Super Dual Auroral Radar Network (SuperDARN), and Defense Meteorological Satellite Program (DMSP) F17 observations. ML-AIM is unique in a sense that it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while other traditional empirical model like Weimer 2005 is designed to provide static ionospheric conditions under steady solar wind/IMF conditions. In future, ML-AIM will include ML-based models of aurora precipitation and ionospheric conductance, improving its performance during active times.

H. K. Connor↗
NASA NTRS2024

Calculating the High-Latitude Ionospheric Electrodynamics Using A Machine Learning-Based Field-Aligned Current Model

We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020, https://doi.org/10.1029/2020JA027908), the FAC-derived auroral conductance model of Robinson et al. (2020, https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993, https://doi.org/10.1029/92gl02109). The ML-AIM inputs are 60-min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML-AIM produces physically accurate ionospheric potential patterns such as the two-cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML-AIM, the Weimer (2005, https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data-assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML-AIM than others. ML-AIM is unique and innovative because it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005, https://doi.org/10.1029/2004ja010884) designed to provide a quasi-static ionospheric condition under quasi-steady solar wind/IMF conditions. Plans are underway to improve ML-AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times.

auroral electrodynamics↗
NASA NTRS1990

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗