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 55 records · Page 3

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗

Improved vortex methods for three-dimensional flows

Robust numerical methods are developed for three-dimensional incompressible vortical flows, using Lagrangian vortex elements. A successful scheme must be able to handle regions of intense vortex stretching and vortex reconnection with reasonable accuracy (without diverging). Here, consideration is given to vortex particles, also commonly called vortons or vortex sticks. The following issues are discussed: (1) use of delta-function elements and weak solutions of the vorticity equation; (2) use of smoothed elements and the choice of the smoothing function; (3) representation of viscous effects and the redistribution of element strength; and (4) conservation laws (are they satisfied?). The various proposed schemes have been tested on flows involving a strong interaction between two vortex rings.

Winckelmans, G.↗

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

Second-order accurate nonoscillatory schemes for scalar conservation laws

Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.

Huynh, Hung T.↗

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

A numerical method for systems of conservation laws of mixed type admitting hyperbolic flux splitting

The present treatment of elliptic regions via hyperbolic flux-splitting and high order methods proposes a flux splitting in which the corresponding Jacobians have real and positive/negative eigenvalues. While resembling the flux splitting used in hyperbolic systems, the present generalization of such splitting to elliptic regions allows the handling of mixed-type systems in a unified and heuristically stable fashion. The van der Waals fluid-dynamics equation is used. Convergence with good resolution to weak solutions for various Riemann problems are observed.

Shu, Chi-Wang↗

Numerical Wave Propagation 211 Based on Wave Primitives

Compact higher order finite difference equations are applied to a sequence of problems in wave propagation and aeroacoustics. Systems of PDE's are reduced to a sequence of simple wave primitives using a local eigenvector decomposition. The wave primitives are first order PDE's in two independent variable and allow natural boundary conditions to be imposed for both single- and multidimensional problems. The method uses a "discrete dispersion relation" approach to obtain high order approximations to the wave primitives on a 3 spatial point / 2 time level computational molecule. The scheme is fourth order accurate for the class of system with constant coefficients, e.g., those that support exponential solutions. Weakly non-linear PDE's are solved in a similar manner using a variant of the "method of frozen coefficients." Experience with the new algorithm for linear, non-linear, and multi-dimensional test problems will be described.

Davis, Sanford S.↗

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗

The calculation of supersonic viscous flows using the parabolized Navier-Stokes equations

Solution of the parabolic Navier-Stokes (PNS) equations for supersonic flows is discussed, and compatibility of the PNS method with the triple-deck theory of Stewartson (1974) is demonstrated. Characteristic and stability analyses show that use of an appropriate filter on the pressure term in the x-momentum equation can suppress the departure solutions, giving the usual desired weak interaction solution. An Alternating Direction Explicit procedure, with minimal computer storage requirements compared to the full Navier-Stokes solvers, is proposed to calculate strongly interacting flows using a global iteration procedure for the PNS equations. The PNS equations are used to solve the hypersonic viscous interaction problem, and good agreement is found with experimental results.

Davis, R. T.↗

An interaction solution algorithm for viscous aerodynamic performance prediction

A weak-interaction solution algorithm is established for aerodynamic flow field prediction about an isolated airfoil. It requires numerical solution of differential equations governing potential flow, viscous and turbulent boundary layer flow, and the turbulent wake flow downstream of a trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. These in turn alter the viscous flow through imposed pressure gradients. Closure for turbulence is accomplished using a second order model. Numerical evaluations assess factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems as solved using a finite element algorithm.

Baker, A. J.↗

Upwind schemes and bifurcating solutions in real gas computations

The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.

Suresh, Ambady↗

Solution-Phase Processes of Macromolecular Crystallization

We have proposed, for the tetragonal form of chicken egg lysozyme, that solution phase assembly processes are needed to form the growth units for crystal nucleation and growth. The starting point for the self-association process is the monomeric protein, and the final crystallographic symmetry is defined by the initial dimerization interactions of the monomers and subsequent n-mers formed, which in turn are a function of the crystallization conditions. It has been suggested that multimeric proteins generally incorporate the underlying multimers symmetry into the final crystallographic symmetry. We posed the question of what happens to a protein that is known to grow as an n-mer when it is placed in solution conditions where it is monomeric. The trypsin-treated, or cut, form of the protein canavalin (CCAN) has been shown to nucleate and grow crystals as a trimer from neutral to slightly acidic solutions. Under these conditions the solution is composed almost wholly of trimers. The insoluble protein can be readily dissolved by weakly basic solution, which results in a solution that is monomeric. There are three possible outcomes to an attempt at crystallization of the protein under monomeric (high pH) conditions: 1) we will obtain the same crystals as under trimer conditions, but at different protein concentrations governed by the self association equilibria; 2) we will obtain crystals having a different symmetry, based upon a monomeric growth unit; 3) we will not obtain crystals. Obtaining the first result would be indicative that the solution-phase self-association process is critical to the crystal nucleation and growth process. The second result would be less clear, as it may also reflect a pH-dependent shift in the trimer-trimer molecular interactions. The third result, particularly for experiments in the transition pH's between trimeric and monomeric CCAN, would indicate that the monomer does not crystallize, and that solution phase self association is not part of the crystal nucleation and growth path. Results are presented for crystallization experiments of CCAN over the pH 6.8 to 9.6 range.

Pusey, Marc L.↗

The Crystallization of Canavalin as a Function of pH and NaCl Concentration

We posed the question of what happens to a protein that is known to grow as an n-mer when it is placed in solution conditions where it is monomeric. The trypsin-treated, or cut, form of the protein canavalin (CCAN) has been shown to nucleate and grow crystals as a trimer from neutral to slightly acidic solutions. Under these conditions the solution is composed almost wholly of trimers. The crystalline protein can be readily dissolved by weakly basic solution, which has been proposed to result in a solution that is monomeric. There are three possible outcomes to an attempt at crystallization of the protein under monomeric (high pH) conditions: 1) we will obtain the same crystals as under trimer conditions, but at different protein concentrations governed by the self association equilibria; 2) we will obtain crystals having a different symmetry, based upon a monomeric growth unit; 3) we will not obtain crystals. Obtaining the first result would be indicative that the solution-phase self-association process is critical to the crystal nucleation and growth process. The second result would be less clear, as it may also reflect a pH-dependent shift in the trimer-trimer molecular interactions. The third result, particularly for experiments in the transition pH's between trimeric and monomeric CCAN, would indicate that the monomer does not crystallize, and that solution phase self association is not part of the crystal nucleation and growth path. Results are presented for crystallization experiments of CCAN over the pH 6.4 to 9.6 range. Fluorescence anisotropy, light scattering, and gel filtration experiments show that the solutions are primarily trimers, with association to form larger species occurring as a function of protein concentration.

Forsythe, Elizabeth L.↗

An approximate calculation of the strong interaction on a transonic airfoil

This paper describes an interactive computational scheme that calculates the transonic flow field over an airfoil. The inviscid flow, the turbulent boundary layer in the weak interaction region, and the elliptic flow in the strong interaction region are computed separately, and the solutions are matched iteratively in a direct, second-order interaction scheme. The pressure distribution, obtained from a transonic, small-perturbation solution, is prescribed on the boundary layer and the Navier-Stokes equations and the ensuing displacement thicknesses are used to modify the inviscid flow. The best overall results are obtained by using a turbulence model that eliminates the eddy viscosity in the strong interaction region. Computed pressure distributions are compared with experimental data and with Navier-Stokes solutions. In weak interactions, the agreement with experimental data is very good, as it is also in strong interactions except for the pressures in the separated-flow region. There the agreement is qualitative only, and the pressure recovery is too high. The results in both cases are superior or equal to those of the Navier-Stokes solutions, and the present method is significantly faster.

Seginer, A.↗

Sound propagation through nonuniform ducts

Methods of determining the transmission and attenuation of sound propagating in nonuniform ducts with and without mean flows are discussed. The approaches reviewed include purely numerical techniques, quasi-one-dimensional approximations, solutions for slowly varying cross sections, solutions for weak wall undulations, approximation of the duct by a series of stepped uniform cross sections, variational methods and solutions for the mode envelopes.

Nayfeh, A. H.↗

A computational study of the topology of vortex breakdown

A fully three-dimensional numerical simulation of vortex breakdown using the unsteady, incompressible Navier-Stokes equations has been performed. Solutions to four distinct types of breakdown are identified and compared with experimental results. The computed solutions include weak helical, double helix, spiral, and bubble-type breakdowns. The topological structure of the various breakdowns as well as their interrelationship are studied. The data reveal that the asymmetric modes of breakdown may be subject to additional breakdowns as the vortex core evolves in the streamwise direction. The solutions also show that the freestream axial velocity distribution has a significant effect on the position and type of vortex breakdown.

Spall, Robert E.↗