Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 631 records · Page 35

Thermophysical properties of lysozyme (protein) solutions

Thermophysical properties of protein solutions composed of the lysozyme crystals with a 0.1 M sodium acetate and 5 percent NaCl solution as the buffer (pH = 4.0) are determined. The properties being measured include specific heat, thermal conductivity, dynamic viscosity, and surface tension. The protein concentrations are varied. Thermal diffusivity is calculated using the measured results. The purpose of the research is to measure thermophysical properties of lysozyme solutions which would serve as the data bank for controlling and modeling the crystal growth process on earth as well as in space.

Liu, Jiaching↗

Finite element solution of optimal control problems with state-control inequality constraints

It is demonstrated that the weak Hamiltonian finite-element formulation is amenable to the solution of optimal control problems with inequality constraints which are functions of both state and control variables. Difficult problems can be treated on account of the ease with which algebraic equations can be generated before having to specify the problem. These equations yield very accurate solutions. Owing to the sparse structure of the resulting Jacobian, computer solutions can be obtained quickly when the sparsity is exploited.

Bless, Robert R.↗

Navier-Stokes solutions with surface catalysis for Martian atmospheric entry

In this study numerical solutions have been obtained for two-dimensional axisymmetric hypersonic nonequilibrium CO2 flow over a high angle blunt cone with appropriate surface boundary conditions to account for energy and mass conservation at the body surface. The flowfield is described by the Navier-Stokes equations and multicomponent conservation laws which account for both translational and internal vibrational nonequilibrium effects. Complete forebody solutions have been obtained for the peak heating point of the Mars entry trajectory specified in the proposed NASA MESUR (Mars Environmental Survey) project. In these solutions, radiative equilibrium wall temperature and surface heating distributions are determined over the MESUR aeroshell forebody for entry velocity equal to 7 km/sec with varying degrees of surface catalysis. The effects of gas kinetics, surface catalysis, transport properties, and vibrational relaxation times on the surface heating are examined. The results identify some important issues in the prediction of surface heating for flows in thermochemical nonequilibrium and show that the Navier-Stokes code used herein is effective for thermal protection system design and materials selection.

Chen, Y.-K.↗

Monte Carlo exploration of Mikheyev-Smirnov-Wolfenstein solutions to the solar neutrino problem

The paper explores the impact of astrophysical uncertainties on the Mikheyev-Smirnov-Wolfenstein (MSW) solution by calculating the allowed MSW solutions for 1000 different solar models with a Monte Carlo selection of solar model input parameters, assuming a full three-family MSW mixing. Applications are made to the chlorine, gallium, Kamiokande, and Borexino experiments. The initial GALLEX result limits the mixing parameters to the upper diagonal and the vertical regions of the MSW triangle. The expected event rates in the Borexino experiment are also calculated, assuming the MSW solutions implied by GALLEX.

Shi, X.↗

Unit-Dose Bags For Formulating Intravenous Solutions

Smaller unit-dose flowthrough bags devised for use with large-volume parenteral (LVP) bags in preparing sterile intravenous solutions. Premeasured amount of solute stored in such unit-dose bag flushed by predetermined amount of water into LVP bag. Relatively small number of LVP bags used in conjunction with smaller unit-dose bags to formulate large number of LVP intravenous solutions in nonsterile environment.

Finley, Mike↗

Combined LAURA-UPS hypersonic solution procedure

A combined solution procedure for hypersonic flowfields around blunted slender bodies was implemented using a thin-layer Navier-Stokes code (LAURA) in the nose region and a parabolized Navier-Stokes code (UPS) on the after body region. Perfect gas, equilibrium air, and non-equilibrium air solutions to sharp cones and a sharp wedge were obtained using UPS alone as a preliminary step. Surface heating rates are presented for two slender bodies with blunted noses, having used LAURA to provide a starting solution to UPS downstream of the sonic line. These are an 8 deg sphere-cone in Mach 5, perfect gas, laminar flow at 0 and 4 deg angles of attack and the Reentry F body at Mach 20, 80,000 ft equilibrium gas conditions for 0 and 0.14 deg angles of attack. The results indicate that this procedure is a timely and accurate method for obtaining aerothermodynamic predictions on slender hypersonic vehicles.

Wood, William A.↗

A hybrid multigrid technique for computing steady-state solutions to supersonic flows

Recently, Li and Sanders have introduced a class of finite difference schemes to approximate generally discontinuous solutions to hyperbolic systems of conservation laws. These equations have the form together with relevant boundary conditions. When modelling hypersonic spacecraft reentry, the differential equations above are frequently given by the compressible Euler equations coupled with a nonequilibrium chemistry model. For these applications, steady state solutions are often sought. Many tens (to hundreds) of super computer hours can be devoted to a single three space dimensional simulation. The primary difficulty is the inability to rapidly and reliably capture the steady state. In these notes, we demonstrate that a particular variant from the schemes presented can be combined with a particular multigrid approach to capture steady state solutions to the compressible Euler equations in one space dimension. We show that the rate of convergence to steady state coming from this multigrid implementation is vastly superior to the traditional approach of artificial time relaxation. Moreover, we demonstrate virtual grid independence. That is, the rate of convergence does not depend on the degree of spatial grid refinement.

Sanders, Richard↗

Determination of monomer concentrations in crystallizing lysozyme solutions

We have developed a non-optical technique for the study of aggregation in lysozyme and other protein solutions. By monitoring the rate at which lysozyme traverses a semipermeable membrane it was possible to quantitate the degree of aggregation in supersaturated solutions. Using this technique, we have measured the concentration of monomers and larger aggregates in under- and oversaturated lysozyme solutions, and in the presence of crystals, at pH 4.0 and 3 percent NaCl (0.1M NaAc). Comparison of these concentration profiles with (110) face growth rate data supports the theory that tetragonal lysozyme crystals grow by addition of preformed aggregates and not by monomer addition. The data suggest that a considerable population of aggregates larger than dimers are present at lysozyme concentrations above 22 mg/ml. Determination of dimer concentrations, and equilibrium constants for subsequent aggregation levels, are currently underway.

Wilson, L. J.↗

Numerical solutions for heat flow in adhesive lap joints

The present formulation for the modeling of heat transfer in thin, adhesively bonded lap joints precludes difficulties associated with large aspect ratio grids required by standard FEM formulations. This quasi-static formulation also reduces the problem dimensionality (by one), thereby minimizing computational requirements. The solutions obtained are found to be in good agreement with both analytical solutions and solutions from standard FEM programs. The approach is noted to yield a more accurate representation of heat-flux changes between layers due to a disbond.

Howell, P. A.↗

An approximate solution for interlaminar stresses in composite laminates

An efficient approximate solution for interlaminar stresses in finite width, symmetric and unsymmetric laminated composites subjected to axial and/or bending loads is presented. The solution is based upon statically admissible stress fields which take into consideration local property mismatch effects and global equilibrium requirements. Unknown constants in the assumed stress states are determined through minimization of the laminate complementary energy. Typical results are presented for through-thickness and interlaminar stress distributions for angle-ply and cross-ply laminates subjected to axial loading. It is shown that the present formulation represents an improved, efficient approximate solution for interlaminar stresses.

Rose, Cheryl A.↗

A solution scheme for the Euler equations based on a multi-dimensional wave model

A scheme for the solution of scalar advection on an unstructured mesh has been developed, tested, and extended to the Euler equations. The scheme preserves a linear function exactly, and yields nearly monotone results. The flux function associated with the Euler scheme is based on a discrete 'wave model' for the system of equations. The wave model decomposes the solution gradient at a location into shear waves, entropy waves and acoustic waves and calculates the speeds, strengths and directions associated with the waves. The approach differs from typical flux-difference splitting schemes in that the waves are not assumed to propagate normal to the faces of the control volumes; directions of propagation of the waves are instead computed from solution-gradient information. Results are shown for three test cases, and two different wave models. The results are compared to those from other approaches, including MUSCL and Galerkin least squares schemes.

Powell, Kenneth G.↗

The structure of the gamma-domain for H-infinity Riccati solutions

The authors investigated some properties of the solutions to H-infinity Riccati equations which play a key role in the state approach to H-infinity optimization problem. The structure of the gamma-domain for the Riccati solutions has been described. There are two kinds of lower bounds for the solutions to exist and to be positive semidefinite respectively. A quadratically convergent algorithm for computing the lower bound is discussed.

Li, X. P.↗

Application of fast Fourier transforms to the direct solution of a class of two-dimensional separable elliptic equations on the sphere

An efficient, direct, second-order solver for the discrete solution of a class of two-dimensional separable elliptic equations on the sphere (which generally arise in implicit and semi-implicit atmospheric models) is presented. The method involves a Fourier transformation in longitude and a direct solution of the resulting coupled second-order finite-difference equations in latitude. The solver is made efficient by vectorizing over longitudinal wave-number and by using a vectorized fast Fourier transform routine. It is evaluated using a prescribed solution method and compared with a multigrid solver and the standard direct solver from FISHPAK.

Moorthi, Shrinivas↗

Discrete ordinates solutions of nongray radiative transfer with diffusely reflecting walls

Nongray gas radiation in a plane parallel slab bounded by gray, diffusely reflecting walls is studied using the discrete ordinates method. The spectral equation of transfer is averaged over a narrow wavenumber interval preserving the spectral correlation effect. The governing equations are derived by considering the history of multiple reflections between two reflecting wails. A closure approximation is applied so that only a finite number of reflections have to be explicitly included. The closure solutions express the physics of the problem to a very high degree and show relatively little error. Numerical solutions are obtained by applying a statistical narrow-band model for gas properties and a discrete ordinates code. The net radiative wail heat fluxes and the radiative source distributions are obtained for different temperature profiles. A zeroth-degree formulation, where no wall reflection is handled explicitly, is sufficient to predict the radiative transfer accurately for most cases considered, when compared with increasingly accurate solutions based on explicitly tracing a larger number of wail reflections without any closure approximation applied.

Menart, J. A.↗

Nonlinear vortex trail dynamics. II - Analytic solutions

Spatially periodic large amplitude solutions of the von Karman model are obtained in the neighborhood of singularities. These singularities correspond to vortex clusters in the physical plane. The quasi-periodic and unbounded solutions found analytically confirm earlier numerical work and show qualitative agreement with experimental observations of large-scale phenomena of vortex trails. Separatrices or heteroclinic orbits were explicitly found for an integrable approximate equation, which indicate that the von Karman model itself supports chaotic solutions.

Lim, Chjan C.↗

Automatic partitioning of unstructured meshes for the parallel solution of problems in computational mechanics

Most of the recently proposed computational methods for solving partial differential equations on multiprocessor architectures stem from the 'divide and conquer' paradigm and involve some form of domain decomposition. For those methods which also require grids of points or patches of elements, it is often necessary to explicitly partition the underlying mesh, especially when working with local memory parallel processors. In this paper, a family of cost-effective algorithms for the automatic partitioning of arbitrary two- and three-dimensional finite element and finite difference meshes is presented and discussed in view of a domain decomposed solution procedure and parallel processing. The influence of the algorithmic aspects of a solution method (implicit/explicit computations), and the architectural specifics of a multiprocessor (SIMD/MIMD, startup/transmission time), on the design of a mesh partitioning algorithm are discussed. The impact of the partitioning strategy on load balancing, operation count, operator conditioning, rate of convergence and processor mapping is also addressed. Finally, the proposed mesh decomposition algorithms are demonstrated with realistic examples of finite element, finite volume, and finite difference meshes associated with the parallel solution of solid and fluid mechanics problems on the iPSC/2 and iPSC/860 multiprocessors.

Farhat, Charbel↗

Implicit upwind solution algorithms for three-dimensional unstructured meshes

The development of implicit upwind algorithms for the solution of the three-dimensional, time-dependent Euler equations on unstructured tetrahedral meshes is described. The implicit temporal discretization involves either a two-sweep Gauss-Seide relaxation procedure, a two-sweep Point-Jacobi relaxation procedure, or a single-sweep Point-Implicit procedure; the upwind spatial discretization is based on the flux-difference splitting of Roe. Detailed descriptions of the three implicit solution algorithms are given, and calculations for the Boeing 747 transport configuration are presented to demonstrate the algorithms. Advantages and disadvantages of the implicit algorithms are discussed. A steady-state solution for the 747 configuration, obtained at transonic flow conditions using a mesh of over 100,000 cells, required less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the general capability.

Batina, John T.↗

Light-curve solutions for S Cancri and TT Hydrae with rapid rotation

Physical model light- and velocity-curve solutions for S Cancri and TT Hydrae are obtained, and analyses with incorporation of asynchronous rotation are carried out. A photometric rotation rate for the primary star of TT Hya is determined, and excellent agreement with results from spectral line profiles is found. Both separate light- and velocity-curve solutions and simultaneous light-velocity solutions are listed. The photometric rotation for S Cnc from existing light curves is indeterminate, but is compatible with line profile measures. Evidence for third light from the light curves of S Cnc is found. An explanation for the apparent conflict between the rotational states and mass-transfer activities of the two binaries is suggested.

Van Hamme, W.↗