Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 487 records · Page 27

Update and Status of the Mars Climate Modeling Center at NASA Ames Research Center

NASA Ames Research Center has a long-established history in the study and numerical modeling of Mars’ atmosphere and climate. Such studies began with the late Jim Pollack and have continued under the leaderships of Bob Haberle, Jeff Hollingsworth, and now Melinda Kahre. Our group has grown to more than 10 members, including civil servants, research scientists, postdocs, and students. While we still prioritize science and the model development that supports it, we have recently enhanced efforts to make our codes and output publicly available and boost community engagement through the hosting of modeling tutorials, etc. Here our goal is to present the status of the Ames Mars GCM, our ongoing science projects, the tools we have recently made publicly available, and our plans for releases and continued community engagement.

M. A. Kahre↗

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗

Automated generation of influence functions for planar crack problems

A numerical procedure for the generation of influence functions for Mode I planar problems is described. The resulting influence functions are in a form for convenient evaluation of stress-intensity factors for complex stress distributions. Crack surface displacements are obtained by a least-squares solution of the Williams eigenfunction expansion for displacements in a cracked body. Discrete values of the influence function, evaluated using the crack surface displacements, are curve fit using an assumed functional form. The assumed functional form includes appropriate limit-behavior terms for very deep and very shallow cracks. Continuous representation of the influence function provides a convenient means for evaluating stress-intensity factors for arbitrary stress distributions by numerical integration. The procedure is demonstrated for an edge-cracked strip and a radially cracked disk. Comparisons with available published results demonstrate the accuracy of the procedure.

Sire, Robert A.↗

Discrete-continuous variable structural synthesis using dual methods

Approximation concepts and dual methods are extended to solve structural synthesis problems involving a mix of discrete and continuous sizing type of design variables. Pure discrete and pure continuous variable problems can be handled as special cases. The basic mathematical programming statement of the structural synthesis problem is converted into a sequence of explicit approximate primal problems of separable form. These problems are solved by constructing continuous explicit dual functions, which are maximized subject to simple nonnegativity constraints on the dual variables. A newly devised gradient projection type of algorithm called DUAL 1, which includes special features for handling dual function gradient discontinuities that arise from the discrete primal variables, is used to find the solution of each dual problem. Computational implementation is accomplished by incorporating the DUAL 1 algorithm into the ACCESS 3 program as a new optimizer option. The power of the method set forth is demonstrated by presenting numerical results for several example problems, including a pure discrete variable treatment of a metallic swept wing and a mixed discrete-continuous variable solution for a thin delta wing with fiber composite skins.

Schmit, L. A.↗

Torque Balances on the Taylor Cylinders in the Geomagnetic Data Assimilation

In this presentation we report on our continuing effort in geomagnetic data assimilation, aiming at understanding and predicting geomagnetic secular variation on decadal time scales. In particular, we focus on the effect of the torque balances on the cylindrical surfaces in the core co-axial with the Earth's rotation axis (the Taylor cylinders) on the time evolution of assimilated solutions. We use our MoSST core dynamics,model and observed geomagnetic field at the Earth's surface derived via Comprehensive Field Model (CFM) for the geomagnetic data assimilation. In our earlier studies, a model solution is selected randomly from our numerical database. It is then assimilated with the observations such that the poloidal field possesses the same field tomography on the core-mantel boundary (CMB) continued downward from surface observations. This tomography change is assumed to be effective through out the outer core. While this approach allows rapid convergence between model solutions and the observations, it also generates sevee numerical instabilities: the delicate balance between weak fluid inertia and the magnetic torques on the Taylor cylinders are completely altered. Consequently, the assimilated solution diverges quickly (in approximately 10% of the magnetic free-decay time in the core). To improve the assimilation, we propose a partial penetration of the assimilation from the CMB: The full-scale modification at the CMB decreases linearly and vanish at an interior radius r(sub a). We shall examine from our assimilation tests possible relationships between the convergence rate of the model solutions to observations and the cut-off radius r(sub a). A better assimilation shall serve our nudging tests in near future.

Kuang, Weijia↗

Chart system simplifies identification of complex design assemblies

Identification breakdown chart that lists the component parts required for any specific end item is used to identify rapidly and accurately, from numerous drawings, all the component parts of a complex design assembly. Cylindrical and complex configurations are depicted as continuous flat surfaces for ready identification.

Morin, H. P.↗

Numerical analyses for treating diffusion in single-, two- and three-phase binary alloy systems

Numerical solutions were applicable for planar, cylindrical, or spherical geometries with any diffusion-zone size and any continuous variation of the diffusion coefficient with concentration. Special techniques were included in the analyses to account for differences in molal volumes, initiation and growth of an intermediate phase, disappearance of a phase, and the presence of an initial composition profile in the specimen. A major improvement in solution accuracy was achieved in the two phase analysis by employing a mass conservation criterion to establish the location of the interface rather than the conventional interface-flux-balance criterion. In the three phase analysis, computation time was minimized without sacrificing solution accuracy by treating the three phase problem as a two phase problem when the thickness of the intermediate phase was less than a preset small value. Three computer codes were developed to perform these analyses.

Tenney, D. R.↗

Study of hypersonic inlet flow fields with a three-dimensional Navier-Stokes code

A three-dimensional Navier-Stokes code is used to study the flow through an inlet with a rectangular capture area, but with circular or essentially circular internal cross-sections. The inlet may have a center body extending partially or fully over the length of the inlet. Calculations are made under a variety of geometrical constraints imposed on the inlet cross-sections to study the flow quality in the inlet. It is observed that for this class of inlets, the boundary layer has a tendency to thicken and separate readily due to the adverse pressure gradient caused by the continuous compression in the inlet. Since no experimental results are yet available, this paper has presented only the numerical results from the present analysis.

Kumar, A.↗

A new model of cometary ionospheres

The coupled continuity, momentum, and energy equations were solved for ionospheric conditions appropriate for Comet Halley at 1 AU. The numerical scheme used is such that any shock transition appears naturally in the solution and no a priori assumptions are necessary. Solutions were obtained for a number of different assumptions concerning electron heating rates, but all showed that the electron temperatures increase rapidly and significantly at a distance from the nucleus where collisional electron-neutral cooling becomes unimportant. This temperature increase is accompanied by a sharp increase in both the plasma pressure and its associated polarization electric field, causing the supersonic plasma flow to go subsonic. It is not clear at this time whether or not this sonic transition is accompanied by a shock.

Korosmezey, A.↗

Superbubbles in disk galaxies

Correlated supernovae from an OB association create a superbubble: a large, thin, shell of cold gas surrounding a hot pressurized interior. Because supernova blast waves usually become subsonic before reaching the walls of the shell or cooling radiatively, the energy input from supernovae may be reasonably approximated as a continuous luminosity. Using the Kompaneets (thin-shell) approximation, the growth of superbubbles in various stratified atmospheres is numerically modeled. A dimensionless quantity predicts whether a superbubble will blow out of the H I disk of a spiral galaxy (and begin to accelerate upward) or collapse. Superbubbles blow out of the H I layer when they have a radius in the plane between one and two scale heights. They blow out only one side of a disk galaxy if their centers are more than 50-60 p above the plane and the gas layer has density and scale height typical of the Milky Way. Fingers of warm interstellar gas intrude into the hot interior when the superbubble overtakes dense clouds.

Mac Low, Mordecai-Mark↗

Optimal placement of actuators and sensors in control augmented structural optimization

A control-augmented structural synthesis methodology is presented in which actuator and sensor placement is treated in terms of (0,1) variables. Structural member sizes and control variables are treated simultaneously as design variables. A multiobjective utopian approach is used to obtain a compromise solution for inherently conflicting objective functions such as strucutal mass control effort and number of actuators. Constraints are imposed on transient displacements, natural frequencies, actuator forces and dynamic stability as well as controllability and observability of the system. The combinatorial aspects of the mixed - (0,1) continuous variable design optimization problem are made tractable by combining approximation concepts with branch and bound techniques. Some numerical results for example problems are presented to illustrate the efficacy of the design procedure set forth.

Sepulveda, A. E.↗

Pitch angle diffusion of newborn ions due to intrinsic turbulence in the solar wind

The objective of the present study is to understand the interaction of the solar wind with newborn ions in far upstream regions of a comet where the level of intrinsic turbulence is moderately low. Based on the assumption that quasi-linear theory is adequate and applicable, the pitch angle diffusion process and the time evolution of the newborn ion distribution function are investigated. Numerical solutions to the quasi-linear diffusion equation, including the effect of resonance broadening and that of continuous creation of newborn ions due to the ionization process, are obtained under several assumptions and approximations. It is found that theoretical results are consistent with the Giotto observations recently reported by Neugebauer et al. (1989).

Ziebell, L. F.↗

The delayed-detonation model of Type Ia supernovae. 2: The detonation phase

The investigation, by use of two-dimensional numerical hydrodynamics simulations, of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is continued. Previously we found that the deflagration is insufficient to unbind the star. Expansion shuts off the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. The burning was allowed to develop into a detonation in these nonspherical models. The detonation grows toward spherical symmetry at late times. At these densities (rho approx. 10(exp 7) to 10(exp 8) g cm(exp -3)), either Ni-56 or nuclei of the Si-Ca group are the dominant products of the burning. The bulk yields are sensitive to the density of the star when the transition to detonation occurs. The relevance of the abundances, velocities, mixing, and total energy release to the theory and interpretation of Type Ia supernovae is discussed.

Arnett, David↗

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗

Development and Implementation of Real-Time Information Delivery Systems for Emergency Management

The disaster management community has an on-going need for real-time data and information, especially during catastrophic events. Currently, twin engine or jet aircraft with limited altitude and duration capabilities collect much of the data. Flight safety is also an issue. Clearly, much of the needed data could be delivered via over-the-horizon transfer through a uninhabited aerial vehicles (UAV) platform to mission managers at various locations on the ground. In fact, because of the ability to stay aloft for long periods of time, and to fly above dangerous situations, UAV's are ideally suited for disaster missions. There are numerous situations that can be considered disastrous for the human population. Some, such as fire or flood, can continue over a period of days. Disaster management officials rely on data from the site to respond in an optimum way with warnings, evacuations, rescue, relief, and to the extent possible, damage control. Although different types of disasters call for different types of response, most situations can be improved by having visual images and other remotely sensed data available. "Disaster Management" is actually made up of a number of activities, including: - Disaster Prevention and Mitigation - Emergency Response Planning - Disaster Management (real-time deployment of resources, during an event) - Disaster / Risk Modeling All of these activities could benefit from real-time information, but a major focus for UAV-based technology is in real-time deployment of resources (i.e., emergency response teams), based on changing conditions at the location of the event. With all these potential benefits, it is desirable to demonstrate to user agencies the ability to perform disaster management missions as described. The following demonstration project is the first in a program designed to prove the feasibility of supporting disaster missions with UAV technology and suitable communications packages on-board. A several-year program is envisioned, in which a broad range of disaster-related activities are demonstrated to the appropriate user communities.

Wegener, Steve↗

Multimodal Neuroelectric Interface Development

This project aims to improve performance of NASA missions by developing multimodal neuroelectric technologies for augmented human-system interaction. Neuroelectric technologies will add completely new modes of interaction that operate in parallel with keyboards, speech, or other manual controls, thereby increasing the bandwidth of human-system interaction. We recently demonstrated the feasibility of real-time electromyographic (EMG) pattern recognition for a direct neuroelectric human-computer interface. We recorded EMG signals from an elastic sleeve with dry electrodes, while a human subject performed a range of discrete gestures. A machine-teaming algorithm was trained to recognize the EMG patterns associated with the gestures and map them to control signals. Successful applications now include piloting two Class 4 aircraft simulations (F-15 and 757) and entering data with a "virtual" numeric keyboard. Current research focuses on on-line adaptation of EMG sensing and processing and recognition of continuous gestures. We are also extending this on-line pattern recognition methodology to electroencephalographic (EEG) signals. This will allow us to bypass muscle activity and draw control signals directly from the human brain. Our system can reliably detect P-rhythm (a periodic EEG signal from motor cortex in the 10 Hz range) with a lightweight headset containing saline-soaked sponge electrodes. The data show that EEG p-rhythm can be modulated by real and imaginary motions. Current research focuses on using biofeedback to train of human subjects to modulate EEG rhythms on demand, and to examine interactions of EEG-based control with EMG-based and manual control. Viewgraphs on these neuroelectric technologies are also included.

Trejo, Leonard J.↗

Bumper: A Tool for Analyzing Spacecraft Micrometeoroid and Orbital Debris Risk

“Bumper” is NASA’s computer program for analyzing spacecraft micrometeoroid and orbital debris (MMOD) risk. Bumper was developed in the late-1980s and has been continuously used and maintained since. The user base has grown from a few government entities to now include numerous commercial entities as well. The NASA Johnson Space Center (JSC) Hypervelocity Impact Technology (HVIT) Team is responsible for all aspects of the Bumper software. Bumper has been used to characterize MMOD risk on hundreds of spacecraft. All of the International Space Station (ISS) modules, visiting vehicles and numerous external components and systems have been analyzed. Bumper was used to analyze each of the Space Shuttle missions since STS-50. The Orion Multi-Purpose Crew Vehicle (MPCV) MMOD shielding is being developed using Bumper as well. Bumper has also been used on numerous telescopes (Hubble, James Webb, and Fermi Gamma-ray Space Telescopes), scientific probes (Stardust, New Horizons, Parker Solar Probe), and Earth observation satellites (Landsat, Joint Polar Satellite System). Bumper is also being used to analyze the micrometeoroid risk and support design of the Deep Space Gateway (DSG) and Mars Sample Return (MSR) missions. The HVIT Bumper Configuration Control Board (CCB) ensures that all changes to the code are approved, reviewed, and documented. Most of the changes are made to add new MMOD damage “ballistic limit equations” (BLEs). BLEs are typically added in response to completion of a hypervelocity impact (HVI) test series and development of an associated BLE. Other less frequent changes include updates of the debris or meteoroid environment models, feature enhancements, and feature retirement. Some BLEs are commercially sensitive and/or proprietary, so the CCB also manages code user-version control and software distribution. The current version – “Bumper 3” – is a FORTRAN executable that utilizes a 64-bit architecture. Bumper 3 has numerous features that make it a powerful tool for analyzing spacecraft MMOD risk. Bumper uses the latest orbital debris and micrometeoroid environment models. Bumper also easily processes large spacecraft geometry models, recognizes hidden surfaces, permits BLE assignment by name or number, and conducts quality checks of the spacecraft geometry model. Bumper 3 can also be used to estimate the effects of particle penetration through thin, high-standoff distance hardware components such as solar arrays and radiators. This is done using a special HVIT-developed technique know as the “3-Part Analysis.” The paper introduces the Bumper 3 MMOD risk analysis code and provides an example MMOD risk assessment showing Bumper’s role in the overall MMOD protection design process.

Lear, Dana M.↗

Students Planting Research Of Use To Space (SPROUTS) Networking and Engagement of Interns at Kennedy Space Center

Beginning in 2016, the Kennedy Space Center (KSC) internship program for life sciences was re-envisioned. Rather than hosting interns in the summer only with numerous innovative small projects, the current model is to offer internship opportunities every semester, so that research is continuous and sustained. Since 2016 we have had 64 internship opportunities for undergraduate and graduate students in space crop production and microgravity simulation. 52 separate individuals have participated in this program, with some participating in more than one term. While interns primarily worked with individual mentors on a main project, most also worked on group projects or developed smaller secondary projects on their own initiative. Additionally, interns participated in a weekly Lunch, Learn, and Discuss seminar series where they presented their research, learned from scientists at KSC, and engaged with members of the space life sciences community. When the face-to-face spring 2020 internship ended prematurely, the weekly seminar series continued as a virtual event, which led to the concept of engaging former interns. Current and former interns were invited to participate in the on-going virtual series, and 46 asked to be included on the mailing list. Our current seminar series involves remote seminars from contacts within and outside of the agency, and engages a large network of current scientists and former interns with weekly participation averaging ~40 people. Speakers include experts from industry, academia, and other space agencies working on topics of interest to KSC researchers. This series has also served as a brainstorming platform for KSC scientists to discuss new ideas and concepts and invite feedback. Remote work and virtual platforms have allowed engagement in a new way that will continue as future interns are integrated into this network of individuals interested in space research. This work is funded by NASA’s Space Biology Program.

Space Crop Production↗