Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuous 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Resolution of the COBE Earth sensor anomaly

Since its launch on November 18, 1989, the Earth sensors on the Cosmic Background Explorer (COBE) have shown much greater noise than expected. The problem was traced to an error in Earth horizon acquisition-of-signal (AOS) times. Due to this error, the AOS timing correction was ignored, causing Earth sensor split-to-index (SI) angles to be incorrectly time-tagged to minor frame synchronization times. Resulting Earth sensor residuals, based on gyro-propagated fine attitude solutions, were as large as plus or minus 0.45 deg (much greater than plus or minus 0.10 deg from scanner specifications (Reference 1)). Also, discontinuities in single-frame coarse attitude pitch and roll angles (as large as 0.80 and 0.30 deg, respectively) were noted several times during each orbit. However, over the course of the mission, each Earth sensor was observed to independently and unexpectedly reset and then reactivate into a new configuration. Although the telemetered AOS timing corrections are still in error, a procedure has been developed to approximate and apply these corrections. This paper describes the approach, analysis, and results of approximating and applying AOS timing adjustments to correct Earth scanner data. Furthermore, due to the continuing degradation of COBE's gyroscopes, gyro-propagated fine attitude solutions may soon become unavailable, requiring an alternative method for attitude determination. By correcting Earth scanner AOS telemetry, as described in this paper, more accurate single-frame attitude solutions are obtained. All aforementioned pitch and roll discontinuities are removed. When proper AOS corrections are applied, the standard deviation of pitch residuals between coarse attitude and gyro-propagated fine attitude solutions decrease by a factor of 3. Also, the overall standard deviation of SI residuals from fine attitude solutions decrease by a factor of 4 (meeting sensor specifications) when AOS corrections are applied.

Sedler, J.↗
NASA NTRS1993

State-constrained booster trajectory solutions via finite elements and shooting

This paper presents an extension of a FEM formulation based on variational principles. A general formulation for handling internal boundary conditions and discontinuities in the state equations is presented, and the general formulation is modified for optimal control problems subject to state-variable inequality constraints. Solutions which only touch the state constraint and solutions which have a boundary arc of finite length are considered. Suitable shape and test functions are chosen for a FEM discretization. All element quadrature (equivalent to one-point Gaussian quadrature over each element) may be done in closed form. The final form of the algebraic equations is then derived. A simple state-constrained problem is solved. Then, for a practical application of the use of the FEM formulation, a launch vehicle subject to a dynamic pressure constraint (a first-order state inequality constraint) is solved. The results presented for the launch-vehicle trajectory have some interesting features, including a touch-point solution.

Bless, Robert R.↗
NASA NTRS1974

A finite-element method for lifting surfaces in steady incompressible subsonic flow

The problem of potential steady subsonic flow for lifting surfaces is considered. This problem requires the solution of an integral equation relating the value of the potential discontinuity on the lifting surface and its wake to the values of the normal derivative of the potential which are known from the boundary conditions. The lifting surface is divided into small (quadrilateral hyperboloidal) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential discontinuity and the normal derivative of the potential are assumed to be constant within each element and equal to their values at the centroids of the elements. This yields a set of linear algebraic equations. Numerical results are in good agreement with existing ones.

Morino, L.↗
NASA NTRS1984

Numerical calculations of complex Mach reflection

Numerical simulations of the interaction of a planar blast wave with a compression ramp are presented. The split coefficient matrix (SCM) method in conjunction with boundary shock and floating discontinuity-fitting procedures was employed to obtain the time-asymptotic solutions of the two-dimensional, unsteady Euler equations. The solutions were computed for the complex Mach reflection (CMR) regime of the shock diffraction problem in an attempt to explore the basic physical process governing the evolution of an incipient second Mach stem and the associated topological changes. Numerical results were obtained for shock diffraction over a 40 degree ramp with varying incident shock Mach numbers. The validity of the present approach has been substantiated by experimental observations and earlier numerical calculations.

Yamamoto, O.↗
NASA NTRS2010

Development of Metallic Sensory Alloys

Existing nondestructive evaluation (NDE) technologies are inherently limited by the physical response of the structural material being inspected and are therefore not generally effective at the identification of small discontinuities, making the detection of incipient damage extremely difficult. One innovative solution to this problem is to enhance or complement the NDE signature of structural materials to dramatically improve the ability of existing NDE tools to detect damage. To address this need, a multifunctional metallic material has been developed that can be used in structural applications. The material is processed to contain second phase sensory particles that significantly improve the NDE response, enhancing the ability of conventional NDE techniques to detect incipient damage both during and after flight. Ferromagnetic shape-memory alloys (FSMAs) are an ideal material for these sensory particles as they undergo a uniform and repeatable change in both magnetic properties and crystallographic structure (martensitic transformation) when subjected to strain and/or temperature changes which can be detected using conventional NDE techniques. In this study, the use of a ferromagnetic shape memory alloy (FSMA) as the sensory particles was investigated.

Wallace Terryl A.↗
NASA NTRS1987

Shocked relativistic magnetohydrodynamic flows with application to pulsar winds

The time-dependent behavior of a shocked spherically symmetric relativistic fluid with tangential magnetic field is investigated, considering the case where the boundaries of the shocked fluid move at constant velocity so that self-similar solutions exist. The behavior of the fluid in the ultrarelativistic regime is compared to that in the nonrelativistic regime; there is a smooth transition between these limits. If a magnetic field is present, the magnetic pressure becomes increasingly important with distance from the shock wave; the gas pressure vanishes at the contact discontinuity that bounds the flow. Analytic expressions are given which describe the flow. The solutions can be applied to the evolution of shocked relativistic pulsar winds, which are probably observed as Crab-like supernova remnants. A model for the Crab Nebula, based on the steady-state model of Kennel and Coroniti (1984), indicates that sigma = 0.0016, where sigma is twice the ratio of magnetic to particle energy in the wind as measured in the fluid frame. This is about half the value suggested by Kennel and Coroniti and is much smaller than the value that might be expected for a pulsar wind.

Emmering, Robert T.↗
DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multicentered black hole saddles for supersymmetric indices

The supersymmetric index in string theory can sometimes have a discontinuous integer-valued jump at co-dimension one surfaces in moduli space called walls of marginal stability. When the index counts black hole microstates, crossing such walls of marginal stability amounts to the appearance or disappearance of a large number of such states. While wall-crossing has been understood in string theory and through the disappearance of extremal Lorentzian supergravity solutions as the moduli are varied, there has been no understanding about how the discontinuous changes in the index occur at the level of the gravitational path integral. In this paper, we find the finite-temperature saddles in $4d$ flatspace supergravity in which fermionic fields are periodic when going around the thermal circle that correspond to the multi-center black hole contributions to the index. By analyzing these saddles, we can explain how wall-crossing occurs: as the scalar moduli in supergravity are varied at the asymptotic boundary, for a given split of the charges, the saddle point equations can no longer be solved and, consequently, the corresponding multi-center saddle no longer contributes to the index. While the values of the scalars and the jump in the index when a wall is crossed all agree with the prediction from previously found Lorentzian supergravity solutions, the saddles in the index exhibit a much richer moduli space, which we analyze in detail.

FOS: Physical sciences↗
NASA NTRS1992

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗
DOE OSTI2025

On nonlocal problems with Neumann boundary conditions: scaling and convergence for nonlocal operators and solutions

Formulations of Neumann-type boundary conditions for boundary value problems in the nonlocal framework are beset with difficulties, some related to the choice of a proper scaling. Here we identify a space-dependent scaling for a nonlocal Neumann operator, for which we prove linear in δ (δ being the radius for the support for the kernel) convergence of the Neumann operator and $\mathcal{O}$(δ 2 ) convergence of solutions to their classical counterparts. The pointwise-like convergence of the nonlocal normal operator is cast as a new type of two-scale operator-point convergence, which we call condensated convergence . The results hold for general integrable kernels, a setting which is favored in numerical simulations. We support this analysis with numerical convergence studies using a piecewise linear discontinuous Galerkin discretization and show an $\mathcal{O}$(δ 2 ) rate of convergence of solutions, also exhibiting an $\mathcal{O}$(h 2 ) convergence, where h is the mesh size.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Viscous stability theory for thermally stratified flow - Discontinuous jets and shear layers.

The method used by Drazin to investigate the stability of unbounded, viscous, homogeneous, parallel shear flow to small wavenumber disturbances is extended to study the effect of thermal stratification on the stability of unbounded jets and shear layers. By this method the stability characteristics of continuous profiles are inferred from the stability characteristics of discontinuous profiles. The characteristic value problem for discontinuous jet and shear layers is posed by the requirement that the solutions of the governing differential equation satisfy the matching conditions and boundedness conditions for layers that extend to infinity. The analysis leads to a characteristic determinant which is required to vanish for the characteristic values of the parameters: the Reynolds number, the wavenumber, and the wave speed. The stabilizing effect of the thermal stratification as parameterized by the Richardson number was found to be most stabilizing for small wavenumber (large-scale) disturbances.

Gage, K. S.↗
NASA NTRS1997

The effects of programmed administration of human parathyroid hormone fragment (1-34) on bone histomorphometry and serum chemistry in rats

PTH treatment can result in dramatic increases in cancellous bone volume in normal and osteopenic rats. However, this potentially beneficial response is only observed after pulsatile treatment; continuous infusion of PTH leads to hypercalcemia and bone abnormalities. The purpose of these studies was to determine the optimal duration of the PTH pulses. A preliminary study revealed that human PTH-(1-34) (hPTH) is cleared from circulation within 6 h after sc administration of an anabolic dose of the hormone (80 microg/kg). To establish the effects of gradually extending the duration of exposure to hPTH without increasing the daily dose, we programmed implanted Alzet osmotic pumps to deliver the 80 microg/kg x day dose of the hormone during pulses of 1, 2, and 6 h/day, or 40 microg/kg x day continuously. Discontinuous infusion was accomplished by alternate spacing of external tubing with hPTH solution and sesame oil. After 6 days of treatment, we evaluated serum chemistry and bone histomorphometry. As negative and positive controls, groups of rats received pumps that delivered vehicle only and 80 microg/kg x day hPTH by daily sc injection, respectively. Dynamic and static bone histomorphometry revealed that the daily sc injection and 1 h/day infusion dramatically increased osteoblast number and bone formation in the proximal tibial metaphysis, whereas longer infusion resulted in systemic side-effects, including up to a 10% loss in body weight, hypercalcemia, and histological changes in the proximal tibia resembling abnormalities observed in patients with chronic primary hyperparathyroidism, including peritrabecular marrow fibrosis and focal bone resorption. Infusion for as little as 2 h/day resulted in minor weight loss and changes in bone histology that were intermediate between sc and continuous administration. The results demonstrate that the therapeutic interval for hPTH exposure is brief, but that programmed administration of implanted hormone is a feasible alternative to daily injection as a route for administration of the hormone.

Non-NASA Center↗
NASA NTRS2012

Elastic Model Transitions Using Quadratic Inequality Constrained Least Squares

A technique is presented for initializing multiple discrete finite element model (FEM) mode sets for certain types of flight dynamics formulations that rely on superposition of orthogonal modes for modeling the elastic response. Such approaches are commonly used for modeling launch vehicle dynamics, and challenges arise due to the rapidly time-varying nature of the rigid-body and elastic characteristics. By way of an energy argument, a quadratic inequality constrained least squares (LSQI) algorithm is employed to e ect a smooth transition from one set of FEM eigenvectors to another with no requirement that the models be of similar dimension or that the eigenvectors be correlated in any particular way. The physically unrealistic and controversial method of eigenvector interpolation is completely avoided, and the discrete solution approximates that of the continuously varying system. The real-time computational burden is shown to be negligible due to convenient features of the solution method. Simulation results are presented, and applications to staging and other discontinuous mass changes are discussed

Orr, Jeb S.↗
NASA NTRS1993

Liquid propellant rocket engine combustion simulation with a time-accurate CFD method

Time-accurate computational fluid dynamics (CFD) algorithms are among the basic requirements as an engineering or research tool for realistic simulations of transient combustion phenomena, such as combustion instability, transient start-up, etc., inside the rocket engine combustion chamber. A time-accurate pressure based method is employed in the FDNS code for combustion model development. This is in connection with other program development activities such as spray combustion model development and efficient finite-rate chemistry solution method implementation. In the present study, a second-order time-accurate time-marching scheme is employed. For better spatial resolutions near discontinuities (e.g., shocks, contact discontinuities), a 3rd-order accurate TVD scheme for modeling the convection terms is implemented in the FDNS code. Necessary modification to the predictor/multi-corrector solution algorithm in order to maintain time-accurate wave propagation is also investigated. Benchmark 1-D and multidimensional test cases, which include the classical shock tube wave propagation problems, resonant pipe test case, unsteady flow development of a blast tube test case, and H2/O2 rocket engine chamber combustion start-up transient simulation, etc., are investigated to validate and demonstrate the accuracy and robustness of the present numerical scheme and solution algorithm.

Chen, Y. S.↗
NASA NTRS1989

Mixed formulation for frictionless contact problems

Simple mixed finite element models and a computational precedure are presented for the solution of frictionless contact problems. The analytical formulation is based on a form of Reissner's large rotation theory of the structure with the effects of transverse shear deformation included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the internal forces (stress resultants), the generalized displacements, and the Lagrange multipliers associated with the contact conditions. The element characteristic array are obtained by using a modified form of the two-field Hellinger-Reissner mixed variational principle. The internal forces and the Lagrange multipliers are allowed to be discontinuous at interelement boundaries. The Newton-Raphson iterative scheme is used for the solution of the nonlinear algebraic equations, and the determination of the contact area and the contact pressures.

Noor, Ahmed K.↗
NASA NTRS2000

An Application of the Quadrature-Free Discontinuous Galerkin Method

The process of generating a block-structured mesh with the smoothness required for high-accuracy schemes is still a time-consuming process often measured in weeks or months. Unstructured grids about complex geometries are more easily generated, and for this reason, methods using unstructured grids have gained favor for aerodynamic analyses. The discontinuous Galerkin (DG) method is a compact finite-element projection method that provides a practical framework for the development of a high-order method using unstructured grids. Higher-order accuracy is obtained by representing the solution as a high-degree polynomial whose time evolution is governed by a local Galerkin projection. The traditional implementation of the discontinuous Galerkin uses quadrature for the evaluation of the integral projections and is prohibitively expensive. Atkins and Shu introduced the quadrature-free formulation in which the integrals are evaluated a-priori and exactly for a similarity element. The approach has been demonstrated to possess the accuracy required for acoustics even in cases where the grid is not smooth. Other issues such as boundary conditions and the treatment of non-linear fluxes have also been studied in earlier work This paper describes the application of the quadrature-free discontinuous Galerkin method to a two-dimensional shear layer problem. First, a brief description of the method is given. Next, the problem is described and the solution is presented. Finally, the resources required to perform the calculations are given.

Lockard, David P.↗
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.↗