Search NASA⌕ Search

SEARCH · Search NASA

Results for “derivative discontinuity”

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

Gap Excitations and Series Loads in Microstrip Lines: Equivalent Network Characterization with Application to THz Circuits

At submillimeter wavelengths typical gap discontinuities in microstrip, CPW lines or at antenna terminals, which might contain diodes or active elements, cannot be viewed as simple quasi statically evaluated lumped elements. Planar Schottky diodes at 2.5 THz, for example, have a footprint that is comparable to a wavelength. Thus, apart from modelling the diodes themselves, the connection with their exciting elements (antennas or microstrip) gives rise to parasitics. Full wave or strictly numeric approaches can be used to account for these parasitics but at the expense of generality of the solution and the CPU time of the calculation. In this paper an equivalent network is derived that accurately accounts for large gap discontinuities (with respect to a wavelength) without suffering from the limitations of available numeric techniques.

Neto, Andrea↗
NASA NTRS1977

Nonlinear steady incompressible lifting-surface analysis with wake roll-up

The problem of lifting surfaces in steady incompressible flow is considered in terms of an integral equation relating the potential discontinuity on wing and wake to the normal derivative of the potential (normal wash) on the lifting surface. The integral equation is approximated by a system of linear algebraic equations obtained by dividing the surfaces into small quadrilateral elements and by assuming the potential discontinuity and the normal wash to be constant within each element. The wake geometry is obtained through iteration by satisfying the condition that the velocity be tangent to the surface of the wake and that the potential discontinuity be constant along the streamlines. Numerical results are in good agreement with existing ones

Suciu, E. O.↗
NASA NTRS1986

Application of NCOREL to aircraft configurations

The NCOREL computer program and methodology (Grossman, 1979) modified for computation of flows over complex geometries is described. The newly introduced features include a more flexible grid-generation package, capable of generating grids for realistic wing-body cross sections, and a more general numerical method, developed for the capture of highly three-dimensional embedded, oblique shocks. In addition, an inlet capability is included, by means of which the cross section is allowed to change discontinuously with the addition of the inlet geometry. The new cross section is regridded and the potentials and their derivatives are interpolated for the marching to continue. Hence, a discontinuity in the grid or mappings is taken into account at the inlet station. The results of the flow-field computation for a realistic fighter configuration are discussed.

Siclari, M. J.↗
NASA NTRS2023

Analytical model of foreshock ion interaction with a discontinuity: A statistical study

When a solar wind discontinuity interacts with foreshock ions, foreshock transients such as hot flow anomalies and foreshock bubbles can form. These create significant dynamic pressure perturbations disturbing the bow shock, magnetopause, and magnetosphere-ionosphere system. However, presently these phenomena are not predictable. In the accompanying paper, we derived analytical equations of foreshock ion partial gyration around a discontinuity and the resultant current density. In this study, we utilize the derived current density strength to model the energy conversion from the foreshock ions, which drives the outward motion or expansion of the solar wind plasma away from the discontinuity. We show that the model expansion speeds match those from local hybrid simulations for varying foreshock ion parameters. Using MMS, we conduct a statistical study showing that the model expansion speeds are moderately correlated with the magnetic field strength variations and the dynamic pressure decreases around discontinuities with correlation coefficients larger than 0.5. We use conjunctions between ARTEMIS and MMS to show that the model expansion speeds are typically large for those already-formed foreshock transients. Our results show that our model can be reasonably successful in predicting significant dynamic pressure disturbances caused by foreshock ion-discontinuity interactions. We discuss ways to improve the model in the future.

Terry Z. Liu↗
NASA NTRS1991

Derivative interface conditions for multiblock grids

Methods are developed for computing numerical solutions along block boundaries, even when there is a discontinuity in the grid lines or slopes. The technique is based on matching derivatives and does not require overlapping and interpolation of solution values at block boundaries. The comparison of block boundary values is implicit and has proven to be stable for both implicit and explicit numerical algorithms. Examples are included for the numerical solution of the Euler equations for compressible flow on grids with both grid line discontinuities and discontinuous slopes at block boundaries.

Mastin, C. W.↗
NASA NTRS2023

Foreshock Ion Motion Across Discontinuities: Formation of Foreshock Transients

In the ion foreshock, hot flow anomalies (HFAs) and foreshock bubbles (FBs) are two types of foreshock transients that have the strongest fluctuations, which can disturb the magnetosphere-ionosphere system and increase shock acceleration efficiency. They form due to interaction between the foreshock ions and solar wind discontinuities: the direction of the foreshock ion-driven current and whether it decreases or increases the magnetic field strength behind the discontinuity determine whether the transient's formation can be promoted or suppressed. Thus, to predict the HFA and FB formation and forecast their space weather effects, it is necessary to predict the foreshock ion-driven current direction. In this study, we derive analytical equations of foreshock ion velocities within discontinuities to estimate foreshock ion-driven current direction, which provides a quantitative criterion of HFA and FB formation. To validate the criterion, we use Acceleration Reconnection Turbulence & Electrodynamics of Moon's Interaction with the Sun to observe pristine solar wind discontinuities and calculate discontinuity parameters. We use Magnetospheric Multiscale to observe the foreshock ion motion around the discontinuities and show that the data support our model. This study is another step toward a predictive model of HFA and FB formation so that we can forecast their space weather effects at Earth using solar wind observations at lunar orbit or L1.

Terry Z. Liu↗
NASA NTRS2013

Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces

Nonlinear entropy stability and a summation-by-parts framework are used to derive provably stable, polynomial-based spectral collocation methods of arbitrary order. The new methods are closely related to discontinuous Galerkin spectral collocation methods commonly known as DGFEM, but exhibit a more general entropy stability property. Although the new schemes are applicable to a broad class of linear and nonlinear conservation laws, emphasis herein is placed on the entropy stability of the compressible Navier-Stokes equations.

Carpenter, Mark H.↗
NASA NTRS1980

Magnetic field directional discontinuities. I - Minimum variance errors

The paper deals with a statistical analysis of the errors associated with a minimum variance analysis of directional discontinuities by use of an idealized model of these discontinuities and various simulations, and also by an examination of actual Mariner 10 interplanetary magnetic field data. An empirical expression is derived for the magnitude of the error in an estimated discontinuity normal component, relative to the total field across the directional discontinuity. The analysis was performed primarily to aid in differentiating between interplanetary tangential and rotational discontinuities observed by Mariner 10.

Lepping, R. P.↗
DOE OSTI2022

Quantitative relationship between structural orthorhombicity, shear modulus, and heat capacity anomaly of the nematic transition in iron-based superconductors

Electronic nematicity in iron pnictide materials has been extensively studied by various experimental techniques, yet its heat capacity anomaly at the phase transition has not been examined quantitatively. In this paper, we review the thermodynamic description of nematicity in Ba(Fe 1- x Cox) 2 As 2 using the Landau free energy, which defines the behavior of three thermodynamic quantities: the structural orthorhombicity that develops below the nematic transition, the softening shear modulus above the transition, and the discontinuous heat capacity at the transition. We derive a quantitative relationship between these three quantities, which is found to hold for a range of dopings. We report this result shows that the nematic transition is exceedingly well described by a mean-field model in the underdoped regime of the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1988

JPL pulsar timing observations. V - Macro and microjumps in the Vela pulsar 0833-45

In this paper, 14.5 yr of JPL timing data are analyzed to study the spindown of the Vela pulsar. Discontinuities in the spin frequency and its derivatives have bimodal amplitude distributions. Macrojumps are characterized by relative changes in spin rate of about 10 to the -6th with a derivative of 0.01, while microjumps have relative changes in spin rate of less than about 10 to the -9th with derivative of less than about 0.0001. Macrojumps display a (+,-) signature with a well-defined ratio of the amplitudes. Microjumps, however, show all possible signs of events and a larger range in signature. Six macrojumps and their associated transients are quantified using a 10-parameter model for each jump. Amplitude and decay time parameters vary by a factor of about five over the set of jumps. The results are discussed in terms of superfluid vortex models. The macrojumps are consistent with a model in which angular momentum is suddenly transferred to the crust from a more rapidly rotating superfluid core, followed by relaxation of the spin-down rate.

Cordes, J. M.↗
DOE OSTI2021

Sequential discontinuities of Feynman integrals and the monodromy group

We generalize the relation between discontinuities of scattering amplitudes and cut diagrams to cover sequential discontinuities (discontinuities of discontinuities) in arbitrary momentum channels. The new relations are derived using time-ordered perturbation theory, and hold at phase-space points where all cut momentum channels are simultaneously accessible. As part of this analysis, we explain how to compute sequential discontinuities as monodromies and explore the use of the monodromy group in characterizing the analytic properties of Feynman integrals. We carry out a number of cross-checks of our new formulas in polylogarithmic examples, in some cases to all loop orders.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Discontinuous Galerkin and Related Methods for ODE

Starting from the standard integral formulation, the DG method is derived here in differential form. The key ingredient is a polynomial called the correction function, which helps ‘correct’ the discontinuous solution by approximating the jump and yields a continuous one. Under the right Radau quadrature, this continuous solution is identical to the solutions by the right Radau collocation and the continuous Galerkin (CG) methods. Next, the correction function facilitates the construction of the associated implicit Runge-Kutta schemes (IRK-DG). Different quadratures for DG result in different IRK-DG methods: left Radau quadrature in Radau IA, right Radau quadrature in Radau IIA or right Radau collocation, and Gauss quadrature in a method called DG-Gauss. The construction of IRK-DG clarifies the meaning and facilitates the proofs of various 𝐵(𝑝), 𝐶(𝜂), and 𝐷(𝜁) conditions for accuracy. The two consequences of these conditions are that all 𝑠-stage IRK-DG methods are accurate to order 2𝑠− 1, and the IRK-DG methods of Radau type are unique. Numerical examples showing the behavior of the DG solutions are provided. In all, the correction function plays a key role and helps establish the relations among the DG, IRK DG, collocation, and CG methods.

Numerical Methods for Ordinary Differential Equati↗
NASA NTRS1974

Stresses in dome-shaped shells of revolution with discontinuities at the apex

Asymptotic solutions of Novozhilov's equations of shells of revolution are derived for axisymmetric and first harmonic loadings. The solutions obtained are valid throughout the shallow and nonshallow regions. Stresses in dome shaped shells of revolution with a discontinuity in the form of a circular hole, a circular rigid insert, or a nozzle at the apex have been investigated. Numerical results are obtained for spheres, ellipsoids, and paraboloids, containing a discontinuity under a internal pressure and moment. Curves depicting stress distributions are given. The influence of three types of discontinuity on the stresses of the shells is also given.

Chen, C. H.↗
NASA NTRS1972

Influence of discontinuity stresses on main propellant tankage of a space shuttle orbiter.

Derivation of a general closed-form solution for the calculation of local stresses in cylindrical axial load-carrying pressure vessels which are heavily stiffened in the axial direction and circumferentially restrained by closely spaced ring stiffeners. The ring stiffeners are treated as discrete stiffeners and are allowed to have a temperature variation throughout the radial extension of the stiffener and different material properties compared with those of the vessel.

Meyer, H.-R.↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
NASA NTRSDate not supplied

Shock Capturing via Limiting for High-Order Methods including Discontinuous Galerkin

High-order methods, such as discontinuous Galerkin (DG), spectral, and flux reconstruction (FR), are prone to generating unwanted oscillations near shocks and discontinuities. Conventional limiting techniques, while effective in suppressing oscillations near shocks, often compromise accuracy near extrema, where the solution is only first-order accurate. This paper introduces a novel limiting technique for these high-order schemes, aimed at effectively managing shocks while preserving accuracy. The key idea is to expand the standard monotonicity limits to provide “room” near smooth extrema, ensuring that limiting has no effect and thus preserving accuracy. Near a discontinuity, these expanded limits effectively reduce to the original monotonicity limits, suppressing oscillations. Additional motivation is drawn from a formula for the derivative of Radau polynomials, which depicts the behavior of oscillations resulting from discontinuities. This behavior leads to a simplification by applying the limits to the sum of magnitudes of all modes, linear and higher degree. Unlike typical approaches, which rely on successful detection to activate limiting, our limiter depends continuously on the data, there by avoiding potential issues if detection fails. To reduce computing time, efficient criteria for detecting smooth regions where limiting is unnecessary are presented. Combined with detection, the continuous dependence on the data is lost, but the method is more economical. A notable characteristic of the entire process is its simplicity in both concept and implementation. Numerical tests for advection and Euler equations are conducted to demonstrate the effectiveness of the proposed method.

numerical methods↗
DOE OSTI2022

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗