Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 451 records · Page 25

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗

Typical motions in multiple systems

In very old times, people counted - one, two, many. The author wants to show that they were right. Consider the motions of isolated bodies: (1) N = 1 - simple motion; (2) N = 2 - Keplerian orbits; and (3) N = 3 - this is the difficult problem. In general, this problem can be studied only by computer simulations. The author studied this problem over many years (see, e.g., Agekian and Anosova, 1967; Anosova, 1986, 1989 a,b). The principal result is that two basic types of dynamics take place in triple systems. The first special type is the stable hierarchical systems with two almost Keplerian orbits. The second general type is the unstable triple systems with complicated motions of the bodies. By random choice of the initial conditions, by the Monte-Carlo method, the stable systems comprised about approx. 10% of the examined cases; the unstable systems comprised the other approx. 90% of cases under consideration. In N greater than 3, the studies of dynamics of such systems by computer simulations show that we have in general also the motions roughly as at the cases 1 - 3 with the relative negative or positive energies of the bodies. In the author's picture, the typical trajectories of the bodies in unstable triple systems of the general type of dynamics are seen. Such systems are disrupted always after close triple approaches of the bodies. These approaches play a role like the gravitational slingshot. Often, the velocities of escapers are very large. On the other hand, the movie also shows the dynamical processes of a formation, dynamical evolution and disruption of the temporary wide binaries in triples and a formation of final hard massive binaries in the final evolution of triples.

Anosova, Joanna P.↗

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup -alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

Heat Diffusion in Gases, Including Effects of Chemical Reaction

The diffusion of heat through gases is treated where the coefficients of thermal conductivity and diffusivity are functions of temperature. The diffusivity is taken proportional to the integral of thermal conductivity, where the gas is ideal, and is considered constant over the temperature interval in which a chemical reaction occurs. The heat diffusion equation is then solved numerically for a semi-infinite gas medium with constant initial and boundary conditions. These solutions are in a dimensionless form applicable to gases in general, and they are used, along with measured shock velocity and heat flux through a shock reflecting surface, to evaluate the integral of thermal conductivity for air up to 5000 degrees Kelvin. This integral has the properties of a heat flux potential and replaces temperature as the dependent variable for problems of heat diffusion in media with variable coefficients. Examples are given in which the heat flux at the stagnation region of blunt hypersonic bodies is expressed in terms of this potential.

Hansen, C. Frederick↗

Relativistic tidal interaction of a white dwarf with a massive black hole

We compute encounters of a realistic white dwarf model with a massive black hole in the regime where relativistic effects are important, using a three-dimensional, finite-difference, Eulerian, piecewise parabolic method (PPM) hydrodynamical code. Both disruptive and nondisruptive encounters are considered. We identify and discuss relativistic effects important for the problem: relativistic shift of the pericenter distance, time delay, relativistic precession, and the tensorial structure of the tidal forces. In the nondisruptive case, stripping of matter takes place. In the surface layers of the surviving core, complicated hydrodynamical phenomena are revealed. In both disruptive and nondispruptive encounters, material flows out in the form of two thin, S-shaped, supersonic jets. Our results provide realistic initial conditions for the subsequent investigation of the dynamics of the debris in the field of the black hole. We evaluate the critical conditions for complete disruption of the white dwarf, and compare our results with the corresponding results for nonrelativistic encounters.

Frolov, V. P.↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Designing an Optimal Ensemble Strategy for GMAO S2S Forecast System

The NASA Global Modeling and Assimilation Office (GMAO) Sub-seasonal to Seasonal (S2S) prediction system is being readied for a major upgrade. An important factor in successful extended range forecasting is the definition of the ensemble. Our overall strategy is to run a relatively large ensemble of about 40 members up to 3 months (focusing on the sub-seasonal forecast problem), after which we sub-sample the ensemble, and continue the forecast with about 10 members (up to 12 months). Here we present the results of our testing of various ways to generate the initial perturbations and the validation of a stratified sampling approach for choosing the members of the smaller ensemble. For the initialization of the ensemble we propose a combination of lagged and burst initial conditions. To generate perturbations for the burst ensemble members we used scaled differences of pairs of analysis states (chosen randomly from the corresponding season) separated by 1-10 days. We consider perturbing separately the atmosphere and the ocean, or both. By varying the separation times between the analysis states, we are able to produce perturbations that resemble well-known modes of variability. Focusing on the ENSO SST indices, we found that all types of perturbations are important for the ensemble spread with, however, considerable differences in the timing of the impacts on spread for the atmospheric and oceanic perturbations.Our initial (larger) ensemble size was determined so as to maximize the skill of predicting some of the leading modes of boreal winter atmospheric modes (namely the NAO, PNA and AO). Since it is not feasible for us to run with the larger ensemble beyond about 3 months, we employ a stratified sampling procedure that identifies the emerging directions of error growth to subset the ensemble. By comparing the results from the stratified ensemble with that of the randomly sampled ensemble of the same size, we find that the former provides substantially better estimates the mean of the original large ensemble.

Borovikov, Anna↗

Developing Deep Learning Models for System Remaining Useful Life Predictions: Application to Aircraft Engines

Prognostics and health management (PHM) is an important part of ensuring reliable operations of complex safety- critical systems. System-level remaining useful life (RUL) estimation is a much more complex problem than making estimations at the component level, and system-level RUL methodologies remain sparse in the literature. Model-based approaches have traditionally worked in the past for components such as capacitors, MOSFETs, batteries, or hard-drives (to name a few examples), but developing high fidelity dynamics models of cyber physical systems that can be used to study the effects of multiple degrading components in the system remains a challenging task. Some initial work on model-based System RUL predictions was demonstrated in Khorasgani, et al [1], but, to generalize the system-level prognostics problem, we have to resort to pure data driven and hybrid approaches. In this work, we propose an end-to-end data- driven framework for developing deep learning models to predict remaining useful life of cyber physical systems operating under unknown faulty conditions. The raw data is organized with a data schema that improves the model development process and down stream data analysis tasks. Due to the unknown faulty conditions, the raw sensor data is transformed into signals that expose the underlying degradation processes, which are then used for model development. Bayesian Optimization is used to tune the model parameters prior to training and validation. We show that this approach results in accurate predictions within 3 cycles to end of life (EOL). We demonstrate the effectiveness of our approach by applying it to the N-CMAPSS turbofan engine dataset recently released by NASA, which includes high fidelity degradation modeling, real world operating conditions, and a large set of fault operating modes.

Prognostics↗

Three-dimensional hybrid-stress finite element analysis of composite laminates with cracks and cutouts

A three-dimensional hybrid-stress finite element analysis of composite laminates containing cutouts and cracks is presented. Fully three-dimensional, hexahedral isoparametric elements of the hybrid-stress model are formulated on the basis of the Hellinger-Reissner variational principle. Traction-free edges, cutouts, and crack surfaces are modeled by imposition of exact traction boundary conditions along element surfaces. Special boundary and surface elements are constructed by introducing proper constraints on assumed stress functions. The Lagrangian multiplier technique is used to enforce ply-interface continuity conditions in hybrid bimaterial composite elements for modeling the interface region in a composite laminate. Two examples are given to illustrate the capability of the present method of approach: (1) the well-known delamination problem in an angle-ply laminate, and (2) the important problem of a composite laminate containing a circular hole. Results are presented in detail for each case. Implications of interlaminar and intralaminar crack initiation, growth and fracture in composites containing cracks and cutouts are discussed.

Wang, S. S.↗

First- and second-order sensitivity analysis of linear and nonlinear structures

This paper employs the principle of virtual work to derive sensitivity derivatives of structural response with respect to stiffness parameters using both direct and adjoint approaches. The computations required are based on additional load conditions characterized by imposed initial strains, body forces, or surface tractions. As such, they are equally applicable to numerical or analytical solution techniques. The relative efficiency of various approaches for calculating first and second derivatives is assessed. It is shown that for the evaluation of second derivatives the most efficient approach is one that makes use of both the first-order sensitivities and adjoint vectors. Two example problems are used for demonstrating the various approaches.

Haftka, R. T.↗

TEMPUS: First results

The electromagnetic levitation facility, TEMPUS, developed by Dornier is designed to operate under microgravity conditions. Compared to terrestrial levitation, microgravity offers the possibility to melt and undercool in ultra high vacuum, thereby, providing an ultra clean environment. The technical concept of the TEMPUS facility was tested on two KC 135 flights and in the Texus 22 mission. Preparative investigations concerning the coil system and the heating and positioning efficiencies were carried out in the TEMPUS laboratory version. Furthermore, temperature-time profiles were determined under various boundary conditions. As a consequence of processing liquid metals under UHV, correct temperature measurement arises as the most critical problem. Experiences with experiments in the TEMPUS laboratory module show that due to the evaporation losses of the sample, the transmission of the CaF2 shielding windows changes drastically during the processing time. The investigation of the effect of contamination on pyrometry and the development of alternative evaporation shielding methods were initiated. During the second KC 135 flight, it was possible to heat up and melt an FeNi sample under He atmosphere. Oscillations of the molten sample, which were excited by switching out the magnetic heating field, could be detected and afterwards analyzed. From the frequency of these oscillations the surface tension of the sample material could be derived. The measurement of the surface tension and viscosity of an undercooled metal is proposed for TEMPUS on IML-2. This document is presented in view graph form.

Neuhaus, P.↗

TPSAS-NF1676L-12159-DND

Flight tests have demonstrated that in many operational environments, insect residue buildup on the leading edge of an aircraft wing can cause the loss of laminar flow. The airborne insect threat occurs between ground level and 152.4 m with the highest population concentration occurring during light winds (2.6 – 5.1 m/s), high humidity, and temperatures between 21-29°C. Under these conditions, mitigation is needed during takeoff, initial climb, approach, and landing. When insect fouling reaches a critical height, laminar flow is disrupted leading to increased drag and hence additional fuel consumption for the duration of the flight. This critical height is highly dependent upon the airfoil and Reynold’s number. Several techniques (e.g. paper, scrapers, surfactants, flexible surfaces, etc.) investigated to date have the potential to mitigate this problem; however, implementation has proven to be impractical thus limiting their utility. For example, the Kreuger flap (illustrated above) has been shown to negate insect excrescence from the wing leading edge while also improving lift during takeoff. However its complexity and physical discontinuities may induce early boundary layer transition on the lower surface.

T M Lorenzi↗

Advanced ceramic material for high temperature turbine tip seals

Ceramic material systems are being considered for potential use as turbine blade tip gas path seals at temperatures up to 1370 1/4 C. Silicon carbide and silicon nitride structures were selected for study since an initial analysis of the problem gave these materials the greatest potential for development into a successful materials system. Segments of silicon nitride and silicon carbide materials over a range of densities, processed by various methods, a honeycomb structure of silicon nitride and ceramic blade tip inserts fabricated from both materials by hot pressing were tested singly and in combination. The evaluations included wear under simulated engine blade tip rub conditions, thermal stability, impact resistance, machinability, hot gas erosion and feasibility of fabrication into engine components. The silicon nitride honeycomb and low-density silicon carbide using a selected grain size distribution gave the most promising results as rub-tolerant shroud liners. Ceramic blade tip inserts made from hot-pressed silicon nitride gave excellent test results. Their behavior closely simulated metal tips. Wear was similar to that of metals but reduced by a factor of six.

Solomon, N. G.↗

Demand for large freighter aircraft as projected by the NASA cargo/logistics airlift systems studies

This paper examines the market conditions up through the year 2008 to provide a preliminary assessment of the potential for and the characteristics of an advanced, all-cargo transport aircraft. Any new freighter must compete with current wide-body aircraft and their derivatives. Aircraft larger than the wide-bodies may incur economic penalties and operational problems. A lower direct operating cost is not a sufficient criterion to base a decision for the initiation of a new aircraft development or to select aircraft characteristics. Other factors of equal importance that are reviewed in this paper include considerations of the system infrastructure, the economics of the airlines, and the aircraft manufacturer return on investment. The results of the market forecast and a computer simulation show that an advanced long range aircraft with a payload between 68 to 181 tonnes (75 to 200 tons) could generate a solid foothold beginning around 1994.

Whitehead, A. H., Jr.↗

Interactions of light and gravity on growth, orientation, and lignin biosynthesis in mung beans

Mung beans (Vigna radiata L.) seedlings grown on the third Space Transport Mission (STS-3) showed marked orientation problems (some of the stems elongated horizontally and many of the roots were growing upward) and had a lower lignin content than the ground based controls. This research was initiated to determine if the atypical growth characteristics of mung beans grown in microgravity could be simulated using horizontal clinostats. Most of the effort focused on the design, construction and testing of the clinostats. In order to closely approximate the growth conditions of the plants grown in the plant growth unit on STS-3, cylindrical lexan minichambers were constructed. Results showed that plants grown using these clinostats in the horizontal position exhibit similar growth characteristics to the plants grown on STS-3 (disorientation of both stems and roots), while the vertical stationary and vertical rotating controls exhibit normal growth.

Jahns, G. C.↗

Existence and non-uniqueness of similarity solutions of a boundary layer problem

A Blasius boundary value problem with inhomogeneous lower boundary conditions f(0) = 0 and f'(0) = - lambda with lambda strictly positive was considered. The Crocco variable formulation of this problem has a key term which changes sign in the interval of interest. It is shown that solutions of the boundary value problem do not exist for values of lambda larger than a positive critical value lambda. The existence of solutions is proven for 0 lambda lambda by considering an equivalent initial value problem. It is found however that for 0 lambda lambda, solutions of the boundary value problem are nonunique. Physically, this nonuniqueness is related to multiple values of the skin friction.

Hussaini, M. Y.↗