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 595 records · Page 33

A Novel Method to Train Classification Models for Structure Detection in In Situ Spacecraft Data

We present a method for creating spacecraft-like data which can be used to train Machine Learning (ML) models to detect and classify structures in in situ spacecraft data. First, we use the Grad-Shafranov equation to numerically solve for several magnetohydrostatic equilibria which are variations on a known analytic equilibrium. These equilibria are then used as the initial conditions for Particle-In-Cell simulations in which the structures of interest are observed and labeled. We then take one-dimensional slices through the simulations to replicate what a spacecraft collecting data from the simulation would observe. This sliced data then can be used as training data for the initial training of ML models intended for use on spacecraft data. We demonstrate the method applied to the problem of detecting small-scale plasmoids in the magnetotail, which is important for understanding complex magnetotail reconnection dynamics. The simple 1D classifier we train is able to detect more than 70% of the plasmoid points in the data set but also produces a large number of false positives. Our further work on this example problem is detailed, and further potential uses of the method are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

Online Model-Free DER Dispatch Via Adaptive Voltage Sensitivity Estimation and Chance Constrained Programming

This paper proposes an online data-driven distributed energy resource management system (DERMS) for distribution system optimal DER dispatch as well as voltage regulation. Here, the key innovation is to leverage the Local Sensitivity Factor (LSF) for transforming the DER control into a computationally efficient linear programming (LP) problem. By taking real-time measurements, the estimation of LSF eliminates the need for an accurate distribution system model as well as full nodal load information, which is difficult to achieve in practice. A robust recursive least squares method is also developed to ensure the robust estimation of LSF, which is initialized using reasonable values from model-derived LSFs. This allows the system to adapt to changing operational conditions effectively. A scenario-based, chance-constrained framework is further employed to ensure voltage remains within acceptable limits in the presence of measurement and estimation uncertainties. Test results on a real-world, 759-node distribution network located in western Colorado, U.S., validate the effectiveness and robustness of the proposed control approach and demonstrate its superior performance as compared to alternative methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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

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

Quantum statistics of Raman scattering model with Stokes mode generation

The model describing three coupled quantum oscillators with decay of Rayleigh mode into the Stokes and vibration (phonon) modes is examined. Due to the Manley-Rowe relations the problem of exact eigenvalues and eigenstates is reduced to the calculation of new orthogonal polynomials defined both by the difference and differential equations. The quantum statistical properties are examined in the case when initially: the Stokes mode is in the vacuum state; the Rayleigh mode is in the number state; and the vibration mode is in the number of or squeezed states. The collapses and revivals are obtained for different initial conditions as well as the change in time the sub-Poisson distribution by the super-Poisson distribution and vice versa.

Tanatar, Bilal↗

Challenges in simulating ground interacting nuclear explosions

This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computation of unsteady transonic flows through rotating and stationary cascades. 3: Acoustic far-field analysis

A small perturbation type analysis has been developed for the acoustic far field in an infinite duct extending upstream and downstream of an axial turbomachinery stage. The analysis is designed to interface with a numerical solution of the near field of the blade rows and, thereby, to provide the necessary closure condition to complete the statement of infinite duct boundary conditions for the subject problem. The present analysis differs from conventional inlet duct analyses in that a simple harmonic time dependence was not assumed, since a transient signal is generated by the numerical near-field solution and periodicity is attained only asymptotically. A description of the computer code developed to carry out the necessary convolutions numerically is included, as well as the results of a sample application using an impulsively initiated harmonic signal.

Slutsky, S.↗

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

Structural analysis of a bolted joint concept for the space shuttle's solid rocket motor casing

The Space Shuttle Challenger accident is thought to have been caused by the failure of one of the tang-clevis joints joining together the casing segments of the Solid Rocket Motors (SRM). Excessive displacement between the tang and clevis, possibly unseating the O-ring seals, may have initiated the resulting accident. An effort was made at NASA Langley Research Center to design an alternative concept for mating the casing segments. A bolted flange joint concept was designed and analyzed to determine if the concept would effectively maintain a seal while minimizing joint weight and controlling stress levels. It is shown that under the loading conditions analyzed the seal area of the joint remains seated. The only potential stress problem is a stress concentration in the flange at the edge of the bolt hole, which is highly localized. While heavier than the existing joint, this concept does have some advantages making the bolted joint an attractive alternative.

Lindell, Michael C.↗

Terrestrial planet and asteroid formation in the presence of giant planets. I. Relative velocities of planetesimals subject to Jupiter and Saturn perturbations

We investigate the orbital evolution of 10(13)- to 10(25) -g planetesimals near 1 AU and in the asteroid belt (near 2.6 AU) prior to the stage of evolution when the mutual perturbations between the planetesimals become important. We include nebular gas drag and the effects of Jupiter and Saturn at their present masses and in their present orbits. Gas drag introduces a size-dependent phasing of the secular perturbations, which leads to a pronounced dip in encounter velocities (Venc) between bodies of similar mass. Plantesimals of identical mass have Venc approximately 1 and approximately 10 m s-1 (near 1 and 2.6 AU, respectively) while bodies differing by approximately 10 in mass have Venc approximately 10 and approximately 100 m s-1 (near 1 and 2.6 AU, respectively). Under these conditions, growth, rather than erosion, will occur only by collisions of bodies of nearly the same mass. There will be essentially no gravitational focusing between bodies less than 10(22) to 10(25) g, allowing growth of planetary embryos in the terrestrial planet region to proceed in a slower nonrunaway fashion. The environment in the asteroid belt will be even more forbidding and it is uncertain whether even the severely depleted present asteroid belt could form under these conditions. The perturbations of Jupiter and Saturn are quite sensitive to their semi-major axes and decrease when the planets' heliocentric distances are increased to allow for protoplanet migration. It is possible, though not clearly demonstrated, that this could produce a depleted asteroid belt but permit formation of a system of terrestrial planet embryos on a approximately 10(6)-year timescale, initially by nonrunaway growth and transitioning to runaway growth after approximately 10(5) years. The calculations reported here are valid under the condition that the relative velocities of the bodies are determined only by Jupiter and Saturn perturbations and by gas drag, with no mutual perturbations between planetesimals. If, while subject to these conditions, the bodies become large enough for their mutual perturbations to influence their velocity and size evolution significantly, the problem becomes much more complex. This problem is under investigation.

NASA Discipline Exobiology↗

Role for Lower Extremity Interstitial Fluid Volume Changes in the Development of Orthostasis after Simulated Microgravity

Reentry orthostasis after exposure to the conditions of spaceflight is a persistent problem among astronauts. In a previous study, a computer model systems analysis was used to examine the physiologic mechanisms involved in this phenomenon. In this analysis, it was determined that an augmented capacitance of lower extremity veins due to a fluid volume contracture of the surrounding interstitial spaces during spaceflight results in an increase in sequestered blood volume upon standing and appears to be the initiating mechanism responsible for reentry orthostasis. In this study, we attempt to validate the central premise of this hypothesis using a ground-based spaceflight analog. 10 healthy subjects were placed at bed rest in a 6 head down tilt position for 60 days of bed rest. The impact of adaptations in interstitial fluid volume and venous capacitance in the lower extremities were then observed during a standard tilt test protocol performed before and after the confinement period. The interstitial thickness superficial to the calcaneous immediately below the lateral malleolus was measured using ultrasound with a 17-5 MHz linear array transducer. Measurements of the changes in anterior tibial vein diameter during tilt were obtained by similar methods. The measurements were taken while the subjects were supine and then during upright tilt (80') for thirty minutes, or until the subject had signs of presyncope. Additional measurements of the superficial left tibia interstitial thickness and stroke volume by standard echocardiographic methods were also recorded. In addition, calf compliance was measured over a pressure range of 10-60 mmHg, using plethysmography, in a subset of these subjects (n = 5). There was a average of 6% diminution in the size of the lower extremity interstitial space as compared to measurements acquired prior to bed rest. This contracture of the interstitial space coincided with a subsequent relative increase in the percentage change in tibial vein diameter and stroke volume upon tilting in contrast to the observations made before bed rest (54 vs 23% respectively). Compliance in the calf increased by an average of 36% by day 27 of bedrest. A systems analysis using a computer model of cardiovascular physiology suggests that microgravity induced interstitial volume depletion results in an accentuation of venous blood volume sequestration and is the initiating event in reentry orthostasis. This hypothesis was tested in volunteer subjects using a ground-based spaceflight analog model that simulated the body fluid redistribution induced by microgravity exposure. Measurements of changes in the interstitial spaces and observed responses of the anterior tibial vein with tilt, together with the increase in calf compliance, were consistent with our proposed mechanism for the initiation of postflight orthostasis often seen in astronauts.

Platts, Steven H.↗

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗