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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

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↗
NASA NTRS1990

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

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

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

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↗
DOE OSTI2020

CPAP Ventilators Needed for Rapid Response to COVID-19 by Modification of CPAP Equipment

Early on in the COVID-19 pandemic, potential ventilator shortages were a critical issue identified by national health care providers. Capacity modeling at the time suggested patient demand may exceed ventilator supply. Thus, the challenge became finding an urgent interim solution to meet health care needs. Our initial hypothesis was that CPAP technology could be modified to provide similar functionality to a ventilator, relieving demand and allowing physicians to decide which patients need high end machines, ultimately saving lives. In conjunction with medical experts and pulmonologists, we were able to identify three key thrusts associated with this research problem: (1) modification of CPAP technology to allow for 02 input that would be capable of providing ventilation; (2) development of an alarming function that would provide real-time audible alarms to alert medical personnel to critical conditions, which would be used inline with CPAP technology; and (3) a method of sterilizing expiratory air from such a system in order to protect medical personnel from biohazard, since CPAPs vent to the atmosphere. We were unable to realize results for thrust 1 (CPAP modification for 02); we identified potential safety issues associated with utilizing medical grade oxygen with a common CPAP device. In order to characterize and mitigate these issues, we would need to partner closely with a device manufacturer; such a partnership could not be achieved in the timeframe needed for this rapid response work. However, we determined that some medical grade BiPAP devices do not need this modification and that the significant progress on thrusts 2 and 3 would be sufficient to buy down risk of a massive ventilator shortage. Our team built a prototype alarm system that can be utilized with any assistive respiratory device to alert on all key conditions identified by medical personnel (high pressure, low pressure, apnea, loss of power, low battery). Finally, our team made significant progress in the rapid prototyping and demonstration of an inline UV air purifier device. The device is cost efficient and can be manufactured at scale with both commercially available and additively manufactured parts. Initial tests with SARS-CoV-2 analog bacteriophage MS2 show 99% efficacy at reducing bioburden. Following a successful demonstration of the prototype device with medical personnel, we were able to obtain follow-on (non-LDRD) funding to provide additional device characterization, validation, and production in order to respond to an immediate regional need.

42 ENGINEERING↗
NASA NTRS1994

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

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↗
NASA NTRS2004

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↗
NASA NTRS2020

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↗
NASA NTRSDate not supplied

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↗
DOE OSTI2022

SRF Cavity Tuners for 3.9 GHz Cryomodules for LCLS-II Project

Fermilab conducted testing of three 3.9 GHz cryomodules for the LCLS-II project that will operate in continuous wave mode. A fast/fine tuning component was added to the LCLS-II 3.9 GHz tuner design due to the cavity bandwidth of 130 Hz which consists of two encapsulated piezos. Several cavities faced problems with fast-tuner operations after cooldown to 2 K and tuning the cavities to 3.9 GHz in cryomodule 2. All the piezo actuators were in working conditions but the slow tuner ranges required to stretch some of the cavities to the operational 3.9 GHz frequency were too small to deliver the required preload on the piezos. This behavior can be attributed to several factors: setting the initial warm cavity frequency during production too high, pressure tests of the warm cryomodule could have changed cavity frequency; and the small bending and twisting of the cavity-tuner system during the cooldown and warmup of the cavities. A decision was made to inelastically retune the warm cavities to decrease the unrestrained frequency by 200-300 kHz, this was done via the slow tuner. The results for this retuning method of three 3.9GHz cryomodules will be discussed.

43 PARTICLE ACCELERATORS↗
NASA NTRS1985

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

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

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

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↗
DOE OSTI2022

Modeling Combustion Reaction ODEs with Neural Networks

The chemistry of combustion reactions is complex as it involves many species and reactions. In practice, such a system is often modeled computationally using an empirically derived chemical mechanism. Given the large range of reaction rates, the system is then time evolved using a stiff ODE solver. However, even when reduced chemical mechanisms are employed, the system can become computationally expensive for two-dimensional or three-dimensional systems. Such systems can also face problems with instability. As such, it is desirable to find a cheaper, stable alternative to solving the reaction system. Neural networks offer the potential to learn these reaction ODEs and time evolve a combustion reaction in a more cost-efficient manner than stiff ODE solvers. In this study, a variety of neural networks are trained on zero-dimensional Cantera simulations of methane combustion with varying initial conditions. Several predictive approaches as well as several neural network architectures (artificial neural network with dropout, ResNet, and Neural ODE) are compared in their ability to predict combustion trajectories. Promising models are then identified.

combustion kinetics↗
DOE OSTI2022

VERA neutronics high-fidelity benchmark for a modern PWR core design

Nearly 20 years since the origin of the infamous 'Kord Smith Challenge', and after 10 years and hundreds of millions of dollars invested by the US Department of Energy into the development of the Virtual Environment for Reactor Applications (VERA), the capability to accurately simulate a modern pressurized water reactor (PWR) fuel cycle at the fuel rod level with high-fidelity Monte Carlo (MC) stochastic transport methods still evades us. The deterministic methods in VERA make approximations that allow it to be successfully benchmarked against measured data from hundreds of PWR fuel cycles, but none of these benchmarks provided information at the local fuel rod level or allowed quantification of these approximations' impact for modern PWR reload designs. Using MPACT to establish the estimated isotopic distributions and thermal hydraulic conditions in a recent and challenging reactor design, a 2D quarter-core benchmark problem was established and solved consistently with both MPACT, the 51-energy group deterministic transport solver in VERA, and Shift, the continuous-energy MC transport solver in VERA. Impacts of energy groups, ray spacing, and scattering treatment were evaluated. Though MPACT has been previously shown to perform well compared to MC methods for initial Cycle 1 cores with only fresh fuel, this is the first benchmark of MPACT for a challenging modern reload core with depleted fuel and burnable absorbers. This work establishes additional credibility for the methods in MPACT and highlights the efficiency of these methods relative to those in use by MC-based tools. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗