Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical dispersion”

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 199 records · Page 11

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗

Evolution of Local Microstructures: Spatial Instabilities of Coarsening Clusters

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are often microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is diffusion-limited and arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. Competitive kinetics of these processes dictates that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub V)=0 of the dispersed phase. Since the publication of TLS theory, there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that predicted by TLS is proportional to V(exp 1/2)(sub v) whereas others suggest V(exp 1/3)(sub v).This issue has been addressed recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team.

Glicksman, Martin E.↗

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗

NASA-JSC Protocol for the Characterization of Single Wall Carbon Nanotube Material Quality

It is well known that the raw as well as purified single wall carbon nanotube (SWCNT) material always contain certain amount of impurities of varying composition (mostly metal catalyst and non-tubular carbon). Particular purification method also creates defects and/or functional groups in the SWCNT material and therefore affects the its dispersability in solvents (important to subsequent application development). A number of analytical characterization tools have been used successfully in the past years to assess various properties of nanotube materials, but lack of standards makes it difficult to compare these measurements across the board. In this work we report the protocol developed at NASA-JSC which standardizes measurements using TEM, SEM, TGA, Raman and UV-Vis-NIR absorption techniques. Numerical measures are established for parameters such as metal content, homogeneity, thermal stability and dispersability, to allow easy comparison of SWCNT materials. We will also report on the recent progress in quantitative measurement of non-tubular carbon impurities and a possible purity standard for SWCNT materials.

Arepalli, Sivaram↗

Numerical experiments on the stability of disklike galaxies

It is noted that the stability properties of available galactic models are not consistent with presumptions based on observation and that axisymmetric disk systems with velocity dispersions like those of the Galaxy display major changes in form on a dynamical time scale. A report is given on a series of numerical experiments carried out as part of a systematic search for purely self-consistent disk galaxy models that might undergo little change over the time of several galactic rotations. The stability problem is reviewed, the method and calculations used are described, and initial conditions are outlined. The problem of estimating growth rates is discussed, and growth rates are evaluated for various disturbances. Experimental results are presented concerning growths of axisymmetric disturbances in 'cold' systems as well as disturbances with m = 1, 2, and 3. It is found that m = 2 disturbances are the most virulent and should dominate when all m values are allowed, but that the most rapidly growing disturbances are inhibited when several m values are allowed.

Miller, R. H.↗

The interstellar line spectra of zeta Ophiuchi and zeta Persei and their relation to the short wavelength microwave background radiation

Thirty-one high dispersion Coude spectrograms of zeta Ophiuchi and seven of zeta Persei were numerically synthesized to produce high resolution, low noise spectra in the interval 3650 A to 4350 that yield data on atomic and molecular absorption in well-defined regions of the interstellar medium. The detection threshold is improved by as much as a factor 5 over single plates. Several interstellar lines were discovered in the zeta Oph - 15km/sec cloud and the zeta Per + 13 km/sec cloud.

Bortolot, V. J., Jr.↗

Propagation of plane waves in flow through a variable area duct

Duct geometry is an important factor that influences the transmission of sound in a duct. Nonconstant area produces variations in steady flow quantities through the duct, which can cause reflections of acoustic disturbances and the creation of standing wave patterns, and which can attenuate and disperse propagating waves. To study the problem, a second-order-accurate numerical method has been developed. The method of characteristics is used to solve the acoustic equations with the assumption of quasi-one-dimensional flow in the duct. Numerical results are presented over a broad range of frequencies. Low-frequency results are compared with approximate solutions for compact nozzles and high-frequency calculations are compared with results from the short-wave theory. The numerical results are also compared with an analytical solution of the acoustic equations for exponential horn sections with zero mean flow to further establish the validity of the numerical method.

King, L. S.↗

Gauged Q-balls

Classical non-topological soliton configurations are considered within the theory of a complex scalar field with a gauged U symmetry. Their existence and stability against dispersion are demonstrated and some of their properties are investigated analytically and numerically. The soliton configuration is such that inside the soliton the local U symmetry is broken, the gauge field becomes massive and for a range of values of the coupling constants the soliton becomes a superconductor pushing the charge to the surface. Furthermore, because of the repulsive Coulomb force, there is a maximum size for these objects, making impossible the existence of Q-matter in bulk form. Also briefly discussed are solitons with fermions in a U gauge theory.

Lee, Kimyeong↗

The stability of a compressible stratified shear layer

The stability of a shear layer under the effect of gravity is investigated using the compressible magnetohydrodynamic (MHD) equations, including an effective gravity term to represent the curvature effects of the flow and magnetic field line geometry. A general eigenmode equation is derived for a two-dimensional MHD fluid, and an energy-principle analysis to explain the effect of compressibility on the critical Richardson number is presented. For the case of a hyperbolic tangent shear flow and exponential density profile, it was found that, in the Boussinesq approximation, the compressibility raises the critical Richardson number from 1/4 to as much as 1/2, with the exact value depending on the value of the magnetic field at infinity. Under approximation of a strong asymptotic magnetic field, without invoking the Boussinesq approximation, it is shown both analytically and numerically that the density gradient terms cause the shear instability to be dispersive. The long-wavelength stability boundary for the Richardson number J = 0 is characterized by a normalized phase velocity c =

Wang, Z.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: High beta

The nonlinear dynamics following saturation of the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation in the case of high plasma beta and no wave dispersion. The magnetohydrodynamic (MHD) code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Turbulent cascades develop after saturation of two coupled oblique three-wave parametric instabilities; one of which is an oblique filamentationlike instability reported earlier. Remnants of the parametric processes, as well as of the original Alfven pump wave, persist during late nonlinear times. Nearly incompressible MHD features such as spectral anisotropies appear as well.

Ghosh, S.↗

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: Low beta

The nature of turbulent cascades arising from the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation for the case of low plasma Beta and no wave dispersion. The magnetohydrodynamic code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Anisotropic turbulent cascades, similar to those found in early incompressible two-dimensional simulations, occur after nonlinear saturation of the parallel propagating decay instability. The turbulent spectrum can be divided into three regimes: the lowest wave numbers are dominated by lower sideband remnants of the parametric process, intermediate wave numbers display nearly incompressible dynamics, and the highest wave numbers are dominated by acoustic turbulence.

Ghosh, S.↗

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

Absolute tracer dye concentration using airborne laser-induced water Raman backscatter

The use of simultaneous airborne-laser-induced dye fluorescence and water Raman backscatter to measure the absolute concentration of an ocean-dispersed tracer dye is discussed. Theoretical considerations of the calculation of dye concentration by the numerical comparison of airborne laser-induced fluorescence spectra with laboratory spectra for known dye concentrations using the 3400/cm OH-stretch water Raman scatter as a calibration signal are presented which show that minimum errors are obtained and no data concerning water mass transmission properties are required when the laser wavelength is chosen to yield a Raman signal near the dye emission band. Results of field experiments conducted with an airborne conical scan lidar over a site in New York Bight into which rhodamine dye had been injected in a study of oil spill dispersion are then indicated which resulted in a contour map of dye concentrations, with a minimum detectable dye concentration of approximately 2 ppb by weight.

Hoge, F. E.↗

Status of aerial applications research in the Langley vortex research facility and the Langley full-scale wind tunnel

Small scale models of agricultural airplanes were tested and numerical methods were utilized to study interactions between the airplane wake and the dispersed spray and granular materials. Methods were developed to measure and predict dispersal transport and wake characteristics and dispersal techniques to obtain interactions more favorable to wide, uniform deposition patterns and reduced drift. In the full scale wind tunnel, full scale agricultural airplanes and dispersal systems for both liquid and solid applications were evaluated to improve aircraft aerodynamics and dispersal systems efficiency. The program status in these two facilities is reported with emphasis on wake interactions and dispersal systems research.

Jordan, F. L., Jr.↗

Ion-acoustic instabilities driven by an ion velocity ring

A ring distribution of ions in velocity space can generate electrostatic waves which propagate predominantly along an ambient magnetic field at frequencies comparable with the ion plasma frequency. A dispersion equation which accounts for these waves is presented, and solved analytically and numerically. It was found that ion-acoustic-like waves are excited in a plasma even if the electron temperature is comparable with the in temperature under the assumption of an anisotropic ion distribution.

Akimoto, K.↗