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 415 records · Page 23

A novel multireceiver communications system configuration based on optimal estimation theory

A multireceiver configuration for the purpose of carrier arraying and/or signal arraying is presented. Such a problem arises for example, in the NASA Deep Space Network where the same data-modulated signal from a spacecraft is received by a number of geographically separated antennas and the data detection must be efficiently performed on the basis of the various received signals. The proposed configuration is arrived at by formulating the carrier and/or signal arraying problem as an optimal estimation problem. Two specific solutions are proposed. The first solution is to simultaneously and optimally estimate the various phase processes received at different receivers with coupled phase locked loops (PLLs) wherein the individual PLLs acquire and track their respective receivers' phase processes, but are aided by each other in an optimal manner. However, when the phase processes are relatively weakly correlated, and for the case of relatively high values of symbol energy-to-noise spectral density ratio, a novel configuration for combining the data modulated, loop-output signals is proposed. The scheme can be extended to the case of low symbol energy-to-noise case by performing the combining/detection process over a multisymbol period. Such a configuration results in the minimization of the effective radio loss at the combiner output, and thus a maximization of energy per bit to noise-power spectral density ration is achieved.

Kumar, R.↗

Nonplanar interface morphologies during unidirectional solidification of a binary alloy. II - Three-dimensional computations

A finite difference method is used to obtain three-dimensional steady-state solutions for nonplanar interface morphologies in order to study the situation of equal thermal properties in the crystal and melt with negligible latent heat release. Stable steady-state solutions corresponding to two-dimensional bands and three-dimensional hexagonal nodes, as well as to rectangular interface planiforms, are found using a model of an aluminum-chromium alloy with a distribution coefficient of greater than one. Hexagonal nodes are predicted near the onset of instability, in agreement with weakly nonlinear theory.

Mcfadden, G. B.↗

Effect of an arcjet plume on satellite reflector performance

The effect of an arcjet plume on the performance of satellite reflector antennas is studied. The arcjet plume is modeled as a weakly ionized plasma. The spatial permittivity distribution of the plume is approximated using the measured electron density profile and a cold plasma model. Geometrical optics is applied to determine the ray paths as well as the transmitted fields through the inhomogeneous plume. The ray optics results are compared against several exact solutions for scattering from inhomogeneous dielectrics, and good agreement is observed for sufficiently large scatterer size. The far-field antenna patterns of the reflector in the presence of the plume are calculated from the transmitted ray fields using a ray-tube integration scheme. For arcjet prototypes in the 1-kW class, the plume effect on the antenna performance is small. As the electron density increases, the main beam and sidelobe level gradually degrade. The main beam also tends to squint away from the plume region.

Ling, Hao↗

Gamma-ray transfer and energy deposition in supernovae

Solutions to the energy-independent (gray) radiative transfer equations are compared to results of Monte Carlo simulations of the Ni-56 and Co-56 decay gamma-ray energy deposition in supernovae. The comparison shows that an effective, purely absorptive, gray opacity, kappa(sub gamma) approximately (0. 06 +/- 0.01)Y(sub e) sq cm/g, where Y is the total number of electrons per baryon, accurately describes the interaction of gamma-rays with the cool supernova gas and the local gamma-ray energy deposition within the gas. The nature of the gamma-ray interaction process (dominated by Compton scattering in the relativistic regime) creates a weak dependence of kappa(sub gamma) on the optical thickness of the (spherically symmetric) supernova atmosphere: The maximum value of kappa(sub gamma) applies during optically thick conditions when individual gamma-rays undergo multiple scattering encounters and the lower bound is reached at the phase characterized by a total Thomson optical depth to the center of the atmosphere tau(sub e) approximately less than 1. Gamma-ray deposition for Type Ia supernova models to within 10% for the epoch from maximum light to t = 1200 days. Our results quantitatively confirm that the quick and efficient solution to the gray transfer problem provides an accurate representation of gamma-ray energy deposition for a broad range of supernova conditions.

Swartz, Douglas A.↗

Inhomogeneous conditions at open boundaries for wave propagation problems

Absorbing boundary conditions contain differential operators even for first-order systems. There is a fundamental difficulty with this, since the conditions are applied on the ingoing variables, and the approximations necessarily become weakly unstable. This difficulty is more pronounced for inhomogeneous boundary conditions which occur if there is a source outside the computational domain D, or if the initial data are nonzero outside D. In this paper this is further investigated and it is shown that reasonable estimates can still be obtained if the solution is smooth. However, it is demonstrated that the approximations are less robust. A cure for this is proposed and an implementation of high-order conditions for first-order systems is described.

Gustafsson, Bertil↗

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗

Numerical simulations of plasma double layers

The results of analytical studies of quasi-static electric fields along geomagnetic field lines are discussed. The calculations were targeted at the structure, generation mechanisms and stability parameters. The field consists of two oppositely charged layers, either weakly or strongly charged, with an electric field between. Existence conditions are defined for the double layer field and balancing requirements are explored. Details of the simulation techniques, i.e., particle in cell and Vlasov simulations, for studying the double layer are outlined, noting that both periodic and quasi-periodic simulations are used. Solutions to Poisson's equation for fixed and floating point boundary conditions are generated. Finally, attention is also given to oblique and two-dimensional magnetic double layers.

Goertz, C. K.↗

Current distribution in a three-dimensional IC analyzed by a perturbation method. Part 1: A simple steady state theory

The steady state current distribution in a three dimensional integrated circuit is presented. A device physics approach, based on a perturbation method rather than an equivalent lumped circuit approach, is used. The perturbation method allows the various currents to be expressed in terms of elementary solutions which are solutions to very simple boundary value problems. A Simple Steady State Theory is the subtitle because the most obvious limitation of the present version of the analysis is that all depletion region boundary surfaces are treated as equipotential surfaces. This may be an adequate approximation in some applications but it is an obvious weakness in the theory when applied to latched states. Examples that illustrate the use of these analytical methods are not given because they will be presented in detail in the future.

Edmonds, Larry D.↗

Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions

We examine the nonlinear stability of the Wisdom-Holman (WH) symplectic mapping applied to the integration of perturbed, highly eccentric (e-0.9) two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories for this class of problems, unless the step size chosen 1s small enough that PeriaPse is always resolved, in which case the method is generically stable. This 'radial orbit instability' persists even for weakly perturbed systems. Using the Stark problem as a fiducial test case, we investigate the dynamical origin of this instability and argue that the numerical chaos results from the overlap of step-size resonances; interestingly, for the Stark-problem many of these resonances appear to be absolutely stable. We similarly examine the robustness of several alternative integration methods: a time-regularized version of the WH mapping suggested by Mikkola; the potential-splitting (PS) method of Duncan, Levison, Lee; and two original methods incorporating approximations based on Stark motion instead of Keplerian motion. The two fixed point problem and a related, more general problem are used to conduct a comparative test of the various methods for several types of motion. Among the algorithms tested, the time-transformed WH mapping is clearly the most efficient and stable method of integrating eccentric, nearly Keplerian orbits in the absence of close encounters. For test particles subject to both high eccentricities and very close encounters, we find an enhanced version of the PS method-incorporating time regularization, force-center switching, and an improved kernel function-to be both economical and highly versatile. We conclude that Stark-based methods are of marginal utility in N-body type integrations. Additional implications for the symplectic integration of N-body systems are discussed.

Rauch, Kevin P.↗

Atmospheric propagation at larger lateral distances from the flight track

Sonic booms received on the ground tend to be restricted to a region of finite lateral extent below the flight track. This occurs because of refraction and because the effective speed of sound, even with winds taken into account, decreases with altitude in the lower atmosphere. Not all rays proceeding initially downwards from the flight track within an allowable range of initial directions will reach the ground. The restricted region which can be reached by rays impacting the ground is known as the primary carpet. However, weak rumbles are heard in the nominal shadow zone beyond the edge of this carpet. A full wave theory is necessary for explaining waveforms in that region, and the present paper gives a matched asymptotic expansion technique for a suitable approximate full wave theory that involves a relatively small number of parameters. The outer solution is derived from the structure of the system of rays that impact near the corridor edge; the inner solution involves a solution of the parabolic equation and results in the special functions encountered in the diffraction of sound over the tops of hills.

Pierce, Allan D.↗

A Hybrid RANS-LES Perspective for the High Lift Common Research Model Using LAVA

An assessment of a Hybrid RANS/LES (HRLES) approach for𝐶𝐿,max prediction is presented for the NASA High-Lift Common Research Model (CRM-HL). Both the free air and the wind tunnel configuration of the CRM-HL are investigated and the results are compared to the QinetiQ wind tunnel experiments and to two other numerical approaches: Reynolds Averaged Navier-Stokes (RANS) and Wall-Modeled Large Eddy Simulations (WMLES). For the free-air configuration, HRLES was shown to address some of the known shortcomings in RANS methods and prevent inboard and outboard flow separation particularly in the region of𝐶𝐿,max and post-stall. To achieve these improvements over RANS, LES-appropriate grids and numerical discretizations are required. HRLES predicts a weak pitch break at the highest angle-of-attack due to onset of wing-root corner flow separation whereas the free-air corrected experiment values indicate an occurrence of a much stronger pitch break. The improvements of HRLES over a URANS approach has been objectively shown by computing a set of solutions with the same grid, same numerics and time-step size and comparing the solutions. It was also found that when applying HRLES to a RANS best practice grid and numerics that the HRLES method significantly under performed RANS. For the in tunnel configuration, HRLES showed good agreement with the loads, surface pressure and oil-flow photographs obtained in the experiment. HRLES was able to improve upon the RANS simulations, which showed a sharp loss of lift at the two highest angles-of-attack due to large scale inboard and outboard separation on the wing, by correctly predicting the corner flow separation and showing remarkably close agreement in the flow topologies with the experiment.

TTT↗

Probing Pre-Galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature approximately greater than10(exp 4) K.We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) with values of about a few cm(exp −3). In more massive halos, corresponding to the first galaxies, the density may be larger, n approximately greater than100 cm(exp −3). The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, respectively, out to redshift z approximately greater than 20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages (C II, O I, Si II and Fe II). We show that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Metal↗

Short-time Lyapunov exponent analysis and the transition to chaos in Taylor-Couette flow

The physical mechanism driving the weakly chaotic Taylor-Couette flow is investigated using the short-time Liapunov exponent analysis. In this procedure, the transition from quasi-periodicity to chaos is studied using direct numerical 3D simulations of axially periodic Taylor-Couette flow, and a partial Liapunov exponent spectrum for the flow is computed by simultaneously advancing the full solution and a set of perturbations. It is shown that the short-time Liapunov exponent analysis yields more information on the exponents and dimension than that obtained from the common Liapunov exponent calculations. Results show that the chaotic state studied here is caused by a Kelvin-Helmholtz-type instability of the outflow boundary jet of Taylor vortices.

Vastano, John A.↗

The Lunar GNSS Receiver Experiment (LuGRE)

The Lunar GNSS Receiver Experiment (LuGRE) is a joint NASA-Italian Space Agency (ASI) payload on the Firefly Blue Ghost Mission 1 (BGM1) with the goal to demonstrate GNSS-based positioning, navigation, and timing at the Moon. LuGRE was chosen by the NASA Commercial Lunar Payload Services (CLPS) program as one of ten payloads on its “19D” task order for delivery to the lunar surface in 2023. The LuGRE payload consists of a weak-signal GNSS receiver, a high-gain L-band patch antenna, a low-noise amplifier, and an RF filter. The receiver will track GPS L1 C/A and L5, and Galileo E1 and E5a signals and will return pseudorange, carrier phase, and Doppler measurements to the ground. It will also calculate least-squares point solutions and Kalman-filter based navigation solutions onboard. In addition, the receiver features the capability to record raw I/Q baseband samples for downlink and ground processing. LuGRE will build on the legacy of prior missions in the Space Service Volume (SSV) including the initial experiments by AMSAT-OSCAR 40 and others, the GOES-R series of geostationary weather satellites, and the NASA Magnetospheric Multiscale (MMS) mission currently operating on GPS-based navigation at nearly 50% of lunar distance. Further, LuGRE will be one of the very first demonstrations of GNSS signal reception and navigation in the lunar environment and on the lunar surface, paving the way for operational use by future lunar missions such as Orion, Gateway, robotic and human landers, and surface rovers. Ultimately, all LuGRE science data will be released to a public data archive for the benefit of the GNSS and space communities. This paper provides a detailed overview of the LuGRE payload, including its design, concept of operations, and its predicted ability to meet its core science objectives. The baseline science investigations and priorities are outlined. Simulated performance results are shown based on the latest calibrated models including signal strength, signal availability, onboard navigation performance and convergence properties, and ground-based post-processed navigation performance.

LuGRE↗

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗

Computation of viscous blast wave flowfields

A method to determine unsteady solutions of the Navier-Stokes equations was developed and applied. The structural finite-volume, approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the interaction of blast-waves with stationary targets. The inviscid flux is evaluated using MacCormack's modified Steger-Warming flux or Roe flux difference splittings with total variation diminishing limiters, while the viscous flux is computed using central differences. The use of implicit boundary conditions in conjunction with a telescoping in time and space method permitted solutions to this strongly unsteady class of problems. Comparisons of numerical, analytical, and experimental results were made in two and three dimensions. These comparisons revealed accurate wave speed resolution with nonoscillatory discontinuity capturing. The purpose of this effort was to address the three-dimensional, viscous blast-wave problem. Test cases were undertaken to reveal these methods' weaknesses in three regimes: (1) viscous-dominated flow; (2) complex unsteady flow; and (3) three-dimensional flow. Comparisons of these computations to analytic and experimental results provided initial validation of the resultant code. Addition details on the numerical method and on the validation can be found in the appendix. Presently, the code is capable of single zone computations with selection of any permutation of solid wall or flow-through boundaries.

Atwood, Christopher A.↗

Polar Motion Constraints on Models of the Fortnightly Tide

Estimates of the near-fortnightly Mf ocean tide from Topex/Poseidon satellite altimetry and from numerical solutions to the shallow water equations agree reasonably well, at least in their basin-scale features. For example, both show that the Pacific Ocean tide lags the Atlantic tide by roughly 30 degrees. There are hints of finer scale agreements in the elevation fields, but noise levels are high. In contrast, estimates of Mf currents are only weakly constrained by the TP data, because high-wavenumber Rossby waves (with intense currents) are associated with relatively small perturbations in surface elevation. As a result, a wide range of Mf current fields are consistent with both the TP data and the hydrodynamic equations within a priori plausible misfit bounds. We find that a useful constraint on the Mf currents is provided by independent estimates of the Earth's polar motion. At the Mf period polar motion shows a weak signal (both prograde and retrograde) which must be almost entirely caused by the ocean tide. We have estimated this signal from the SPACE2000 time series, after applying a broad-band correction for atmospheric angular momentum. Although the polar motion estimates have relatively large uncertainties, they are sufficiently precise to fix optimum data weights in a global ocean inverse model of Mf. These weights control the tradeoff between fitting a prior hydrodynamic model of Mf and fitting the relatively noisy T/P measurements of Mf. The predicted polar motion from the final inverse model agrees remarkably well with the Mf polar motion observations. The preferred model is also consistent with noise levels suggested by island gauges, and it is marginally consistent with differences observed by subsetting the altimetry (to the small extent that this is possible). In turn, this new model of the Mf ocean tide allows the ocean component to be removed from Mf estimates of length of day, thus yielding estimates of complex Love numbers less contaminated by oceanic effects than has hitherto been possible.

Ray, Richard D.↗

A powerful local shear instability in weakly magnetized disks. III - Long-term evolution in a shearing sheet. IV - Nonaxisymmetric perturbations

The nonlinear evolution of the recently identified accretion disk magnetic shear instability is investigated through a series of numerical simulations. Finite-difference computations of the equations of compressible MHD are carried out on an axisymmetric shearing sheet system with periodic boundary conditions designed to approximate a local region within an accretion disk. Initial field configurations that include some net vertical component evolve into a nonlinear, exponentially growing solution with large poloidal velocities and magnetic fields with energies comparable to the thermal energy density. The stability of a purely azimuthal field configuration is examined, and it is found that nonaxisymmetric instability is present, but with a growth time measured in tens of orbital periods. In general, the most rapid growth occurs for very small radial and azimuthal wavenumbers, leading to coherent magnetic field structure in planes parallel to the disk. It is suggested that this instability is a key ingredient for the generation of magnetic fields in disks.

Hawley, John F.↗