Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NONLINEAR PROGRAMMING”

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 1,693 records · Page 94

Guest investigator program study: Physics of equatorial plasma bubbles

Plasma bubbles are large-scale (10 to 100 km) depletions in plasma density found in the night-time equatorial ionosphere. Their formation has been found to entail the upward transport of plasma over hundreds of kilometers in altitude, suggesting that bubbles play significant roles in the physics of many of the diverse and unique features found in the low-latitude ionosphere. In the simplest scenario, plasma bubbles appear first as perturbations in the bottomside F layer, which is linearly unstable to the gravitationally driven Rayleigh-Taylor instability. Once initiated, bubbles develop upward through the peak of the F layer into its topside (sometimes to altitudes in excess of 1000 km), a behavior predicted by the nonlinear form of the same instability. While good general agreement has been found between theory and observations, little is known about the detailed physics associated with plasma bubbles. Our research activity centered around two topics: the shape of plasma bubbles and associated electric fields, and the day-to-day variability in the occurrence of plasma bubbles. The first topic was pursued because of a divergence in view regarding the nonlinear physics associated with plasma bubble development. While the development of perturbations in isodensity contours in the bottomside F layer into plasma bubbles is well accepted, some believed bubbles to be cylinder-like closed regions of depleted plasma density that floated upward leaving a turbulent wake behind them (e.g., Woodman and LaHoz, 1976; Ott, 1978; Kelley and Ott, 1978). Our results, summarized in a paper submitted to the Journal of Geophysical Research, consisted of incoherent scatter radar measurements that showed unambiguously that the depleted region is wedgelike and not cylinderlike, and a case study and modeling of SM-D electric field instrument (EFI) measurements that showed that the absence of electric-field perturbations outside the plasma-depleted region is a distinct signature of wedge-shaped plasma bubbles. The second topic was pursued because the inability to predict the day-to-day occurrence of plasma bubbles indicated inadequate knowledge of the physics of plasma bubbles. An understanding of bubble formation requires an understanding of the roles of the various terms in the linearized growth rate of the collisional Rayleigh-Taylor instability. In our study, we examined electric-field perturbations found in SM-D EFI data and found that the seeding is more likely to be produced in the E region rather than the F region. The results of this investigation are presented in the Appendix of this report and will be submitted for publication in the Journal of Geophysical Research.

Tsunoda, Roland T.↗

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

A study of prediction methods for the high angle-of-attack aerodynamics of straight wings and fighter aircraft

Work is described dealing with two areas which are dominated by the nonlinear effects of vortex flows. The first area concerns the stall/spin characteristics of a general aviation wing with a modified leading edge. The second area concerns the high-angle-of-attack characteristics of high performance military aircraft. For each area, the governing phenomena are described as identified with the aid of existing experimental data. Existing analytical methods are reviewed, and the most promising method for each area used to perform some preliminary calculations. Based on these results, the strengths and weaknesses of the methods are defined, and research programs recommended to improve the methods as a result of better understanding of the flow mechanisms involved.

Mcmillan, O. J.↗

Research on the Problem of High-Precision Deployment for Large-Aperture Space-Based Science Instruments

The present paper summarizes results from an ongoing research program conducted jointly by the University of Colorado and NASA Langley Research Center since 1994. This program has resulted in general guidelines for the design of high-precision deployment mechanisms, and tests of prototype deployable structures incorporating these mechanisms have shown microdynamically stable behavior (i.e., dimensional stability to parts per million). These advancements have resulted from the identification of numerous heretofore unknown microdynamic and micromechanical response phenomena, and the development of new test techniques and instrumentation systems to interrogate these phenomena. In addition, recent tests have begun to interrogate nanomechanical response of materials and joints and have been used to develop an understanding of nonlinear nanodynamic behavior in microdynamically stable structures. The ultimate goal of these efforts is to enable nano-precision active control of micro-precision deployable structures (i.e., active control to a resolution of parts per billion).

Lake, Mark S.↗

A Nonlinear Propulsion System Simulation Technique for Piloted Simulators

In the past, propulsion system simulations used in flight simulators have been extremely simple. This resulted in a loss of simulation realism since significant engine and aircraft interactions were neglected and important internal engine parameters were not computed. More detailed propulsion system simulators are needed to permit evaluations of modern aircraft propulsion systems in a simulated flight environment. A real time digital simulation technique has been developed which provides the capabilities needed to evaluate propulsion system performance and aircraft system interaction on manned flight simulators. A parameter correlation technique is used with real and pseudo dynamics in a stable integration convergence loop. The technique has been applied to a multivariable propulsion system for use in a piloted NASA flight simulator program. Cycle time is 2.0 ms on a Univac 1110 computer and 5.7 ms on the simulator computer, a Xerox Sigma 8. The model is stable and accurate with time steps up to 50 ms. The program evaluated the simulation technique and the propulsion system digital control. The simulation technique and model used in that program are described and results from the simulation are presented.

Mihaloew, J. R.↗

Simultaneous Observations of Surface Layer Profiles of Humidity, Temperature, and Wind Using Scanning Lidar Instruments

We demonstrate the combination of three Doppler lidars (DLs), a water vapor differential absorption lidar and a temperature rotational Raman lidar for the investigation of the interactions between the land-surface and the atmospheric boundary layer. This combination of scanning lidars was operated for the first time during the Land-Atmosphere Feedback Experiment at the Atmospheric Radiation Measurement program's Southern Great Plain site, Oklahoma, USA, in August 2017, and provided simultaneous surface layer profiles of horizontal wind, humidity and temperature. The horizontal wind profiles were determined using the dual-Doppler method with two DLs. The scans were performed above four towers providing atmospheric variables and fluxes. These combined lidar data allowed for the estimation of the friction velocity as well as the surface latent and sensible heat fluxes. For this purpose, profiles calculated with Monin-Obukhov similarity theory were fitted with the Levenberg-Marquardt nonlinear least squares curve-fitting method to the measured surface-layer lidar profiles. We present case studies of three 50-min periods to illustrate the new method. The measurements resulted in friction velocities of 0.31, 0.29, and 0.38 ms –1 , sensible surface heat fluxes of 312.1, 234.1, and 183.3 Wm –2 and latent surface heat fluxes of 251.8, 227.3, and 274.4 Wm –2 . These values were compared with the in-situ measurements of the towers. Considering the sampling differences and the error analyses, the results agree within 25.1%, 42.5%, and 28.8%, respectively, demonstrating that this synergy of scanning active remote sensing systems can be used to derive surface fluxes with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗

Microwave Imaging Reflectometry for the study of Edge Harmonic Oscillations on DIII-D [Microwave Imaging Reflectometry (MIR) for the study of Edge Harmonic Oscillations (EHOs) on DIII-D]

Quiescent H-mode (QH) is an ELM free mode of operation in which edge-localized harmonic oscillations (EHOs) are believed to enhance particle transport, thereby stabilizing ELMs and preventing damage to the divertor and plasma facing components. Microwave Imaging Reflectometer (MIR) enabling direct comparison between the measured and simulated 2D images of density fluctuations near the edge can determine the 2D structure of density oscillation which can help to explain the physics behind EHO modes. MIR data sometimes indicates a counter-propagation between higher (n>1) and dominant (n=1) harmonics of coherent EHOs in the steep gradient regions of the pedestal. To preclude diagnostic artifacts, we have performed forward modeling that includes possible optical misalignments to show that offsets between transmitting and receiving antennas do not account for this feature. We have also simulated the non-uniform rotation of the EHO structure, which induces multiple harmonics that are properly characterized in the synthetic diagnostic. Excluding these possible explanations for the data, the counter-propagation observed in MIR data, which is not corroborated by external Mirnov coil array measurements, may be due to subtleties of the eigenmode structure, such as an inversion radius consistent with a magnetic island. Similar effects are observed in analysis of internal ECE-Imaging and BES data. Furthermore, the identification of a non-ideal structure motivates further exploration of nonlinear models of this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Registering Articulated Objects With Human-in-the-loop Corrections

Remotely programming robots to execute tasks often relies on registering objects of interest in the robot’s environment. Frequently, these tasks involve articulating objects such as opening or closing a valve. However, existing human-in-the-loop methods for registering objects do not consider articulations and the corresponding impact to the geometry of the object, which can cause the methods to fail. In this work, we present an approach where the registration system attempts to automatically determine the object model, pose, and articulation for user-selected points using a nonlinear iterative closest point algorithm. When the automated fitting is incorrect, the operator can iteratively intervene with corrections after which the system will refit the object. We present an implementation of our fitting procedure and evaluate it with a user study that shows that it can improve user performance, in measures of time on task and task load, ease of use, and usefulness compared to a manual registration approach. We also present a situated example that demonstrates the integration of our method in an end-to-end system for articulating a remote valve.

teleoperation↗

On current aspects of finite element computational fluid mechanics for turbulent flows

A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Baker, A. J.↗

The development of advanced technology blades for tilt-rotor aircraft

The paper discusses the development and ground testing of blades for the XV-15 tilt-rotor demonstrator aircraft. This work was performed under contract NAS2-11250 with NASA Ames Research Center. These blades, known as the Advanced Technology Blades (ATB), replace the rectangular, steel blades which were part of the XV-15 original design. The materials used in the primary structure of the ATB are fiberglass and high strain graphite epoxy laminates. This facilitates the use of 43 deg of nonlinear twist, a nonuniform tapered planform and thin airfoils required for aerodynamic efficiency. Instrumentation life is extended by encapsulating gages and wiring in the composite structure. Tip shells and cuff fairings are removable to provide access to tip weights and retention hardware; they are also replaceable with alternate research configurations. Extensive laboratory testing has validated predicted strength characteristics. Hover testing has demonstrated performance significantly superior to that predicted by contemporary methodology. Key elements of the test rig used for rotor performance measurement were developed as an ancillary part of the present program. The performance testing included measurement of near- and far-field noise. Induced inflow velocity distributions were also determined and photographs of tip vortex condensation trails were taken. These are providing guidance for modifications to hover peformance codes.

Alexander, Harold R.↗

Life assessment of structural components using inelastic finite element analyses

The need for enhanced and improved performance of structural components subject to severe cyclic thermal/mechanical loadings, such as in the aerospace industry, requires development of appropriate solution technologies involving time-dependent inelastic analyses. Such analyses are mandatory to predict local stress-strain response and to assess more accurately the cyclic life time of structural components. The NASA-Lewis Research Center is cognizant of this need. As a result of concerted efforts at Lewis during the last few years, several such finite element solution technologies (in conjunction with the finite element program MARC) were developed and successfully applied to numerous uniaxial and multiaxial problems. These solution technologies, although developed for use with MARC program, are general in nature and can easily be extended for adaptation with other finite element programs such as ABAQUS, ANSYS, etc. The description and results obtained from two such inelastic finite element solution technologies are presented. The first employs a classical (non-unified) creep-plasticity model. An application of this technology is presented for a hypersonic inlet cowl-lip problem. The second of these technologies uses a unified creep-plasticity model put forth by Freed. The structural component for which this finite element solution technology is illustrated, is a cylindrical rocket engine thrust chamber. The advantages of employing a viscoplastic model for nonlinear time-dependent structural analyses are demonstrated. The life analyses for cowl-lip and cylindrical thrust chambers are presented. These analyses are conducted by using the stress-strain response of these components obtained from the corresponding finite element analyses.

Arya, Vinod K.↗

A Comparison of PETSC Library and HPF Implementations of an Archetypal PDE Computation

Two paradigms for distributed-memory parallel computation that free the application programmer from the details of message passing are compared for an archetypal structured scientific computation a nonlinear, structured-grid partial differential equation boundary value problem using the same algorithm on the same hardware. Both paradigms, parallel libraries represented by Argonne's PETSC, and parallel languages represented by the Portland Group's HPF, are found to be easy to use for this problem class, and both are reasonably effective in exploiting concurrency after a short learning curve. The level of involvement required by the application programmer under either paradigm includes specification of the data partitioning (corresponding to a geometrically simple decomposition of the domain of the PDE). Programming in SPAM style for the PETSC library requires writing the routines that discretize the PDE and its Jacobian, managing subdomain-to-processor mappings (affine global- to-local index mappings), and interfacing to library solver routines. Programming for HPF requires a complete sequential implementation of the same algorithm, introducing concurrency through subdomain blocking (an effort similar to the index mapping), and modest experimentation with rewriting loops to elucidate to the compiler the latent concurrency. Correctness and scalability are cross-validated on up to 32 nodes of an IBM SP2.

Hayder, M. Ehtesham↗

A Survey Of Earth-Moon Libration Orbits: Stationkeeping Strategies And Intra-Orbit Transfers

Cislunar space is a readily accessible region that may well develop into a prime staging area in the effort to colonize space near Earth or to colonize the Moon. While there have been statements made by various NASA programs regarding placement of resources in orbit about the Earth-Moon Lagrangian locations, there is no survey of the total cost associated with attaining and maintaining these unique orbits in an operational fashion. Transfer trajectories between these orbits required for assembly, servicing, and positioning of these resources have not been extensively investigated. These orbits are dynamically similar to those used for the Sun-Earth missions, but differences in governing gravitational ratios and perturbation sources result in unique characteristics. We implement numerical computations using high fidelity models and linear and nonlinear targeting techniques to compute the various maneuver (Delta)V and temporal costs associated with orbits about each of the Earth-Moon Lagrangian locations (L1, L2, L3, L4, and L5). From a dynamical system standpoint, we speak to the nature of these orbits and their stability. We address the cost of transfers between each pair of Lagrangian locations.

Folta, David↗

Application of Variance-Based Sensitivity Analysis to a Large System Dynamics Model

Variance-based sensitivity methods can provide insights into large computational models. We present a novel application of sensitivity analysis to the Biomass Scenario Model (BSM) a large and complex system dynamics model of the developing biofuels industry in the United States. We apply a two-stage sensitivity approach consisting of an initial sensitivity screening, followed by a variance decomposition approach. Identifying key system levers and quantifying their strength is not straightforward in complex system dynamics models that have numerous feedbacks and nonlinear results. Variance-based sensitivity analysis (VBSA) offers a systematic, global approach to assessing system dynamics models because it addresses nonlinear responses and interactive effects. Especially when a large model's size makes manual exploration of the input space difficult and time-consuming, the approach can help to provide a comprehensive understanding of interactions that drive model behaviors.

97 MATHEMATICS AND COMPUTING↗

The 2019 Nonlinear Mechanics and Dynamics Research Institute

The 2019 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 17 to August 1, 2019. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 20 students came to Albuquerque, New Mexico to participate in the seven-week long program held at the Mechanical Engineering building on the University of New Mexico campus. The students collaborated on one of seven research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and academic institutions. In addition to the research activities, the students attended weekly technical seminars, various tours, and socialized at various off-hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗

The 2018 Nonlinear Mechanics and Dynamics Research Institute

The 2018 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 18 to August 2, 2018. NOMAD brings together participants with diverse technical backgrounds to work in small teams to cultivate new ideas and approaches in engineering mechanics and dynamics research. NOMAD provides an opportunity for researchers -- especially early career researchers - to develop lasting collaborations that go beyond what can be established from the limited interactions at their institutions or at annual conferences. A total of 17 students came to Albuquerque, New Mexico to participate in the seven-week long program held at the Mechanical Engineering building on the University of New Mexico campus. The students collaborated on one of six research projects that were developed by various mentors from Sandia National Laboratories, University of New Mexico, and academic institutions. In addition to the research activities, the students attended weekly technical seminars, various tours, and socialized at various off-hour events including an Albuquerque Isotopes baseball game. At the end of the summer, the students gave a final technical presentation on their research findings. Many of the research discoveries made at NOMAD are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗

Jim Starnes' Contributions to Residual Strength Analysis Methods for Metallic Structures

A summary of advances in residual strength analyses methods for metallic structures that were realized under the leadership of Dr. James H. Starnes, Jr., is presented. The majority of research led by Dr. Starnes in this area was conducted in the 1990's under the NASA Airframe Structural Integrity Program (NASIP). Dr. Starnes, respectfully referred to herein as Jim, had a passion for studying complex response phenomena and dedicated a significant amount of research effort toward advancing damage tolerance and residual strength analysis methods for metallic structures. Jim's efforts were focused on understanding damage propagation in built-up fuselage structure with widespread fatigue damage, with the goal of ensuring safety in the aging international commercial transport fleet. Jim's major contributions in this research area were in identifying the effects of combined internal pressure and mechanical loads, and geometric nonlinearity, on the response of built-up structures with damage. Analytical and experimental technical results are presented to demonstrate the breadth and rigor of the research conducted in this technical area. Technical results presented herein are drawn exclusively from papers where Jim was a co-author.

Young, Richard D.↗

Material Model Evaluation of a Composite Honeycomb Energy Absorber

A study was conducted to evaluate four different material models in predicting the dynamic crushing response of solid-element-based models of a composite honeycomb energy absorber, designated the Deployable Energy Absorber (DEA). Dynamic crush tests of three DEA components were simulated using the nonlinear, explicit transient dynamic code, LS-DYNA . In addition, a full-scale crash test of an MD-500 helicopter, retrofitted with DEA blocks, was simulated. The four material models used to represent the DEA included: *MAT_CRUSHABLE_FOAM (Mat 63), *MAT_HONEYCOMB (Mat 26), *MAT_SIMPLIFIED_RUBBER/FOAM (Mat 181), and *MAT_TRANSVERSELY_ANISOTROPIC_CRUSHABLE_FOAM (Mat 142). Test-analysis calibration metrics included simple percentage error comparisons of initial peak acceleration, sustained crush stress, and peak compaction acceleration of the DEA components. In addition, the Roadside Safety Verification and Validation Program (RSVVP) was used to assess similarities and differences between the experimental and analytical curves for the full-scale crash test.

Jackson, Karen E.↗