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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Systematic generation of nonlinear discretized dynamic equilibrium equations of spinning cantilevers

General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.

El-Essawi, M.↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1977

Displacement incrementation in non-linear structural analysis by the self-correcting method

The prediction of nonlinear structural behavior by the finite element method wherein buckling does not occur has received considerable attention and, with it, reasonable success has been achieved. However, the post-buckling problem has been less actively pursued probably because of the inherent numerical difficulties encountered. This note reviews very briefly the numerical methods currently being used for pre- and post-buckling analysis and presents a self-correcting approach based on load and displacement incrementation which is shown to be efficient, reliable, and easy to program. Numerical solutions are presented which demonstrate the effectiveness of the method.

Haisler, W. E.↗
NASA NTRS1988

An initial investigation into methods of computing transonic aerodynamics sensitivity coefficients

Based upon the program outlined in the modified proposal, the initial effort has concentrated on developing the quasi-analytical approach for determining aerodynamic influence coefficients. In addition, in order to keep these initial efforts computationally efficient and straightforward, it was decided to first consider only two-dimensional flow and to model the problem using the transonic small perturbation equation with surface boundary conditions applied along the slit. The small perturbation equation was selected instead of the full potential equation because it permits the capturing of important nonlinear physical phenomena while being by comparison easy to handle from a coding standpoint. Consequently, the initial objectives were: (1) investigate the feasibility of using the quasi-analytical approach for determining aerodynamic sensitivity coefficients in the transonic regime; (2) determine the accuracy of the quasi-analytical approach by comparing results to those computed using finite difference techniques; and (3) investigate efficient methods for solving the large set of equations associated with the quasi-analytical approach.

Carlson, Leland A.↗
NASA NTRS1989

Model and simulation development for turbopump health monitoring by parameter estimation

The nonlinear equations describing the flow, pressure and temperature dynamics of the High Pressure Fuel Turbopump of the Space Shuttle Main Engine are summarized and implemented in a digital simulation. Some closed form approximations to the operating curves of the turbopump are made so that closed form derivatives exist. These equations and approximations are then used to develop an extended Kalman filter that estimates the dynamic states of the turbopump and a parameter that can reflect turbopump degradation. The filter is implemented as a digital computer program and applied to data generated by the digital simulation of the turbopump with substantial noise added to the output signals. The results indicate that the parameter estimate converges quickly and accurately. This suggests that the parameter estimate can be used to monitor for turbopump degradations affecting its value.

Walker, Bruce K.↗
NASA NTRS1990

A superlinear interior points algorithm for engineering design optimization

We present a quasi-Newton interior points algorithm for nonlinear constrained optimization. It is based on a general approach consisting of the iterative solution in the primal and dual spaces of the equalities in Karush-Kuhn-Tucker optimality conditions. This is done in such a way to have primal and dual feasibility at each iteration, which ensures satisfaction of those optimality conditions at the limit points. This approach is very strong and efficient, since at each iteration it only requires the solution of two linear systems with the same matrix, instead of quadratic programming subproblems. It is also particularly appropriate for engineering design optimization inasmuch at each iteration a feasible design is obtained. The present algorithm uses a quasi-Newton approximation of the second derivative of the Lagrangian function in order to have superlinear asymptotic convergence. We discuss theoretical aspects of the algorithm and its computer implementation.

Herskovits, J.↗
NASA NTRS2005

Twist Model Development and Results From the Active Aeroelastic Wing F/A-18 Aircraft

Understanding the wing twist of the active aeroelastic wing F/A-18 aircraft is a fundamental research objective for the program and offers numerous benefits. In order to clearly understand the wing flexibility characteristics, a model was created to predict real-time wing twist. A reliable twist model allows the prediction of twist for flight simulation, provides insight into aircraft performance uncertainties, and assists with computational fluid dynamic and aeroelastic issues. The left wing of the aircraft was heavily instrumented during the first phase of the active aeroelastic wing program allowing deflection data collection. Traditional data processing steps were taken to reduce flight data, and twist predictions were made using linear regression techniques. The model predictions determined a consistent linear relationship between the measured twist and aircraft parameters, such as surface positions and aircraft state variables. Error in the original model was reduced in some cases by using a dynamic pressure-based assumption and by using neural networks. These techniques produced excellent predictions for flight between the standard test points and accounted for nonlinearities in the data. This report discusses data processing techniques and twist prediction validation, and provides illustrative and quantitative results.

Lizotte, Andrew↗
DOE OSTI2022

Optimization of Superconducting Linac for Proton Improvement Plan-II (PIP-II)

PIP-II is an essential upgrade of the Fermilab complex that will enable the worlds most intense high-energy beam of neutrinos for the international Deep Underground Neutrino Experiment at LBNF and support a broad physics program at Fermilab. Ultimately, the PIP-II superconducting linac will be capable of accelerating the $H-$ CW beam to 800 MeV with an average power of 1.6 MW. To operate the linac with such high power, beam losses and beam emittance growth must be tightly controlled. In this paper, we present the results of global optimization of the Linac options towards a robust and efficient physics design for the superconducting section of the PIP-II linac. We also investigate the impact of the nonlinear field of the dipole correctors on the beam quality and derive the requirement on the field quality using statistical analysis. Finally, we assess the need to correct the quadrupole focusing produced by Half Wave, and Single Spoke accelerating cavities. We assess the feasibility of controlling the beam coupling in the machine by changing the polarity of the field of linac focusing solenoids

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

The 2020 Nonlinear Mechanics and Dynamics Research Institute

The 2020 Nonlinear Mechanics and Dynamics (NOMAD) Research Institute was successfully held from June 15 to July 30, 2020. 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 11 students participated in the seven-week long program held virtually due to the COVID-19 health pandemic. The students collaborated on one of four research projects that were developed by various mentors from Sandia National Laboratories, the University of New Mexico, and other academic and research institutions. In addition to the research activities, the students attended weekly technical seminars, various virtual tours, and socialized at virtual gatherings. 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 2020 are published as proceedings at technical conferences and have direct alignment with the critical mission work performed at Sandia.

42 ENGINEERING↗
NASA NTRS1991

On the theory of the synthesis of offset dual-shaped reflectors - Case examples

In an earlier paper by Galindo-Israel et al. (1987), the geometrical optics (GO) principles, constraints, and requirements of the dual- and single-offset-shaped reflector synthesis problem were collected and developed into a set of nonlinear first-order PDEs. An extension of the methods by which solutions to the PDEs can be obtained is presented, together with several case examples. These examples are independently analyzed by GO and physical optics diffraction methods. The starting point for the integration over each reflector can be taken on the outer rim, at the center, or at an intermediate point-the intermediate starting point being the more general case. The utility of the speed of this synthesis method is demonstrated. For example, the program utilized completes a synthesis in less than 45 s on a 386/20 personal computer and within a fraction of a second on a mainframe Cray. This makes practical the incorporation of the synthesis into a search algorithm that can optimize one or more parameters of the reflector system. As an example, the optimization of the mapping equations for low cross polarization is discussed.

Galindo-Israel, Victor↗
NASA NTRS1977

Ionosphere/microwave beam interaction study

A solar power satellite microwave power density of 20mw sq cm was confirmed as the level where nonlinear interactions may occur in the ionosphere, particularly at 100 km altitude. Radio wave heating at this altitude, produced at the Arecibo Observatory, yielded negative results for radio wave heating of an underdense ionosphere. Overdense heating produced striations in the ionosphere which may cause severe radio frequency interference problems under certain conditions. The effects of thermal self-focusing are shown to be limited severely geographically. The aspect sensitivity of field-aligned striations makes interference-free regions above magnetic latitude about 60 deg. A test program is proposed to simulate the interaction of the SPS beam with the ionosphere, to measure the effects of the interaction on the ionosphere and on communication and navigation systems, and to interpret the results.

Duncan, L. M.↗
NASA NTRS2005

The Role of the Heavy Ions in the Generation of EMIC Waves

The effect of EMIC waves, generated by a positive ion temperature anisotropy on Earth s RC dynamics is one of the best known examples of wave-particle interaction in the magnetosphere and the most controversial mechanism of RC losses. Under certain conditions, relativistic electrons with energy greater than or equal to 1 MeV can be removed from the outer RB by EMIC wave scattering during a magnetic storm much faster than by any other loss mechanisms. That is why the calculation of EMIC waves are very critical part of the LWS program. Systematic studies of magnetosphere-plasmasphere-ionosphere coupling are needed in order to provide EMIC waves forecast on a global scale and include WPI processes in the RB modeling. To quantify the EMIC wave effects on the RC-ion and RB-electron dynamics, a self-consistent theoretical description of the ions, electrons, and EMIC waves should be employed in future RB studies. This talk will emphasize the role of the heavy ions in the EMIC waves formation on a global scale and their nonlinear coupling with lower hybrid waves.

Khazanov, G. V.↗
NASA NTRS2006

Processing TES Level-2 Data

TES Level 2 Subsystem is a set of computer programs that performs functions complementary to those of the program summarized in the immediately preceding article. TES Level-2 data pertain to retrieved species (or temperature) profiles, and errors thereof. Geolocation, quality, and other data (e.g., surface characteristics for nadir observations) are also included. The subsystem processes gridded meteorological information and extracts parameters that can be interpolated to the appropriate latitude, longitude, and pressure level based on the date and time. Radiances are simulated using the aforementioned meteorological information for initial guesses, and spectroscopic-parameter tables are generated. At each step of the retrieval, a nonlinear-least-squares- solving routine is run over multiple iterations, retrieving a subset of atmospheric constituents, and error analysis is performed. Scientific TES Level-2 data products are written in a format known as Hierarchical Data Format Earth Observing System 5 (HDF-EOS 5) for public distribution.

Poosti, Sassaneh↗
NASA NTRS1994

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

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

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

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

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