Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solutes”

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 307 records · Page 17

Electromagnetic fields radiated from a lightning return stroke - Application of an exact solution to Maxwell's equations

A solution is presented for the electromagnetic fields radiated by an arbitrarily oriented current filament over a conducting ground plane in the case where the current propagates along the filament at the speed of light, and this solution is interpreted in terms of radiation from lightning return strokes. The solution is exact in the fullest sense; no mathematical approximations are made, and the governing differential equations and boundary conditions are satisfied. The solution has the additional attribute of being specified in closed form in terms of elementary functions. This solution is discussed from the point of view of deducing lightning current wave forms from measurements of the electromagnetic fields and understanding the effects of channel tortuosity on the radiated fields. In addition, it is compared with two approximate solutions, the traditional moment approximation and the Fraunhofer approximation, and a set of criteria describing their applicability are presented and interpreted.

D M Le Vine↗

Stability of PMR-polyimide monomer solutions

The stability of alcohol solutions of norbornenyl capped PMR-polyimide resins (both PMR-15 and PMR-II) has been monitored during storage at ambient and subambient (5 C and -18 C) temperatures. PMR-15 solutions exhibit nominally two weeks of useful life and PMR-II solutions exhibit nominally two days of useful life at ambient conditions. The limiting factor is precipitation of imide reaction products from the monomer solutions. Both solutions exhibit substantially longer useful lifetimes in subambient storage. PMR-15 shows no precipitation after several months storage at subambient temperatures. PMR-II solutions do exhibit precipitates after extended subambient storage (nominally three months at -18 C), however, the precipitates formed under these conditions can be redissolved. The chemical implications of these observations are discussed.

Lauver, R. W.↗

Solution sensitivity and accuracy study of NASTRAN for large dynamic problems involving structural damping

Large dynamic problems involving NASTRAN SOLUTION 8 (i.e., the steady state dynamic response option wherein all response quantities vary as e sub i omega t, where omega is the driving frequency and t is time) are considered. Using a submerged steel plate with a viscoelastic layer layer as the bench mark sample, the solution sensitivity and solution accuracy are checked. The solution sensitivity is examined by running the same finite element model on different computers, different versions of NASTRAN, and different precision levels. The solution accuracy is evaluated for these same runs by comparing the NASTRAN results with the exact solution of the same problem.

Kalinowski, A. J.↗

Influence of convection on free growth of dendrite crystals from solution

The free growth of dendrites in a uniformly supercooled solution was examined using cine photography with a Schlieren optical system. Crystals were grown in the bulk of the solution from a centrally located capillary tube, nucleated at the interface with a liquid nitrogen cooled wire. Crystals propagated along the tube, the slower growing orientations eliminated, and emerged at the tip, usually growing parallel to the tube direction. For both sodium sulfate decahydrate from its solution and ice from sodium chloride solution, growth rate and fineness of dendrites increased with supercooling. In sodium sulfate, upward convection of the less dense depleted solution occurs; downward convection was observed for the rejected, more concentrated sodium chloride solution. In both cases, there was a spatial and temporal delay in the release of the convective plume from the moving dendrite tip. The role of this convection on the growth characteristics and the production of secondary crystals is examined. A proposed low-g experiment to examine differences in growth rate, crystal texture, and secondary nucleation in a reduced convective regime where molecular diffusion is the dominant transfer process is discussed.

Hallett, J.↗

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗

Euler and Navier-Stokes solutions for supersonic shear flow past a circular cylinder

Euler and Navier-Stokes solutions of the supersonic shear flow past a circular cylinder are obtained. These solutions are used to study the basic flow structure around the cylinder. Both the inviscid and viscous calculations show the formation of a large recirculating flow region around the front stagnation point. The calculations further show that the overall size of the recirculating region is approximately the same for the Euler and Navier-Stokes solutions but the inside structure is quite different. The inviscid flow shows only one vortex whereas the viscous flow shows two vortices inside the recirculating flow region. The inner vortex in the Navier-Stokes solution is formed primarily due to the viscous effects near the body surface and its size depends upon the Reynolds number. It is found that with increasing Reynolds number, the inner vortex diminishes in size and the Navier-Stokes solution asymptotically approaches the Euler solution. These results indicate that the Euler equations may correctly predict certain high Reynolds number separation phenomenon in flows with natural inviscid vorticity source.

Kumar, A.↗

Stresses and deformations in cross-ply composite tubes subjected to a uniform temperature change: Elasticity and Approximate Solutions

The effects of a uniform temperature change on the stresses and deformations of composite tubes are investigated. The accuracy of an approximate solution based on the principle of complementary virtual work is determined. Interest centers on tube response away from the ends and so a planar elasticity approach is used. For the approximate solution a piecewise linear variation of stresses with the radial coordinate is assumed. The results from the approximate solution are compared with the elasticity solution. The stress predictions agree well, particularly peak interlaminar stresses. Surprisingly, the axial deformations also agree well. This, despite the fact that the deformations predicted by the approximate solution do not satisfy the interface displacement continuity conditions required by the elasticity solution. The study shows that the axial thermal expansion coefficient of tubes with a specific number of axial and circumferential layers depends on the stacking sequence. This is in contrast to classical lamination theory which predicts the expansion to be independent of the stacking arrangement. As expected, the sign and magnitude of the peak interlaminar stresses depends on stacking sequence.

Hyer, M. W.↗

On similarity solutions for turbulent and heated round jets

Commonly used empirical correlations for incompressible, heated round jets are shown to represent similarity solutions of the governing jet equations. These solutions give rise to self-similar eddy viscosities. Not all the similarity solutions are physically valid because some lead to zero eddy viscosities at the jet centerline. One physically valid solution is found to correlate best with round jet measurements and it gives a Gaussian error function description for the normalized mean velocity and temperature. Heat and momentum fluxes thus calculated are also in good agreement with measurements. Therefore, in addition to the classical similarity solution obtained by assuming constant eddy viscosity, another similarity solution to the jet equations is found where the eddy viscosity is self-similar.

So, R. M. C.↗

Numerical solution methods for viscoelastic orthotropic materials

Numerical solution methods for viscoelastic orthotropic materials, specifically fiber reinforced composite materials, are examined. The methods include classical lamination theory using time increments, direction solution of the Volterra Integral, Zienkiewicz's linear Prony series method, and a new method called Nonlinear Differential Equation Method (NDEM) which uses a nonlinear Prony series. The criteria used for comparison of the various methods include the stability of the solution technique, time step size stability, computer solution time length, and computer memory storage. The Volterra Integral allowed the implementation of higher order solution techniques but had difficulties solving singular and weakly singular compliance function. The Zienkiewicz solution technique, which requires the viscoelastic response to be modeled by a Prony series, works well for linear viscoelastic isotropic materials and small time steps. The new method, NDEM, uses a modified Prony series which allows nonlinear stress effects to be included and can be used with orthotropic nonlinear viscoelastic materials. The NDEM technique is shown to be accurate and stable for both linear and nonlinear conditions with minimal computer time.

Gramoll, K. C.↗

Raman intensity as a probe of concentration near a crystal growing in solution

The feasibility of using Raman spectral scattering signals for measurements of concentration profiles near a crystal interface during growth or dissolution is discussed. With KH2PO4 (KDP) as a test material, optical multichannel analyzer (OMA) detection of a solute Raman vibrational band provided direct quantification of solute concentration with band intensity. The intersection of incident laser and Raman collection optics provided 3-D selective point measurements of the solution concentration field. Unlike many other techniques, the Raman band intensity is not sensitive to the typical temperature variations. Precision calibration of Raman intensity versus KDP concentration with less than 1 pct standard deviation error levels was demonstrated. A fiber optic, which sampled incident laser intensity and coupled it to the OMA, provided a fully synchronized monitor of fluctuations in laser power to correlate with observed Raman signals. With 1 W of laser power at the sample, good data statistics required eight repeated data collections at approximately 2.5 min collection. The accumulated time represents the concentration measurement time at one spatial location. Photomicroscopy documented a 30 micrometer diameter by 200 micrometer of laser Raman scattering region in the solution near the crystal surface. The laser beam was able to approach up to 25 micrometer from the crystal surface. However, a crystal surface reflected intensity contribution was weakly detectable. Nucleated microcrystals were seen in the crystal-growing solution. These microcrystals convect right up to the crystal surface and indicate no quiet diffusion region under normal gravity conditions. Translation of the solution cell with respect to the optics caused systematic intensity errors.

Wilkinson, R. Allen↗

Crystallization of supercooled solutions

Crystallization of uniformly supercooled solutions (Na2SO4, NaCl, H2SO4, HNO3, HCl) was studied. It is shown how crystal growth velocity and habit depend on solution and concentration. The segregation coefficient for the solute in ice is measured by analysis of ice and solution, separated immediately after initial freezing, at different supercoolings. Subsequent solidification gives ion rejection at a varying rate depending on the geometry of the freezing, and may result in separation of hydrates, particularly when the initial concentration is high, as in haze (inactivated) droplets and low temperatures found in the Antarctic stratosphere. Electrical effects associated with rapid freezing are also investigated. Results suggest that more extensive measurements need to be made in solutions at different supercoolings, and that substantial electrical effects may be present for higher concentrations under these conditions. Damage to vegetation could occur under specific conditions as concentrated solutions (possibly H2SO4) are rejected in the freezing of rime or dew.

Harrison, K.↗

An interpretation and solution of ill-conditioned linear equations

Data insufficiency, poorly conditioned matrices and singularities in equations occur regularly in complex optimization, correlation, and interdisciplinary model studies. This work concerns itself with two methods of obtaining certain physically realistic solutions to ill-conditioned or singular algebraic systems of linear equations arising from such studies. Two efficient computational solution procedures that generally lead to locally unique solutions are presented when there is insufficient data to completely define the model, or a least-squares error formulation of this system results in an ill-conditioned system of equations. If it is assumed that a reasonable estimate of the uncertain data is available in both cases cited above, then we shall show how to obtain realistic solutions efficiently, in spite of the insufficiency of independent data. The proposed methods of solution are more efficient than singular-value decomposition for dealing with such systems, since they do not require solutions for all the non-zero eigenvalues of the coefficient matrix.

Ojalvo, I. U.↗

Long-period perturbations in Starlette orbit and tide solution

The long-period perturbations present in the orbit of the Starlette geodetic satellite is analyzed using a one-year continuous orbit from October 1976 through October 1977 and a three-year orbit during 1983 through 1985. A solution for 66 ocean tidal parameters for 14 low-degree tidal constituents is obtained in a simultaneous least squares estimation procedure with other geodynamic and orbit parameters. The yearly values of the second-degree S(a) tides, which have been observed to produce significant seasonal variations on the Starlette node, are included in the solution. Good agreement is found between the solutions, other satellite solutions, the oceanographic tide solutions, and other investigations using the computation of the earth's tidal response. The results indicate that the Starlette tide solutions produces significantly improved orbit fits for both Lageos and Starlette.

Cheng, M. K.↗

An adaptive grid method for computing time accurate solutions on structured grids

The solution method consists of three parts: a grid movement scheme; an unsteady Euler equation solver; and a temporal coupling routine that links the dynamic grid to the Euler solver. The grid movement scheme is an algebraic method containing grid controls that generate a smooth grid that resolves the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling is performed with a grid prediction correction procedure that is simple to implement and provides a grid that does not lag the solution in time. The adaptive solution method is tested by computing the unsteady inviscid solutions for a one dimensional shock tube and a two dimensional shock vortex iteraction.

Bockelie, Michael J.↗

The effect of viscosity on steady transonic flow with a nodal solution topology

The effect of viscosity on a steady, transonic flow for which the inviscid limit has a nodal solution topology near the critical point is investigated. For the accelerating case, viscous solutions tend to repel each other, so that a very delicate choice of initial conditions is required to prevent them from diverging. Only the two critical solutions extend to arbitrarily large distances into both the subsonic and supersonic flows. For the decelerating case, the solutions tend to attract, and so an entire two-parameter family of solutions now extends over large distances. The general effect of viscosity on the solution degeneracy of a nodal topology is thus to reduce or limit it for the accelerating case and to enhance it for the decelerating case. The astrophysical implications of these findings are addressed.

Owocki, Stanley P.↗

A two-dimensional Euler solution for an unbladed jet engine configuration

A two dimensional nonaxisymmetric Euler solution in a geometric representation of a jet engine configuration without blades is presented. The domain, including internal and external flow, is covered with a multiblock grid. To construct the grid, a domain decomposition technique is used to subdivide the domain and smooth grids are dimensioned and placed in each block. The grid contains 44 blocks which cover the external field, the inlet, bypass duct, core duct and nozzle of the nonaxisymmetric engine configuration. The geometry is symmetric about the meanline of the hub, but the grid is not since there is no symmetry condition applied to the grid between the two halves. With a symmetric grid at zero angle of attack, the measures of the solution would cancel exactly. With an asymmetric grid, the solution will not necessarily be symmetric and the lift coefficient will not necessarily be zero. Thus, grid asymmetry can be exploited to verify the resolution of the solution. The solution may be verified on the basis of five theoretical quantities: conservation of mass and energy, deviation of the lift coefficient from zero, deviation of the drag coefficient from zero, deviations from constant entropy, and deviations in the pressure distributions over the symmetric surfaces of the components. This technique is suitable for obtaining numerical solutions in complex geometries and provides a foundation for complete engine throughflow calculations.

Stewart, Mark E. M.↗

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↗

Iodine addition using triiodide solutions

The study develops: a triiodide solution for use in preparing ground service equipment (GSE) water for Shuttle support, an iodine dissolution method that is reliable and requires minimal time and effort to prepare, and an iodine dissolution agent with a minimal concentration of sodium salt. Sodium iodide and hydriodic acid were both found to dissolve iodine to attain the desired GSE iodine concentrations of 7.5 +/- 2.5 mg/L and 25 +/- 5 mg/L. The 1.75:1 and 2:1 sodium iodide solutions produced higher iodine recoveries than the 1.2:1 hydriodic acid solution. A two-hour preparation time is required for the three sodium iodide solutions. The 1.2:1 hydriodic acid solution can be prepared in less than 5 min. Two sodium iodide stock solutions (2.5:1 and 2:1) were found to dissolve iodine without undergoing precipitation.

Rutz, Jeffrey A.↗