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

Gravity flow of powder in a lunar environment. Part 2: Analysis of flow initiation

A small displacement-small strain finite element technique utilizing the constant strain triangle and incremental constitutive equations for elasticplastic (media nonhardening and obeying a Coulomb yield condition) was applied to the analysis of gravity flow initiation. This was done in a V-shaped hopper containing a powder under lunar environmental conditions. Three methods of loading were examined. Of the three, the method of computing the initial state of stress in a filled hopper prior to drawdown, by adding material to the hopper layer by layer, was the best. Results of the analysis of a typical hopper problem show that the initial state of stress, the elastic moduli, and the strength parameters have an important influence on material response subsequent to the opening of the hopper outlet.

Pariseau, W. G.↗

Application of a radiometric method to the measurement of temperature in self-absorbing gases

A radiometric method for the measurement of gas temperature in self-absorbing gases has been applied in the study of shock-tube-generated flows. This method involves making two absolute-intensity measurements at the same wavelength in a specified gas sample, but using two different path lengths. Experimental results are presented for reflected shock waves in air and in xenon at conditions corresponding to shock velocities from 4 to 10 km/sec and an initial driven tube pressure of 1 torr. These results indicate that, with the technique, temperature measurements with an accuracy of 5% can be carried out. The measurements in air suggest certain facility-related problems, and results obtained in the application of this measurement technique to the investigation of these problems and to the study of radiative cooling phenomena in both air and xenon are presented.

Nerem, R. M.↗

Stability of three-dimensional boundary layers

The computational modeling of the transition process characteristic of flows over swept wings is discussed. Specifically, the crossflow instability and crossflow/Tollmien-Schlichting wave interactions are analyzed through the numerical solution of the full three-dimensional Navier-Stokes equations including unsteadiness, curvature, and sweep. This approach is chosen because of the complexity of the problem and because it appears that regular stability theory is insufficient to explain the discrepancies between experiments and between theory and experiment. The leading edge region of a swept wing will be considered in a three-dimensional spatial simulation with random disturbances as the initial conditions.

Reed, Helen L.↗

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗

Temperature measurements behind reflected shock waves in air

A radiometric method for the measurement of gas temperature in self-absorbing gases has been applied in the study of shock tube generated flows. This method involves making two absolute intensity measurements at identical wavelengths, but for two different pathlengths in the same gas sample. Experimental results are presented for reflected shock waves in air at conditions corresponding to incident shock velocities from 7 to 10 km/s and an initial driven tube pressure of 1 torr. These results indicate that, with this technique, temperature measurements with an accuracy of + or - 5 percent can be carried out. The results also suggest certain facility related problems.

Bader, J. B.↗

A singular perturbation analysis of minimum time long range intercept

The application of singular perturbation theory for deriving algorithms suitable for on board real time computation of optimal aircraft trajectory control was investigated. Minimum time intercept in three dimensions was selected as a pilot problem formulation, and data for an early version F-4 aircraft was used to represent aircraft aerodynamic and propulsion characteristics. Results include the derivation of nonlinear, near optimal feedback control laws that are readily implementable in an aircraft flight computer, performance results for a variety of initial conditions (including optimal intercept at short ranges), and a point mass three dimensional simulation suitable for representing fighter aircraft under closed loop optimal control.

Calise, A.↗

A study of conditional instability of the second kind /CISK/

The theory of interaction between the cumulus clouds and the large-scale field proposed by Arakawa and Schubert (1974) is used to reexamine the CISK mechanism. After a qualitative discussion, a quantitative study is presented, using a three-level approach. The difference between the cases with conditional and unconditional heating formulations are discussed, and growth rate is treated as an eigenvalue problem. To confirm the results obtained and to extend the study to a case where the basic thermodynamic structure varies in the horizontal direction, the CISK mechanism is analyzed by an initial-value-problem approach. One of the findings is that the cumulus heating formulation of Arakawa and Schubert, through its scale dependence, is an improvement over previous CISK studies.

Chao, W. C.↗

A full-implicit-continuous-Eulerian (FICE) scheme for multidimensional transient magnetohydrodynamic (MHD) flows

A full implicit continuous Eulerian (FICE) scheme is developed for solving multidimensional transient MHD flow problems. The physical system under consideration is a general case of a transient MHD flow in which an initial steady state is subject to a finite amplitude disturbance. The governing equations are described, their finite difference formulation is presented and the FICE algorithm is given. The boundary conditions are treated by classifying them into physical and computational ones. The usefulness of the FICE algorithm is demonstrated using a physical example concerning the dynamical response of the static solar atmosphere due to a representative photospheric disturbance.

Wu, S. T.↗

The Impact of Land Initialization and Assimilation on Climate Predictability and Prediction

Analysis will be presented which explores the impact of land conditions on monthly to seasonal climate simulations in a variety of atmospheric general circulation models (AGCMs). In one set of experiments, the Geophysical Fluid Dynamics Laboratory (GDFL) AGCM is used to explore the nature of soil-moisture predictability and associated climate predictability as an initial value problem. For another set of experiments, the Center for Ocean Land Atmosphere (COLA) and the Goddard Earth Observing System 2 (GEOS-2) AGCMs are used to investigate the impact of realistic snow initialization and assimilation in retrospective climate forecasts for the northern hemisphere spring (March-June).

Schlosser, C. Adam↗

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗

Tycho: Realistically Simulating Exoplanets within Stellar Clusters. I. Improving the Monte Carlo Approach

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses >0.5M {sub ⊙} (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale N-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, Tycho, which employs a variety of different computational techniques, including multiple N-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1% of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which that perturb long-period planets lead to 38.3% of solar-system-like planetary systems becoming long-term unstable.

79 ASTRONOMY AND ASTROPHYSICS↗

Unconditional instability of inflow dependent boundary conditions in difference approximations to hyperbolic systems

The stability of finite difference approximations to initial boundary hyperbolic systems is studied. As is well known, a proper specification of boundary conditions for such systems is essential for their solutions to be well defined. A discrete analogue of the above is proved - if the numerical boundary conditions are consistent with an inflow part of the problem, they render the overall computation unstable. An example of the inviscid gasdynamics equations is considered.

Tadmor, E.↗

The unconditional instability of inflow-dependent boundary conditions in difference approximations to hyperbolic systems

The stability of finite difference approximations to initial boundary hyperbolic systems is studied. As is well known, a proper specification of boundary conditions for such systems is essential for their solutions to be well defined. A discrete analogue of the above is proved - if the numerical boundary conditions are consistent with an inflow part of the problem, they render the overall computation unstable. An example of the inviscid gasdynamics equations is considered. Previously announced in STAR as N81-33874

Tadmor, E.↗

The solution of linear systems of equations with a structural analysis code on the NAS CRAY-2

Two methods for solving linear systems of equations on the NAS Cray-2 are described. One is a direct method; the other is an iterative method. Both methods exploit the architecture of the Cray-2, particularly the vectorization, and are aimed at structural analysis applications. To demonstrate and evaluate the methods, they were installed in a finite element structural analysis code denoted the Computational Structural Mechanics (CSM) Testbed. A description of the techniques used to integrate the two solvers into the Testbed is given. Storage schemes, memory requirements, operation counts, and reformatting procedures are discussed. Finally, results from the new methods are compared with results from the initial Testbed sparse Choleski equation solver for three structural analysis problems. The new direct solvers described achieve the highest computational rates of the methods compared. The new iterative methods are not able to achieve as high computation rates as the vectorized direct solvers but are best for well conditioned problems which require fewer iterations to converge to the solution.

Poole, Eugene L.↗

Cascade Optimization Strategy Maximizes Thrust for High-Speed Civil Transport Propulsion System Concept

The design of a High-Speed Civil Transport (HSCT) air-breathing propulsion system for multimission, variable-cycle operations was successfully optimized through a soft coupling of the engine performance analyzer NASA Engine Performance Program (NEPP) to a multidisciplinary optimization tool COMETBOARDS that was developed at the NASA Lewis Research Center. The design optimization of this engine was cast as a nonlinear optimization problem, with engine thrust as the merit function and the bypass ratios, r-values of fans, fuel flow, and other factors as important active design variables. Constraints were specified on factors including the maximum speed of the compressors, the positive surge margins for the compressors with specified safety factors, the discharge temperature, the pressure ratios, and the mixer extreme Mach number. Solving the problem by using the most reliable optimization algorithm available in COMETBOARDS would provide feasible optimum results only for a portion of the aircraft flight regime because of the large number of mission points (defined by altitudes, Mach numbers, flow rates, and other factors), diverse constraint types, and overall poor conditioning of the design space. Only the cascade optimization strategy of COMETBOARDS, which was devised especially for difficult multidisciplinary applications, could successfully solve a number of engine design problems for their flight regimes. Furthermore, the cascade strategy converged to the same global optimum solution even when it was initiated from different design points. Multiple optimizers in a specified sequence, pseudorandom damping, and reduction of the design space distortion via a global scaling scheme are some of the key features of the cascade strategy. HSCT engine concept, optimized solution for HSCT engine concept. A COMETBOARDS solution for an HSCT engine (Mach-2.4 mixed-flow turbofan) along with its configuration is shown. The optimum thrust is normalized with respect to NEPP results. COMETBOARDS added value in the design optimization of the HSCT engine.

Source record↗

Integrated Sustainable Product Design With Warranty and End-of-Use Considerations

The concept of integrated sustainable product design has recently emerged, aiming to incorporate downstream life cycle performance into the initial product design to enhance sustainability. Various sustainable product design tools based on life cycle assessment or quality function deployment have been established while the impact of reliability on circular practices has received limited attention. Recognizing the critical role of product reliability in post-design performance, this paper develops a product design optimization model that considers the warranty performance and the effect of end-of-use options. The model takes into account the effect of uncertain operating conditions on product reliability. Two optimization goals including the minimization of expected unit life cycle cost and environmental impact are achieved by the model. To demonstrate the benefits of the integrated approach, the model is applied to an electric motor design problem. The results highlight that integrating end-of-use options in the early design phase leads to adjustments in component selection and reliability design. Moreover, the circular utilization of used products enables cost savings throughout the product’s life cycle and contributes to environmental impact reduction. Finally, the study analyzes the effects of operating conditions, warranty policies, and take-back prices for used products on design decisions, providing valuable insights for product designers.

Integrated design↗

How Much Reserve Fuel: Quantifying the Maximal Energy Cost of System Disturbances

Motivated by the design question of additional fuel needed to complete a task in an uncertain environment, this paper introduces metrics to quantify the maximal additional energy used by a control system in the presence of bounded disturbances, compared to a nominal, disturbance-free system. In particular, we consider the task of finite-time stabilization for a linear, time-invariant system. We compare the nominal energy required to achieve this task in the disturbance-free system to the worst-case energy over all feasible disturbances. Solving for the worst-case energy over all disturbances first leads to an optimal control problem with a least-squares solution, and then an infinite-dimensional optimization problem where we derive an upper bound on the solution. The comparison of energies is accomplished using additive and multiplicative metrics, for which we derive bounds. Simulation examples on an ADMIRE fighter jet model demonstrate the practicability of these metrics, and their variation with the distance of the initial condition from the origin and the task completion time.

koopman operator, resilience↗