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 127 records · Page 7

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.↗

The block adaptive multigrid method applied to the solution of the Euler equations

In the present study, a scheme capable of solving very fast and robust complex nonlinear systems of equations is presented. The Block Adaptive Multigrid (BAM) solution method offers multigrid acceleration and adaptive grid refinement based on the prediction of the solution error. The proposed solution method was used with an implicit upwind Euler solver for the solution of complex transonic flows around airfoils. Very fast results were obtained (18-fold acceleration of the solution) using one fourth of the volumes of a global grid with the same solution accuracy for two test cases.

Pantelelis, Nikos↗

Combined LAURA-UPS solution procedure for chemically-reacting flows

A new procedure seeks to combine the thin-layer Navier-Stokes solver LAURA with the parabolized Navier-Stokes solver UPS for the aerothermodynamic solution of chemically-reacting air flowfields. The interface protocol is presented and the method is applied to two slender, blunted shapes. Both axisymmetric and three dimensional solutions are included with surface pressure and heat transfer comparisons between the present method and previously published results. The case of Mach 25 flow over an axisymmetric six degree sphere-cone with a noncatalytic wall is considered to 100 nose radii. A stability bound on the marching step size was observed with this case and is attributed to chemistry effects resulting from the noncatalytic wall boundary condition. A second case with Mach 28 flow over a sphere-cone-cylinder-flare configuration is computed at both two and five degree angles of attack with a fully-catalytic wall. Surface pressures are seen to be within five percent with the present method compared to the baseline LAURA solution and heat transfers are within 10 percent. The effect of grid resolution is investigated and the nonequilibrium results are compared with a perfect gas solution, showing that while the surface pressure is relatively unchanged by the inclusion of reacting chemistry the nonequilibrium heating is 25 percent higher. The procedure demonstrates significant, order of magnitude reductions in solution time and required memory for the three dimensional case over an all thin-layer Navier-Stokes solution.

Wood, William A.↗

Solutions of the benchmark problems by the dispersion-relation-preserving scheme

The 7-point stencil Dispersion-Relation-Preserving scheme of Tam and Webb is used to solve all the six categories of the CAA benchmark problems. The purpose is to show that the scheme is capable of solving linear, as well as nonlinear aeroacoustics problems accurately. Nonlinearities, inevitably, lead to the generation of spurious short wave length numerical waves. Often, these spurious waves would overwhelm the entire numerical solution. In this work, the spurious waves are removed by the addition of artificial selective damping terms to the discretized equations. Category 3 problems are for testing radiation and outflow boundary conditions. In solving these problems, the radiation and outflow boundary conditions of Tam and Webb are used. These conditions are derived from the asymptotic solutions of the linearized Euler equations. Category 4 problems involved solid walls. Here, the wall boundary conditions for high-order schemes of Tam and Dong are employed. These conditions require the use of one ghost value per boundary point per physical boundary condition. In the second problem of this category, the governing equations, when written in cylindrical coordinates, are singular along the axis of the radial coordinate. The proper boundary conditions at the axis are derived by applying the limiting process of r approaches 0 to the governing equations. The Category 5 problem deals with the numerical noise issue. In the present approach, the time-independent mean flow solution is computed first. Once the residual drops to the machine noise level, the incident sound wave is turned on gradually. The solution is marched in time until a time-periodic state is reached. No exact solution is known for the Category 6 problem. Because of this, the problem is formulated in two totally different ways, first as a scattering problem then as a direct simulation problem. There is good agreement between the two numerical solutions. This offers confidence in the computed results. Both formulations are solved as initial value problems. As such, no Kutta condition is required at the trailing edge of the airfoil.

Tam, Christopher K. W.↗

Analytic solution of the thickness problem of a rectangular wing in steady subsonic flow

An analytic solution of the thickness problem of a rectangular wing with parabolic airfoil section in three-dimensional flow is presented. The free-air solution is obtained by integrating the equation of the axial perturbation velocity. The Prandtl-Glauert rule can be used to derive the subsonic solution. Parts of the free-air solution are verified by taking the limit of the axial peturbation velocity on the model surface as the wing span goes to infinity. Pressure coefficients are studied for a selected wing geometry in the flow field. The solution of the thickness problem of a rectangular wing in a rectangular wind tunnel is derived using the free-air solution and the method of images.

Ulbrich, N.↗

A simulation to study the feasibility of improving the temporal resolution of LAGEOS geodynamic solutions by using a sequential process noise filter

A key drawback to estimating geodetic and geodynamic parameters over time based on satellite laser ranging (SLR) observations is the inability to accurately model all the forces acting on the satellite. Errors associated with the observations and the measurement model can detract from the estimates as well. These 'model errors' corrupt the solutions obtained from the satellite orbit determination process. Dynamical models for satellite motion utilize known geophysical parameters to mathematically detail the forces acting on the satellite. However, these parameters, while estimated as constants, vary over time. These temporal variations must be accounted for in some fashion to maintain meaningful solutions. The primary goal of this study is to analyze the feasibility of using a sequential process noise filter for estimating geodynamic parameters over time from the Laser Geodynamics Satellite (LAGEOS) SLR data. This evaluation is achieved by first simulating a sequence of realistic LAGEOS laser ranging observations. These observations are generated using models with known temporal variations in several geodynamic parameters (along track drag and the J(sub 2), J(sub 3), J(sub 4), and J(sub 5) geopotential coefficients). A standard (non-stochastic) filter and a stochastic process noise filter are then utilized to estimate the model parameters from the simulated observations. The standard non-stochastic filter estimates these parameters as constants over consecutive fixed time intervals. Thus, the resulting solutions contain constant estimates of parameters that vary in time which limits the temporal resolution and accuracy of the solution. The stochastic process noise filter estimates these parameters as correlated process noise variables. As a result, the stochastic process noise filter has the potential to estimate the temporal variations more accurately since the constraint of estimating the parameters as constants is eliminated. A comparison of the temporal resolution of solutions obtained from standard sequential filtering methods and process noise sequential filtering methods shows that the accuracy is significantly improved using process noise. The results show that the positional accuracy of the orbit is improved as well. The temporal resolution of the resulting solutions are detailed, and conclusions drawn about the results. Benefits and drawbacks of using process noise filtering in this type of scenario are also identified.

Hartman, Brian Davis↗

The Effect of Solution Conditions on the Nucleation Kinetics of Tetragonal Lysozyme Crystals

An understanding of protein crystal nucleation rates and the effect of solution conditions upon them, is fundamental to the preparation of protein crystals of the desired size and shape for X-ray diffraction analysis. The ability to predict the effect of supersaturation, temperature, pH and precipitant concentration on the number and size of crystals formed is of great benefit in the pursuit of protein structure analysis. In this study we experimentally examine the effect of supersaturation, temperature, pH and sodium chloride concentration on the nucleation rate of tetragonal chicken egg white lysozyme crystals. In order to do this batch crystallization plates were prepared at given solution concentrations and incubated at three different temperatures over the period of one week. The number of crystals per well with their size and dimensions were recorded and correlated against solution conditions. Duplicate experiments indicate the reproducibility of the technique. Although it is well known that crystal numbers increase with increasing supersaturation, large changes in crystal number were also correlated against solution conditions of temperature, pH and salt concentration over the same supersaturation ranges. Analysis of these results enhance our understanding of the effect of solution conditions such as the dramatic effect that small changes in charge and ionic strength can have on the number of tetragonal lysozyme crystals that form and grow in solution.

Judge, Russell A.↗

A Low Order and a High Order Solution for a Converging-Diverging Nozzle Problem

The Computational Aeroacoustics Workshops on Benchmark problems are conducted in order to generate solutions with a variety of methods for problems that are chosen both to be representative of significant and relevant acoustic phenomena. It is generally recognized that CAA applications require the resolution in time and in space of solution details over a significantly broader range of scales than typical computational fluid dynamics applications. These requirements are forcing the development of new methodologies for CAA applications. This paper presents numerical results from a third and a seventh order algorithm for the propagation of an acoustic signal through a converging-diverging nozzle. These two algorithms are members of a new family of high accuracy methods that have the same order accuracy in both space and time, and are an extension of previous work for linearized Euler equations to fully nonlinear time dependent problems. The simulations are all with the fully nonlinear quasi-1D Euler equations for the total solution, which includes both the fluid dynamics and the acoustics. The acoustic solution is obtained from the time dependent nonlinear solution by subtracting the steady solution.

Goodrich, John W.↗

Tetragonal Lysozyme Nucleation and Crystal Growth: The Role of the Solution Phase

Lysozyme, and most particularly the tetragonal form of the protein, has become the default standard protein for use in macromolecule crystal nucleation and growth studies. There is a substantial body of experimental evidence, from this and other laboratories, that strongly suggests this proteins crystal nucleation and growth is by addition of associated species that are preformed by standard reversible concentration-driven self association processes in the bulk solution. The evidence includes high resolution AFM studies of the surface packing and of growth unit size at incorporation, fluorescence resonance energy transfer measurements of intermolecular distances in dilute solution, dialysis kinetics, and modeling of the growth rate data. We have developed a selfassociation model for the proteins crystal nucleation and growth. The model accounts for the obtained crystal symmetry, explains the observed surface structures, and shows the importance of the symmetry obtained by self-association in solution to the process as a whole. Further, it indicates that nucleation and crystal growth are not distinct mechanistically, but identical, with the primary difference being the probability that the particle will continue to grow or dissolve. This model also offers a possible mechanism for fluid flow effects on the growth process and how microgravity may affect it. While a single lysozyme molecule is relatively small (M.W. = 14,400), a structured octamer in the 4(sub 3) helix configuration (the proposed average sized growth unit) would have a M.W. = 115,000 and dimensions of 5.6 x 5.6 x 7.6 nm. Direct AFM measurements of growth unit incorporation indicate that units as wide as 11.2 nm and as long as 11.4 nm commonly attach to the crystal. These measurements were made at approximately saturation conditions, and they reflect the sizes of species that both added or desorbed from the crystal surface. The larger and less isotropic the associated species the more likely that it will be oriented to some degree in a flowing boundary layer, even at the low flow velocities measured about macromolecule crystals. Flow-driven effects resulting in misorientation upon addition to and incorporation into the crystal need only be a small fraction of a percentage to significantly affect the resulting crystal. One Earth, concentration gradient driven flow will maintain a high interfacial concentration, i.e., a high level (essentially that of the bulk solution) of solute association at the interface and higher growth rate. Higher growth rates mean an increased probability that misaligned growth units are trapped by subsequent growth layers before they can be desorbed and try again, or that the desorbing species will be smaller than the adsorbing species. In microgravity the extended diffusive boundary layer will lower the interfacial concentration. This results in a net dissociation of aggregated species that diffuse in from the bulk solution, i.e., smaller associated species, which are more likely able to make multiple attempts to correctly bind, yielding higher quality crystals.

Pusey, Marc L.↗

A New Global Mascon Solution Tuned for High-Latitude Ice Studies

A new global mascon solution has been developed with I-arc-degree spatial and IO-day temporal sampling. The global mas cons are estimated from the reduction of nearly 8 years of GRACE K-band range-rate data. Temporal and anisotropic spatial constraints have been applied for land, ocean and ice regions. The solution construction and tuning is focused towards the Greenland and Antarctic ice sheets (GIS and AIS) as well as the Gulf of Alaska mountain glaciers (GoA). Details of the solution development will be discussed, including the mascon parameter definitions, constraints, and the tuning of the constraint damping factor. Results will be presented, exploring the spatial and temporal variability of the ice sheets and GoA regions. A detailed error analysis will be discussed, including solution dependence on iteration, damping factor, forward modeling, and multitechnique comparisons. We also investigate the fundamental resolution of the solution and the spatial correlation of ice sheet inter-annual change. Finally, we discuss future improvements, including specific constraint application for the rest of the major land ice regions and improvements in solution regularization.

Luthcke, S. B.↗

An Iterated Global Mascon Solution with Focus on Land Ice Mass Evolution

Land ice mass evolution is determined from a new GRACE global mascon solution. The solution is estimated directly from the reduction of the inter-satellite K-band range rate observations taking into account the full noise covariance, and formally iterating the solution. The new solution increases signal recovery while reducing the GRACE KBRR observation residuals. The mascons are estimated with 10-day and 1-arc-degree equal area sampling, applying anisotropic constraints for enhanced temporal and spatial resolution of the recovered land ice signal. The details of the solution are presented including error and resolution analysis. An Ensemble Empirical Mode Decomposition (EEMD) adaptive filter is applied to the mascon solution time series to compute timing of balance seasons and annual mass balances. The details and causes of the spatial and temporal variability of the land ice regions studied are discussed.

Luthcke, S. B.↗