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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗
NASA NTRS2016

Tensor-Product Preconditioners for Higher-Order Space-Time Discontinuous Galerkin Methods

space-time discontinuous-Galerkin spectral-element discretization is presented for direct numerical simulation of the compressible Navier-Stokes equat ions. An efficient solution technique based on a matrix-free Newton-Krylov method is developed in order to overcome the stiffness associated with high solution order. The use of tensor-product basis functions is key to maintaining efficiency at high order. Efficient preconditioning methods are presented which can take advantage of the tensor-product formulation. A diagonalized Alternating-Direction-Implicit (ADI) scheme is extended to the space-time discontinuous Galerkin discretization. A new preconditioner for the compressible Euler/Navier-Stokes equations based on the fast-diagonalization method is also presented. Numerical results demonstrate the effectiveness of these preconditioners for the direct numerical simulation of subsonic turbulent flows.

Diosady, Laslo T.↗
NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗
NASA NTRS1990

Long-distance sound propagation over discontinuous impedances

A calculation method is presented for sound propagation over an impedance discontinuity in flat ground with a homogeneous, still atmosphere. The method is based on an approximate solution to a two dimensional boundary integral equation formulation of the problem, which expresses the wave field as the solution for homogeneous ground plus an integral over half of the boundary. Through recognizing this integral as a generalized Fourier integral, asymptotic methods are applied to evaluate the part of the integral most expensive to compute by numerical quadrature. Single frequency excess attenuation results for propagation from a point source above rigid ground to a receiver above absorbing ground are discussed. The results are applied, with air attenuation and A-weighting, to a notional jet engine noise source. Simple trends are noted.

Chandler-Wilde, Simon N.↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
NASA NTRS1979

Laminar and turbulent flows over a spherically blunted cone with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone with massive surface blowing. Time-dependent viscous shock layer equations are used to describe the flowfield. The boundary conditions on the body surface include a prescribed blowing rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the surface blowing smooths out the effect of the curvature discontinuity at the sphere-cone juncture point on the laminar flowfield and results in a negative pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer-like region on the flowfield and heating rates are discussed.

Kumar, A.↗
NASA NTRSDate not supplied

Exploration Helmet Permanent Anti-fog Study

For the current Extravehicular Mobility Unit (EMU) spacesuit, an astronaut applies an anti-fog solution to the interior of the helmet bubble before each EVA. However, the anti-fog solution has been reported to cause eye discomfort during at least seven EMU EVAs when the anti-fog solution contacted the crew member’s eyes. During STS-100, astronaut Chris Hadfield reported the eye irritation temporarily blinded him during his spacewalk. In addition, the wipe on anti-fog solution is a consumable that needs to be accounted for and a supply launched for missions. To solve this, the Exploration EMU (xEMU) pressure bubble investigated using a permanent anti-fog coating during Design Verification and Test (DVT) human in the loop (HITL) events. Chosen for previous use in Constellation Space Suit develop, HTAF-601, a water based permanent anti-fog solution, was tested. However, major issues have arisen with further HITL testing with the coating. Cleaning the helmet has been a challenge to avoid damaging the permanent anti-fog coating. NASA has completed a set of different methods for cleaning the anti-fog to try to document a preferred method. It was found isopropyl alcohol (IPA) cannot be used to clean the interior of the helmet because it will strip and delaminate the permanent anti-fog coating. Even with using a very gentle cleaning method of flushing with distilled or DI water and dabbing at facial oils, the permanent anti-fog starts to delaminate consistently after 50 hours manned pressurized time (MPT). Finally, the HTAF-601 coating is being discontinued by the vendor. Due to these issues, further investigation is being completed evaluating two solvent based anti-fog solutions: Exxene’s HCF-100 and FSI’s Visgard 106-94. Each coating will be evaluated on polycarbonate samples for application consistency, steam cycles (mimicking breath cycles), cleaning durability, haze, and light transmission. This study will determine if either coating is a viable option to pursue as a permanent anti-fog for spacesuit helmet applications.

spacesuit↗
NASA NTRS1977

Computation of the inviscid supersonic flow about cones at large angles of attack by a floating discontinuity approach

The technique of floating shock fitting is adapted to the computation of the inviscid flowfield about circular cones in a supersonic free stream at angles of attack that exceed the cone half-angle. The resulting equations are applicable over the complete range of free-stream Mach numbers, angles of attack and cone half-angles for which the bow shock is attached. A finite difference algorithm is used to obtain the solution by an unsteady relaxation approach. The bow shock, embedded cross-flow shock, and vortical singularity in the leeward symmetry plane are treated as floating discontinuities in a fixed computational mesh. Where possible, the flowfield is partitioned into windward, shoulder, and leeward regions with each region computed separately to achieve maximum computational efficiency. An alternative shock fitting technique which treats the bow shock as a computational boundary is developed and compared with the floating-fitting approach. Several surface boundary condition schemes are also analyzed.

Daywitt, J.↗
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 NTRS1995

Discontinuities formed by steepening the nonlinear Alfven waves

Montgomery (1959) first studied how Alfven waves steepen to form collisionless shocks by use of a two-fluid model. He solved the coupled Lorentz equations of motion for electrons and ions and Maxwell equations by a perturbation technique and showed that linearized polarized Alfven waves in the course of time steepen to form shocks. We have extended this work and have obtained exact simple wave solutions and show that in the limit of slow fluid time scales, we recover the earlier result of Montgomery. For an arbitrary fluid time scale, the solutions include elliptically polarized Alfven waves which also steepen and evolve into shocks. These nonlinear Alfven waves are important in the study of the dynamic solar wind discontinuities.

Lee, N. C.↗
NASA NTRS1995

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗
DOE OSTI2024

Vanishing adsorption limit of Riemann problem solutions for the polymer model

Here we examine the vanishing adsorption limit of solutions of Riemann problems for the Glimm–Isaacson model of chemical flooding of a petroleum reservoir. A contact discontinuity is deemed admissible if it is the limit of traveling waves or rarefaction waves for an augmented system that accounts for weak chemical adsorption onto the rock. We prove that this criterion justifies the admissibility criteria adopted previously by Keyfitz–Kranzer, Isaacson–Temple and de Souza–Marchesin, provided that the fractional flow function depends monotonically on chemical concentration. We also demonstrate that the adsorption criterion selects the undercompressive contact discontinuities required to solve the general Riemann problem in an example model with non-monotone dependence.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Similarity theory of the buoyantly interactive planetary boundary layer with entrainment

A similarity model is developed for the vertical profiles of turbulent flow variables in an entraining turbulent boundary layer of arbitrary buoyant stability. In the general formulation the vertical profiles, internal rotation of the velocity vector, discontinuities or jumps at a capping inversion and bulk aerodynamic coefficients of the boundary layer are given by solutions to a system of ordinary differential equations in the similarity variable. To close the system, a formulation for buoyantly interactive eddy diffusivity in the boundary layer is introduced which recovers Monin-Obukhov similarity near the surface and incorporates a hypothesis accounting for the observed variation of mixing length throughout the boundary layer. The model is tested in simplified versions which depend only on roughness, surface buoyancy, and Coriolis effects by comparison with planetary-boundary-layer wind- and temperature-profile observations, measurements of flat-plate boundary layers in a thermally stratified wind tunnel and observations of profiles of terms in the turbulent kinetic-energy budget of convective planetary boundary layers. On balance, the simplified model reproduced the trend of these various observations and experiments reasonably well, suggesting that the full similarity formulation be pursued further.

Hoffert, M. I.↗
NASA NTRS1978

A Mach line panel method for computing the linearized supersonic flow over planar wings

A method is described for solving the linearized supersonic flow over planar wings using panels bounded by two families of Mach lines. Polynomial distributions of source and doublet strength lead to simple, closed form solutions for the aerodynamic influence coefficients, and a nearly triangular matrix yields rapid solutions for the singularity parameters. The source method was found to be accurate and stable both for analysis and design boundary conditions. Similar results were obtained with the doublet method for analysis boundary conditions on the portion of the wing downstream of the supersonic leading edge, but instabilities in the solution occurred for the region containing a portion of the subsonic leading edge. Research on the method was discontinued before this difficulty was resolved.

Ehlers, F. E.↗
NASA NTRS1986

Aspects and applications of patched grid calculations

Patched grid calculations within the framework of an implicit, flux-vector split upwind/relaxation algorithm for the Euler equations are presented. The effect of a metric-discontinuous interface on the convergence rate of the algorithm is discussed along with the spatial accuracy of the solution and the effect of curvature along an interface. Results are presented and discussed for the free-stream problem, shock reflection problem, supersonic inlet with a 5 degree ramp, aerodynamically choked inlet, and three-dimensional analytic forebody.

Walters, R. W.↗
NASA NTRS1988

Investigation of stress concentration factor - Stress intensity factor interaction for flaws in filleted rods

To predict the fatigue life of machine parts which contain flaws in regions of stress concentration, the engineer needs to treat a series of crack sizes and/or shapes quickly and economically. In this paper the stress intensity factor distributions for circular surface flaws in the roots of fillets in rods loaded by remote tension are presented and examined. The examination reveals insight into the interaction of cracks with other stress concentrating discontinuities, and based on this it is suggested that weight function methods used in conjunction with known crack solutions may provide needed, efficient stress intensity factor prediction methods for such three-dimensional geometries.

Springfield, C. W., Jr.↗
NASA NTRS1993

Contact-impact analysis of deployable space systems

This paper presents simulation results of a flexible, deployable space system with multiple contact-impact events, using three different modeling approaches. In the first two approaches, the equations of motion of the system are written in terms of generalized coordinates of each component, with modal coordinates used to account for flexibility in the components. Two different methods are used to model contact-impact events in the systems. In the momentum balance method, as contact-impact occurs, momentum equations of the systems are solved for the jump discontinuity in velocity of generalized coordinates. Integration of the system is restarted with the new initial conditions from solutions of momentum equations. In the logical function method, a virtual spring connecting the two contact points is assumed. The virtual spring is active only when impact occurs. No special treatment is needed in the integration algorithm using this method. In the third approach, instead of using modal coordinates, a finite element based formulation for the transient dynamics of mechanical system is used to model flexibility of the components. The logical function method is used to model contact-impact in the third approach.

Wu, Shih-Chin↗
NASA NTRS2001

SGBEM-FEM Alternating Method for Analyzing 3D Non-Planar Cracks and Their Growth in Structural Components

The highly accurate and efficient Symmetric Galerkin Boundary Element Method (SGBEM), a Finite Element Method (FEM)-based alternating method, is proposed for analyzing three-dimensional non-planar cracks and their growth. The cracks are modeled using the symmetric Galerkin boundary element method as a distribution of displacement discontinuities, simulating an infinite medium. The finite element method only analyzes the stress for the uncracked body. The solution for the cracked structural component is determined by an iteration procedure. This process alternates between an FEM solution for the uncracked body and the SGBEM solution for a crack in an infinite body. Numerical analysis, and the Java code used, evaluate stress intensity factors and model fatigue crack growth. Examples of non-planar cracks in infinite media and planar cracks in finite bodies, as well as growth under fatigue, show the accuracy of the method.

Nikishkov, G. P.↗