Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Weakly multi-dimensional cosmic-ray-modified MHD shocks

The multi-dimensional structure of weak, energetic-particle-modified shocks is investigated by means of appropriate perturbation techniques. The time-dependent shock-structure equation is found to be a generalized form of the well-known one-dimensional Burgers equation, whose steady state, in the absence of cosmic rays, is shown to be related to an equation modeling steady transonic flow in several dimensions. The time-dependent (1 + 2)- and (1 + 3)-dimensional Burgers equations are integrated exactly by means of Hirota's technique for one-shock solutions. On the basis of the exact solutions, a discussion relating the various length scales associated with the shock is presented.

Zank, G. P.↗
NASA NTRS2016

Modeling Global Water Use for the 21st Century: Water Futures and Solutions (WFaS) Initiative and Its Approaches

To sustain growing food demand and increasing standard of living, global water use increased by nearly 6 times during the last 100 years, and continues to grow. As water demands get closer and closer to the water availability in many regions, each drop of water becomes increasingly valuable and water must be managed more efficiently and intensively. However, soaring water use worsens water scarcity conditions already prevalent in semi-arid and arid regions, increasing uncertainty for sustainable food production and economic development. Planning for future development and investments requires that we prepare water projections for the future. However, estimations are complicated because the future of the world's waters will be influenced by a combination of environmental, social, economic, and political factors, and there is only limited knowledge and data available about freshwater resources and how they are being used. The Water Futures and Solutions (WFaS) initiative coordinates its work with other ongoing scenario efforts for the sake of establishing a consistent set of new global water scenarios based on the shared socio-economic pathways (SSPs) and the representative concentration pathways (RCPs). The WFaS "fast track" assessment uses three global water models, namely H08, PCR-GLOBWB, and WaterGAP. This study assesses the state of the art for estimating and projecting water use regionally and globally in a consistent manner. It provides an overview of different approaches, the uncertainty, strengths and weaknesses of the various estimation methods, types of management and policy decisions for which the current estimation methods are useful. We also discuss additional information most needed to be able to improve water use estimates and be able to assess a greater range of management options across the water-energy-climate nexus.

soaring↗
NASA NTRS1979

Influence of a weak gravitational wave on a bound system of two point-masses

The problem of a weak gravitational wave impinging upon a nonrelativistic bound system of two point masses is considered. The geodesic equation for each mass is expanded in terms of two small parameters, v/c and dimensionless wave amplitude, in a manner similar to the post-Newtonian expansion; the geodesic equations are resolved into orbital and center-of-mass equations of motion. The effect of the wave on the orbit is determined by using Lagrange's planetary equations to calculate the time evolution of the orbital elements. The gauge properties of the solutions and, in particular, the gauge invariance of the secular effects are discussed.

Turner, M. S.↗
NASA NTRS1993

On nonlinear Tollmien-Schlichting/vortex interaction in three-dimensional boundary layers

The instability of an incompressible three-dimensional boundary layer (that is, one with cross-flow) is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two low amplitude, lower-branch Tollmien-Schlichting waves which mutually interact to induce a weak longitudinal vortex flow; the vortex motion, in turn, gives rise to significant wave-modulation via wall-shear forcing. The characteristic Reynolds number is taken as a large parameter and, as a consequence, the waves' and the vortex motion are governed primarily by triple-deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Three distinct possibilities were found to emerge for the nonlinear behavior of the flow solution downstream - an algebraic finite-distance singularity, far downstream saturation or far-downstream wave-decay (leaving pure vortex flow) - depending on the input conditions, the wave angles, and the size of the cross-flow.

Davis, Dominic A. R.↗
NASA NTRS2018

Parameterization of Mixed Layer and Deep-Ocean Mesoscales Including Nonlinearity

In 2011, Chelton et al. carried out a comprehensive census of mesoscales using altimetry data and reached the following conclusions: "essentially all of the observed mesoscale features are nonlinear" and "mesoscales do not move with the mean velocity but with their own drift velocity," which is "the most germane of all the nonlinear metrics." Accounting for these results in a mesoscale parameterization presents conceptual and practical challenges since linear analysis is no longer usable and one needs a model of nonlinearity. A mesoscale parameterization is presented that has the following features: 1) it is based on the solutions of the nonlinear mesoscale dynamical equations, 2) it describes arbitrary tracers, 3) it includes adiabatic (A) and diabatic (D) regimes, 4) the eddy-induced velocity is the sum of a Gent and McWilliams (GM) term plus a new term representing the difference between drift and mean velocities, 5) the new term lowers the transfer of mean potential energy to mesoscales, 6) the isopycnal slopes are not as flat as in the GM case, 7) deep-ocean stratification is enhanced compared to previous parameterizations where being more weakly stratified allowed a large heat uptake that is not observed, 8) the strength of the Deacon cell is reduced. The numerical results are from a stand-alone ocean code with Coordinated Ocean-Ice Reference Experiment I (CORE-I) normal-year forcing.

Baroclinic flows; Mesoscale processes; Mixing; Sub↗
NASA NTRS1996

Active Sonic Boom Control

A theory and simulation code are developed to study nonsteady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of flight path. An aircraft in nonsteady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender nonsteady sources but requires no far field approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two nonsteady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and convert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on the average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical

Crow, Steven C.↗
NASA NTRS1993

Variable-permittivity linear inverse problem for the H(sub z)-polarized case

The H(sub z)-polarized inverse problem has rarely been studied before due to the complicated way in which the unknown permittivity appears in the wave equation. This problem is equivalent to the acoustic inverse problem with variable density. We have recently reported the solution to the nonlinear variable-permittivity H(sub z)-polarized inverse problem using the Born iterative method. Here, the linear inverse problem is solved for permittivity (epsilon) and permeability (mu) using a different approach which is an extension of the basic ideas of diffraction tomography (DT). The key to solving this problem is to utilize frequency diversity to obtain the required independent measurements. The receivers are assumed to be in the far field of the object, and plane wave incidence is also assumed. It is assumed that the scatterer is weak, so that the Born approximation can be used to arrive at a relationship between the measured pressure field and two terms related to the spatial Fourier transform of the two unknowns, epsilon and mu. The term involving permeability corresponds to monopole scattering and that for permittivity to dipole scattering. Measurements at several frequencies are used and a least squares problem is solved to reconstruct epsilon and mu. It is observed that the low spatial frequencies in the spectra of epsilon and mu produce inaccuracies in the results. Hence, a regularization method is devised to remove this problem. Several results are shown. Low contrast objects for which the above analysis holds are used to show that good reconstructions are obtained for both permittivity and permeability after regularization is applied.

Moghaddam, M.↗
NASA NTRS1994

Progress in Grid Generation: From Chimera to DRAGON Grids

Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Liou, Meng-Sing↗
NASA NTRS2000

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗
NASA NTRS1983

Orbiting chains and rings

Two very different, highly flexible, space structures proposed during the last decade which involve analytical solution of certain partial differential equations are discussed. The first structure is an array, or hanging chain of aluminum beads which would serve as a convenient communicator if, under the influence of the Earth's gravity gradient, it assumes a local vertical orientation as it circles the Earth. Two passive schemes have been proposed for damping the rigid modes: twist the wire at the ends to provide non-zero moment of inertia about the vertical, thereby inducing relative motion of the two tips during rigid pitch or roll, and thus exercising a damper; and introduce weak lossy springs between the end sections and the main section; thereby providing linear coupling between the springs and all the in-plane (pitch) modes. The second structure is a complete ring of satellites cabled together at synchronous altitude. By a slight increase in altitude, the cable is in tension. This configuration is, however, unstable, and an active feedback control scheme is required to stabilize it. A possible scheme involves local cable length adjustment based on measurement of local altitude and shape variations and their rates.

Breakwell, J. V.↗
NASA NTRS1995

The adaptive, cut-cell Cartesian approach (warts and all)

Solution-adaptive methods based on cutting bodies out of Cartesian grids are gaining popularity now that the ways of circumventing the accuracy problems associated with small cut cells have been developed. Researchers are applying Cartesian-based schemes to a broad class of problems now, and, although there is still development work to be done, it is becoming clearer which problems are best suited to the approach (and which are not). The purpose of this paper is to give a candid assessment, based on applying Cartesian schemes to a variety of problems, of the strengths and weaknesses of the approach as it is currently implemented.

Powell, Kenneth G.↗
NASA NTRS2001

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗
NASA NTRS1987

Expert system development for commonality analysis in space programs

This report is a combination of foundational mathematics and software design. A mathematical model of the Commonality Analysis problem was developed and some important properties discovered. The complexity of the problem is described herein and techniques, both deterministic and heuristic, for reducing that complexity are presented. Weaknesses are pointed out in the existing software (System Commonality Analysis Tool) and several improvements are recommended. It is recommended that: (1) an expert system for guiding the design of new databases be developed; (2) a distributed knowledge base be created and maintained for the purpose of encoding the commonality relationships between design items in commonality databases; (3) a software module be produced which automatically generates commonality alternative sets from commonality databases using the knowledge associated with those databases; and (4) a more complete commonality analysis module be written which is capable of generating any type of feasible solution.

Yeager, Dorian P.↗
NASA NTRS1985

Review of recent laboratory double layer experiments

Recent laboratory measurements of double layers are reviewed. Most experiments that are considered employed triple plasma devices or Q-machines. It is shown that a variety of both one and three dimensional phenomena have been achieved. Stationary one and three dimensional structures consist of monotonic double layers, very weak, weak, strong and very strong double layers with potential steps e(phi)/T(e) equal to about 1, less than, greater than, and much greater than 10, respectively. Multiple double layers consist of structures with large potential dips on the low potential side and stairstep double layers. Two and three dimensional structures include multiple double layers, which resemble U shaped double layers, and ionization produced strong double layers. Moving double layers include ion acoustic double layers and double layers which are the result of drifting species and mismatches at plasma boundaries. It appears that double layers can be BGK solutions or related to turbulence. Double layer creation can occur in a variety of different ways. Ionization can sometimes be important and sometimes play no role at all.

Hershkowitz, N.↗
NASA NTRS1990

The compressible vortex pair

A numerical solution for the flow field associated with a compressible pair of counter-rotating vortices is developed. The compressible, two-dimensional potential equation is solved utilizing the numerical method of Osher et al. (1985) for flow regions in which a non-zero density exists. Close to the vortex centers, vacuum 'cores' develop owing to the existence of a maximum achievable flow speed in a compressible flow field. A special treatment is required to represent these vacuum cores. Typical streamline patterns and core boundaries are obtained for upstream Mach numbers as high as 0.3, and the formation of weak shocks, predicted by Moore and Pullin (1987), is observed.

Heister, S. D.↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS1994

(abstract) Venus Gravity Data Reduction

The Magellan spacecraft has provided high resolution gravity data to its very end, October 13, 1994, when it was consumed by the Venusian atmosphere. After aerobraking in August of 1993 to attain a near circular orbit, excellent high latitude data were acquired which previously were very weak during the elliptical orbit coverage. There are 1500 orbits during the near circular orbit, supplying redundant coverage at different geometries over many features. This allowed the relaxation of apriori constraints, so true amplitudes are being extracted from the data. In this paper we present the results of a 75(sup th) degree and order field that incorporates all the old Pioneer Venus Orbiter data as well as all the Magellan data to September 1994. The new results reflect even higher correlation with topography, higher amplitude values for the highs and lows, and global results that have essentially very little apriori constraint on the solution parameters. We also correlate our new model with the earlier ones based on 60(sup th) and 40(sup th) degree and order presentations.

Magellan Venus gravity orbits topography surface f↗
NASA NTRS1985

Results from observations of the 15 June 1983 occultation by the Neptune system

Observations of eight Neptune occultations from six sites in the southwestern Pacific on June 15, 1983. The data were used to search for evidence of rings around Neptune down to a distance of 0.03 Neptune radii from the planetary surface, but the results were negative. An astrometric analysis of the timings yielded solution for the equatorial radius a(0) of Neptune at 1 microbar pressure and the oblatness e at this level. The results are: a(0) = 25,295 + or - 50 km; e = 0.022 + or - 0.004, from which a value for a(1) (equatorial radius at 1 bar pressure) of 24,830 + or - 100 km; and a rotation period of P = 15 h (+3, -2 h) is derived. These results are based on a recent determination of Neptune's pole position and mass quadruple moment, and are consistent with the hypothesis that Neptune and Uranus have homologous mass distributions, although the constraint on interior structure is a weak one. If the probable uncertainty of the pole position is taken into account, error bars on e and derived quantities should be increased by about 50 percent. The present determinations of Neptune atmospheric temperatures at 1 microbar are consistent with earlier results and when combined with all previous data give an average value of 156 + or - 10 K. There is only slight evidence for any latitude dependence in the temperatures. Profiles with a high signal-to-noise ratio suggest the possible presence of an absorbing layer at altitudes higher than the 1-microbar level.

Hubbard, W. B.↗