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 325 records · Page 18

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

Evaluation of Linear, Inviscid, Viscous, and Reduced-Order Modeling Aeroelastic Solutions of the AGARD 445.6 Wing Using Root Locus Analysis

Reduced-order modeling (ROM) methods are applied to the CFD-based aeroelastic analysis of the AGARD 445.6 wing in order to gain insight regarding well-known discrepancies between the aeroelastic analyses and the experimental results. The results presented include aeroelastic solutions using the inviscid CAP-TSD code and the FUN3D code (Euler and Navier-Stokes). Full CFD aeroelastic solutions and ROM aeroelastic solutions, computed at several Mach numbers, are presented in the form of root locus plots in order to better reveal the aeroelastic root migrations with increasing dynamic pressure. Important conclusions are drawn from these results including the ability of the linear CAP-TSD code to accurately predict the entire experimental flutter boundary (repeat of analyses performed in the 1980's), that the Euler solutions at supersonic conditions indicate that the third mode is always unstable, and that the FUN3D Navier-Stokes solutions stabilize the unstable third mode seen in the Euler solutions.

Silva, Walter A.↗

Recent Developments and Challenges Implementing New and Improved Stress Intensity Factor (K) Solutions in NASGRO for Damage Tolerance Analyses

Fatigue crack growth analysis software has been available to damage tolerance analysts for many years in either commercial products or via proprietary in-house codes. The NASGRO software has been publicly available since the mid-80s (known as NASA/FLAGRO up to 1999) and since 2000 has been sustained and further developed by a collaborative effort between Southwest Research Institute® (SwRI®), the NASA Johnson Space Center (JSC), and the members of the NASGRO Industrial Consortium. Since the stress intensity factor (K) is the foundation of fracture mechanics and damage tolerance analysis of aircraft structures, a significant focus of development efforts in the past fifteen years has been geared towards enhancing legacy K solutions and developing new and efficient numerical K solutions that can handle the complicated stress gradients computed by today’s analysts using detailed finite element models of fatigue critical locations. This paper provides an overview of K solutions that have been recently implemented or improved for the analysis of geometries such as two unequal through cracks at a hole and two unequal corner cracks at a hole, as well as state-of-the-art weight function models capable of computing K in the presence of univariant and/or bivariant stress gradients and complicated residual stress distributions. Some historical background is provided to review how common K solutions have evolved over the years, including selective examples from the literature and from new research. Challenges and progress in rectifying discrepancies between older legacy solutions and newer models are reviewed as well as approaches and challenges for verification and validation of K solutions. Finally, a summary of current challenges and future research and development needs is presented. A key theme throughout the presentation of this paper will be how members of the aerospace industry have collaborated with software developers to develop a practical analysis tool that is used world-wide to support new design as well as the ongoing sustainment and airworthiness of commercial and military aircraft.

Cardinal, Joseph W.↗

An Analytical Solution for Yaw Maneuver Optimization on the International Space Station and Other Orbiting Space Vehicles

This paper suggests a new method for optimizing yaw maneuvers on the International Space Station (ISS). Yaw rotations are the most common large maneuvers on the ISS often used for docking and undocking operations, as well as for other activities. When maneuver optimization is used, large maneuvers, which were performed on thrusters, could be performed either using control moment gyroscopes (CMG), or with significantly reduced thruster firings. Maneuver optimization helps to save expensive propellant and reduce structural loads ‐ an important factor for the ISS service life. In addition, optimized maneuvers reduce contamination of the critical elements of the vehicle structure, such as solar arrays. This paper presents an analytical solution for optimizing yaw attitude maneuvers. Equations describing pitch and roll motion needed to counteract the major torques during a yaw maneuver are obtained. A yaw rate profile is proposed. Also the paper describes the physical basis of the suggested optimization approach. In the obtained optimized case, the torques are significantly reduced. This torque reduction was compared to the existing optimization method which utilizes the computational solution. It was shown that the attitude profiles and the torque reduction have a good match for these two methods of optimization. The simulations using the ISS flight software showed similar propellant consumption for both methods. The analytical solution proposed in this paper has major benefits with respect to computational approach. In contrast to the current computational solution, which only can be calculated on the ground, the analytical solution does not require extensive computational resources, and can be implemented in the onboard software, thus, making the maneuver execution automatic. The automatic maneuver significantly simplifies the operations and, if necessary, allows to perform a maneuver without communication with the ground. It also reduces the probability of command errors. The suggested analytical solution provides a new method of maneuver optimization which is less complicated, automatic and more universal. A maneuver optimization approach, presented in this paper, can be used not only for the ISS, but for other orbiting space vehicles.

Dobrinskaya, Tatiana↗

Dynamic Beam Solutions for Real-Time Simulation and Control Development of Flexible Rockets

In this study, flexible rockets are structurally represented by linear beams. Both direct and indirect solutions of beam dynamic equations are sought to facilitate real-time simulation and control development for flexible rockets. The direct solution is completed by numerically integrate the beam structural dynamic equation using an explicit Newmark-based scheme, which allows for stable and fast transient solutions to the dynamics of flexile rockets. Furthermore, in the real-time operation, the bending strain of the beam is measured by fiber optical sensors (FOS) at intermittent locations along the span, while both angular velocity and translational acceleration are measured at a single point by the inertial measurement unit (IMU). Another study in this paper is to find the analytical and numerical solutions of the beam dynamics based on the limited measurement data to facilitate the real-time control development. Numerical studies demonstrate the accuracy of these real-time solutions to the beam dynamics. Such analytical and numerical solutions, when integrated with data processing and control algorithms and mechanisms, have the potential to increase launch availability by processing flight data into the flexible launch vehicle's control system.

flexible body dyanmics fiber optic strain sensor p↗

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

Gnoffo, Peter A.↗

Method for Non-Invasive Determination of Chemical Properties of Aqueous Solutions

A method for non-invasively determining a chemical property of an aqueous solution is provided. The method provides the steps of providing a colored solute having a light absorbance spectrum and transmitting light through the colored solute at two different wavelengths. The method further provides the steps of measuring light absorbance of the colored solute at the two different transmitted light wavelengths, and comparing the light absorbance of the colored solute at the two different wavelengths to determine a chemical property of an aqueous solution.

Todd, Paul W.↗

An exact solution of an augmented Burgers equation and amplitude-dependent acoustic propagation speed

Nonlinear sound propagation in the atmosphere is usually modeled using an augmented Burgers equation accounting for a weak nonlinearity and atmospheric absorption. Because the absorption includes the molecular vibrational relaxation, such a Burgers equation is more complex than the regular Burgers equation that only accounts for the thermoviscous dissipation in the absorption. Although an exact solution of the regular Burgers equation has long been derived using the Cole- Hopf transform, an exact solution of the augmented Burgers equation has not been derived previously. Thus, this paper presents an exact solution of the augmented Burgers equation. This novel solution is shown to be equivalent to the solution using the Cole-Hopf transform when the absorption only involves thermoviscous dissipation. It can also be reduced to the known solution of an N-wave when the absorption is ignored. The augmented Burgers equation is an approximation valid for weak nonlinearity. However, this assumption may not be accurate for acoustic signals propagating from the lower atmosphere and which are subsequently refracted downward from the upper atmosphere (e.g., stratosphere and thermosphere) due to the decreasing air density with increasing altitude [Lonzaga, et al., Geophysical Journal International, 200(3), pp.1347-1361]. Consequently, the current paper also discusses the effects of a strong nonlinearity that lead to an amplitude-dependent increase in signal propagation speed. For an impulsive signal such as a sonic boom, these effects cause a dispersion of the signal similar to the observed dispersion of acoustic signals from supersonic Concorde as well as from large explosions.

Joel B Lonzaga↗

Hybrid Analytics Solution for Power Plant Operations

Expert Microsystems, Inc., MapEx Software Inc. and XMPLR Energy LLC have teamed to submit a proposal to the US Department of Energy’s Office of Fossil Energy to improve the efficiency, reliability and flexibility of existing coal-based power plants. The project, titled Hybrid Analytics Solution to Improve Coal Power Plant Operations (the “Hybrid Analytics” project), will develop, demonstrate and commercialize a novel approach that integrates two proven real-time monitoring techniques. The Hybrid Analytics approach integrates into a single, integrated solution: a) an established, advanced data-driven analytics solution that includes artificial intelligence, advanced pattern recognition and machine-learning techniques and b) a well-proven, first principle thermal heat balance model solution. The future role of coal plants will depend on the ability to cycle and follow loads to meet marginal power demands, accommodate renewables and support the grid. Essential to success in this new role is a plant’s ability to maintain high reliability, efficiency and flexibility. However, these changing operational demands put significant stresses on power plant components that directly impact the ability to deliver these requirements. In the past, plants have independently monitored performance using physics-based solutions (heat rate, boiler losses, etc.) and reliability using data-driven analytics (advanced pattern recognition). The Hybrid Analytics solution integrates both approaches to enhance fault detection, provide automated diagnostics and estimates of remaining time to act, and provide critical, enhanced and accurate real-time information to coal plant operators. As a result, coal plants can improve efficiency, reliability, and flexibility and can cost-effectively maintain their key role in the power delivery system.

01 COAL, LIGNITE, AND PEAT↗