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

Distributed reacceleration of cosmic rays

We develop a model in which cosmic rays, in addition to their initial acceleration by a strong shock, are continuously reaccelerated (e.g., by weak shocks) while propagating through the galaxy. The equations describing this acceleration scheme are solved analytically (approximating ionization losses by a cutoff) and numerically. Solutions for the spectra of primary and secondary cosmic rays are given in a closed analytic form, and they allow a rapid search in parameter space for viable propagation models with distributed reacceleration included. The observed boron-to-carbon ratio can be reproduced by the reacceleration theory over a range of escape parameters, some of them quite different from the standard "leaky box" model. It is also shown that even a very modest amount of reacceleration by strong shocks causes the boron-to carbon ratio to level off at sufficiently high energies, and this effect may be observed in the CRNE data. Several other curiosities in the data may be explained naturally if a modest amount of distributed reacceleration is invoked, including (a) the apparent truncation at low energy in the otherwise exponential pathlength distribution associated with the leaky box model, (b) the sub-iron isotopic anomalies and other effects noted by Silberberg et al., and (c) the discrepancy between the reported 10Be lifetime and the lifetime of cosmic rays in the dense strata of the galactic disk.

NASA Discipline Radiation Health↗
NASA NTRS2008

Chromium and Tantalum Site Substitution Patterns in Ni3Al (L1(sub 2))gamma(prime)- Precipitates

The site substitution behavior of Cr and Ta in the Ni3Al (Ll2)-type gamma'-precipitates of a Ni-Al-Cr-Ta alloy is investigated by atom-probe tomography (APT) and first-principles calculations. Measurements of the gamma'-phase composition by APT suggest that Al, Cr, and Ta share the Al sublattice sites of the gamma'-precipitates. The calculated substitutional energies of the solute atoms at the Ni and Al sublattice sites indicate that Ta has a strong preference for the Al sites, while Cr has a weak Al site preference. Furthermore, Ta is shown to replace Cr at the Al sublattice sites of the gamma'-precipitates, altering the elemental phase partitioning behavior of the Ni-Al-Cr-Ta alloy.

Booth-Morrison, Christopher↗
NASA NTRS2009

A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners

A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state-of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they should behave as resonant liners while the construction of others suggest that they should behave as broadband attenuators. In each case the finite element theory is observed to predict the proper attenuation trend.

Watson, W. R.↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗
NASA NTRS1975

Note on the electric splitting of drift shells

A nonrelativistic analysis is presented for the deviation of a particle from its dipole magnetic drift shell in the presence of an electric field. An alternative definition of weak shell splitting is introduced in which the corotation drift velocity need not be small with respect to the azimuthal magnetic drift velocity. Approximate explicit solutions are obtained for the asymmetry of drift shells in the strong shell splitting case, i.e., for the case where the cross-magnetosphere potential drop is not small compared to particle kinetic energy. Upper limits to drift shell asymmetry are obtained from considerations of the Alfven layer.

Kivelson, M. G.↗
NASA NTRS1990

The distortion of the magnetotail equilibrium structure by a net cross-tail magnetic field

Explicit solutions are presented for quiet, quasi-static magnetotail configurations in the presence of a net cross-tail magnetic field component in addition to field components associated with tail flaring. The solutions are derived using the general three-dimensional equilibrium theory of Birn (1987). It is found that the net By(N) field varies only weakly along field lines and across the tail. The presence of the net By(N) in the neutral sheet implies that a part of the cross-tail current becomes field-aligned, flowing from the southern to the northern hemisphere for By(N) larger than zero and from north to south for By(n) less than zero. The deformations of magnetic flux surfaces suggest that there may be a class of field lines extending from the earth into the magnetosheath region through the low-latitude flanks of the tail. A topological connection is found of the low-latitude region near the magnetopause with the plasma sheet boundary layer, separating closed plasma sheet and open lobe field lines.

Birn, J.↗
NASA NTRS1991

Single-agent parallel window search

Parallel window search is applied to single-agent problems by having different processes simultaneously perform iterations of Iterative-Deepening-A(asterisk) (IDA-asterisk) on the same problem but with different cost thresholds. This approach is limited by the time to perform the goal iteration. To overcome this disadvantage, the authors consider node ordering. They discuss how global node ordering by minimum h among nodes with equal f = g + h values can reduce the time complexity of serial IDA-asterisk by reducing the time to perform the iterations prior to the goal iteration. Finally, the two ideas of parallel window search and node ordering are combined to eliminate the weaknesses of each approach while retaining the strengths. The resulting approach, called simply parallel window search, can be used to find a near-optimal solution quickly, improve the solution until it is optimal, and then finally guarantee optimality, depending on the amount of time available.

Powley, Curt↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗
NASA NTRS1974

Hydrodynamic solutions at a sonic-boom focus

Study of the focusing of sonic-boom N waves using a numerical shock-following method. The computer code used in the calculations employs a moving mesh so that weak N wave shocks can be followed for large distances without numerical dissipation. In several sample calculations, overpressures at the focus ranging from 2 to 20 times the nominal overpressure have been computed. It is also found that weak shocks tend to cusp at the axis in accord with the geometric-acoustics description, whereas strong shocks straighten out without cusping.

Parker, L. W.↗
NASA NTRS1990

Inverse methods for assessing ship-of-opportunity networks and estimating circulation and winds from tropical expendable bathythermograph data

Inverse methods for estimating the surface ciculation of the equatorial Pacific by combining a linear reduced-gravity shallow-water model with the Tropical Ocean-Global Atmosphere ship-of-opportunity expendable bathythermograph (TOGA SOP XBT) observing program are examined. It is demonstrated that a simple linear model of the upper circulation of the equatorial Pacific can be successfully used as a weak constraint when smoothing the TOGA SOP XBT data. A circulation is sought as the weighted least squares fit to the dynamics and the data. The solution method is an expansion in representer functions, and the generalized inverse problem is thereby reduced from a functional problem to an algebraic problem for the coefficients of the representer. A specific inverse calculation using synthetic forcing and data is presented.

Bennett, Andrew F.↗
NASA NTRS1996

Damped Oscillator with Delta-Kicked Frequency

Exact solutions of the Schrodinger equation for quantum damped oscillator subject to frequency delta-kick describing squeezed states are obtained. The cases of strong, intermediate, and weak damping are investigated.

Manko, O. V.↗
NASA NTRS1988

Management of hypophosphatemia

The etiology, clinical presentation, and management of hypophosphatemia are reviewed. Phosphorus is a major intracellular anion and plays an important role in many biochemical pathways relating to normal physiologic functions. Approximately 60 to 90% of the 1 to 1.5 g of daily dietary phosphorus intake is absorbed, and of that amount, about two thirds is excreted in the urine. The overall incidence of hypophosphatemia is about 2 to 3% of all hospitalized patients. Factors associated with hypophosphatemia include phosphate-binding antacid therapy, nasogastric suction, liver disease, sepsis, alcoholism, and acidosis associated with diabetic ketoacidosis. Patients receiving parenteral nutrient solutions were also at higher risk for hypophosphatemia before the routine supplementation of these formulations with phosphate. Patients with hypophosphatemia may be asymptomatic or may experience weakness, malaise, anorexia, bone pain, and respiratory arrest. The major systems involved include the neuromuscular, hematologic, and skeletal systems. Phosphorus-containing products used to treat hypophosphatemia are a combination of monobasic and dibasic phosphate salts. Therefore, it is essential to calculate doses in millimoles rather than milligrams or milliequivalents to more accurately reflect the phosphorus concentration and to avoid potentially serious dosage errors. Normal daily requirements are readily maintained by dietary sources of phosphorus such as milk products or may be supplemented by phosphate-containing products administered orally or intravenously. Since phosphorus is a key factor in many organ systems, it is essential to monitor serum phosphorus concentrations in patients at risk for hypophosphatemia.

Review, Tutorial↗
NASA NTRS2021

Gravitational analog of the canonical acoustic black hole in Einstein-scalar-Gauss-Bonnet theory

In this work, in the context of modified gravity, a curved spacetime analogous to the 'canonical acoustic black hole (CABH)' is constructed. The source is a self-interacting scalar field which is non-minimally coupled to gravity through the Gauss–Bonnet invariant. The scalar-Gauss–Bonnet coupling function is characterized by three positive parameters: σ with units of (length), μ with units of (length)4, and a dimensionless parameter s, thus defining a three-parameter model for which the line element of CABH is a solution. The spacetime is equipped with spherical and static symmetry and has a single horizon determined in Schwarzschild coordinates by the region r = μ^(1/4). The solution admits a photon sphere at r = (3μ)^(1/4), and it is shown that in the region (3μ)^(1/4) ⩽ r < ∞ the scalar field satisfies the null, weak, and strong energy conditions. Nonetheless, the model with s = 1 has major physical relevance since for this case the scalar field is well defined in the entire region r ⩾ μ^(1/4), while for s ≠ 1 the scalar field blows up on the horizon.

Pedro Cañate↗
NASA NTRS1972

On stochastic control system design methods for weakly coupled large scale linear systems.

This paper considers the problem of decentralized control of two weakly coupled linear stochastic systems, using quadratic performance indices. The basic idea is to have each controller control independently his own system, based upon noisy measurements of his own output. To compensate for the effects of weak coupling upon the resultant performance, fake white plant noise is introduced to each system. The appropriate intensity of the fake plant noise is obtained through the solution of an off-line deterministic matrix optimal control problem. The effects of this design method upon the overall coupled system performance are analyzed as a function of the degree of intersystem coupling.

Kwong, R.↗
NASA NTRS1978

Numerical simulation of MHD shock waves in the solar wind

The effects of the interplanetary magnetic field on the propagation speed of shock waves through an ambient solar wind are examined by numerical solutions of the time-dependent nonlinear equations of motion. The magnetic field always increases the velocity of strong shocks. Although the field may temporarily slow down weak shocks inside 1 AU, it eventually also causes weak shocks to travel faster than they would without the magnetic field at larger distances. Consistent with the increase in the shock velocity, the gas pressure ratio across a shock is reduced considerably in the presence of the magnetic field. The numerical method is used to simulate (starting at 0.3 AU) the large deceleration of a shock observed in the lower corona by ground-based radio instrumentation and the more gradual deceleration of the shock in the solar wind observed by the Pioneer 9 and Pioneer 10 spacecraft.

Steinolfson, R. S.↗
NASA NTRS2003

Nonlinear Aeroacoustics Computations by the Space-Time CE/SE Method

The Space-Time Conservation Element and Solution Element Method, or CE/SE Method for short, is a recently developed numerical method for conservation laws. Despite its second order accuracy in space and time, it possesses low dispersion errors and low dissipation. The method is robust enough to cover a wide range of compressible flows: from weak linear acoustic waves to strong discontinuous waves (shocks). An outstanding feature of the CE/SE scheme is its truly multi-dimensional, simple but effective non-reflecting boundary condition (NRBC), which is particularly valuable for computational aeroacoustics (CAA). In nature, the 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 careful treatment of the surface fluxes and geometry, it is different from the existing schemes. Currently, the CE/SE scheme has been developed to a matured stage that a 3-D unstructured CE/SE Navier-Stokes solver is already available. However, in the present review paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen and sketched in section 2. Then applications of the 2-D and 3-D CE/SE schemes to linear, and in particular, nonlinear aeroacoustics are depicted in sections 3, 4, and 5 to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS1971

Optimal startup control of a jacketed tubular reactor.

The optimal startup policy of a jacketed tubular reactor, in which a first-order, reversible, exothermic reaction takes place, is presented. A distributed maximum principle is presented for determining weak necessary conditions for optimality of a diffusional distributed parameter system. A numerical technique is developed for practical implementation of the distributed maximum principle. This involves the sequential solution of the state and adjoint equations, in conjunction with a functional gradient technique for iteratively improving the control function.

Hahn, D. R.↗
NASA NTRS1979

An evaluation of four single element airfoil analytic methods

A comparison of four computer codes for the analysis of two-dimensional single element airfoil sections is presented for three classes of section geometries. Two of the computer codes utilize vortex singularities methods to obtain the potential flow solution. The other two codes solve the full inviscid potential flow equation using finite differencing techniques, allowing results to be obtained for transonic flow about an airfoil including weak shocks. Each program incorporates boundary layer routines for computing the boundary layer displacement thickness and boundary layer effects on aerodynamic coefficients. Computational results are given for a symmetrical section represented by an NACA 0012 profile, a conventional section illustrated by an NACA 65A413 profile, and a supercritical type section for general aviation applications typified by a NASA LS(1)-0413 section. The four codes are compared and contrasted in the areas of method of approach, range of applicability, agreement among each other and with experiment, individual advantages and disadvantages, computer run times and memory requirements, and operational idiosyncrasies.

Freuler, R. J.↗