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 1,639 records · Page 91

Optimization-based calibration of hydrodynamic drag coefficients for a semisubmersible platform using experimental data of an irregular sea state

For the simulation of the coupled dynamic response of floating offshore wind turbines, it is crucial to calibrate the hydrodynamic damping with experimental data. The aim of this work is to find a set of hydrodynamic drag coefficients for the semisubmersible platform of the Offshore Code Comparison Collaboration, Continuation, with Correlation and unCertainity (OC6) project which provides suitable results for an irregular sea state. Due to the complex interaction of several degrees of freedom (DOF), it is common to calibrate drag coefficients with the time series of decay tests. However, applying these drag coefficients for the simulation of an irregular sea state results in misprediction of the motions. By using numerical optimization, it is possible to calibrate multiple drag coefficients simultaneously and effectively, while also considering several DOF. This work considers time series of structural displacements from wave tank tests of the OC6 project and from simulations of the same load cases in OpenFAST. Results are transferred into the frequency domain and the deviation between power spectral densities of surge, pitch and heave from experiment and numerical simulation is used as an objective function to obtain the best fitting drag coefficients. This novel numerical optimization approach enables finding one set of drag coefficients for different load cases, which is a major improvement compared to decay-test-tuned drag coefficients.

42 ENGINEERING↗

On the conservation of the first adiabatic invariant in perpendicular shocks

The conservation of the first adiabatic invariant of particle motion, the particle magnetic moment, at the front of perpendicular magnetosonic fast-mode shock waves is considered. Results of numerical simulations are presented which show that preshock and postshock values of the magnetic moment are on the average equal for particles transmitted by infinitesimally thin perpendicular magnetosonic fast-mode shocks, in contrast to extrapolations from adiabatic theory, which suggest that the first invariant is least likely to be conserved in this situation. The observed conservation is thus attributed to the continuity of the flux of total particle and field angular momentum across the shock front as a result of angular momentum conservation.

Pesses, M. E.↗

On the continuous dependence with respect to sampling of the linear quadratic regulator problem for distributed parameter systems

The convergence of solutions to the discrete or sampled time linear quadratic regulator problem and associated Riccati equation for infinite dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero (infinity) is established. Both the finite and infinite time horizon problems are studied. In the finite time horizon case, strong continuity of the operators which define the control system and performance index together with a stability and consistency condition on the sampling scheme are required. For the infinite time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary of delay system, and a flexible beam are presented and discussed.

Rosen, I. G.↗

Direct simulation for the instability and breakup of laminar liquid jets

A direct numerical simulation method is described for predicting the deformation of laminar liquid jets. In the present nonlinear direct simulation, the convective term, which has been discarded in past linear analyses by Rayleigh and others, is included in the hydrodynamic equations. It is shown that only by maintaining full complexity of the nonlinear surface tension term accurate drop formation can be predicted. The continuity and momentum equations in the transient form are integrated on an adaptive grid, conforming the jet and surface wave shape. The equations, which are parabolic in time and elliptic in space, are solved by a TVD scheme with characteristic flux splitting. The results of the present work are discussed and compared with available measurements and other analyses. The comparison shows that among the predictions, the current 1-D direct simulation results agree best with the experimental data. Furthermore, the computer time requirements are much (an order of magnitude) smaller than those of previously reported multidimensional analyses.

Chuech, S. G.↗

Optimal discrete-time dynamic output-feedback design - A w-domain approach

An alternative method for optimal digital control design is described in this paper. The method is based on the usage of the w-transform and has many attractive design features. One of these is its immediate connection with frequency loop-shaping techniques that are now popular and effective for multivariable control synthesis in continuous-time domain. Furthermore, any design algorithms originally developed for continuous-time systems can now be immediately extended to the discrete-time domain. The main results presented in this paper are the exact problem formulation and solution of an optimal discrete-time dynamic output-feedback design in the w-domain involving a quadratic performance index to random disturbances. In addition, necessary conditions for optimality are obtained for the numerical solution of the optimal output-feedback compensator design. A numerical example is presented illustrating its application to the design of a low-order dynamic compensator in a stability augmentation system of a commercial transport.

Ha, Cheolkeun↗

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗

RIACS/USRA

The Research Institute for Advanced Computer Science (RIACS) was established by the Universities Space Research Association (USRA) at the NASA Ames Research Center (ARC) on 6 June 1983. RIACS is privately operated by USRA, a consortium of universities with research programs in the aerospace sciences, under contract with NASA. The primary mission of RIACS is to provide research and expertise in computer science and scientific computing to support the scientific missions of NASA ARC. The research carried out at RIACS must change its emphasis from year to year in response to NASA ARC's changing needs and technological opportunities. A flexible scientific staff is provided through a university faculty visitor program, a post doctoral program, and a student visitor program. Not only does this provide appropriate expertise but it also introduces scientists outside of NASA to NASA problems. A small group of core RIACS staff provides continuity and interacts with an ARC technical monitor and scientific advisory group to determine the RIACS mission. RIACS activities are reviewed and monitored by a USRA advisory council and ARC technical monitor. Research at RIACS is currently being done in the following areas: Parallel Computing, Advanced Methods for Scientific Computing, High Performance Networks and Technology, and Learning Systems. Parallel compiler techniques, adaptive numerical methods for flows in complicated geometries, and optimization were identified as important problems to investigate for ARC's involvement in the Computational Grand Challenges of the next decade.

Oliger, Joseph↗

Time-domain approach for the transient responses in stratified viscoelastic Earth models

We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Hanyk, L.↗

EUV energy distributions of accretion disks in active galactic nuclei

We have calculated UV/EUV (300 A less than or = lambda less than or = 1500 A) continuous energy distributions of accretion disks in the centers of active galactic nuclei (AGNs) for disk luminosities in the range 0.1 L(sub Edd) and central masses ranging from 10(exp 8) solar mass to 10(exp 10) solar mass. The vertical gas pressure structure of the disk is obtained analytically, the temperature stratification and the resulting continuum radiation fields are calculated numerically. We show that weak Lyman edges are an intrinsic feature of such disks. The strength of the H I Lyman edge decreases for increasing accretion rate and fixed mass of the central black hole. It increases for increasing central mass and fixed luminosity in terms of the Eddington luminosity.

Storzer, H.↗

Autonomous Assembly of Modular Structures in Space and on Extraterrestrial Locations

The new U.S. National Vision for Space Exploration requires many new enabling technologies to accomplish the goals of space commercialization and returning humans to the moon and extraterrestrial environments. Traditionally, flight elements are complete subsystems requiring humans to complete the integration and assembly. These bulky structures also require the use of heavy launch vehicles to send the units to a desired location. This philosophy necessitates a high degree of safety, numerous space walks at a significant cost. Future space mission costs must be reduced and safety increased to reasonably achieve exploration goals. One proposed concept is the autonomous assembly of space structures. This concept is an affordable, reliable solution to in-space and extraterrestrial assembly. Assembly is autonomously performed when two components join after determining that specifications are correct. Local sensors continue monitor joint integrity post assembly, which is critical for safety and structural reliability. Achieving this concept requires a change in space structure design philosophy and the development of innovative technologies to perform autonomous assembly. Assembly of large space structures will require significant numbers of integrity sensors. Thus simple, low-cost sensors are integral to the success of this concept. This paper addresses these issues and proposes a novel concept for assembling space structures autonomously. Core technologies required to achieve in space assembly are presented. These core technologies are critical to the goal of utilizing space in a cost efficient and safe manner. Additionally, these novel technologies can be applied to other systems both on earth and extraterrestrial environments.

Alhorn, Dean C.↗

Pulse Detonation Engine Test Bed Developed

A detonation is a supersonic combustion wave. A Pulse Detonation Engine (PDE) repetitively creates a series of detonation waves to take advantage of rapid burning and high peak pressures to efficiently produce thrust. NASA Glenn Research Center's Combustion Branch has developed a PDE test bed that can reproduce the operating conditions that might be encountered in an actual engine. It allows the rapid and cost-efficient evaluation of the technical issues and technologies associated with these engines. The test bed is modular in design. It consists of various length sections of both 2- and 2.6- in. internal-diameter combustor tubes. These tubes can be bolted together to create a variety of combustor configurations. A series of bosses allow instrumentation to be inserted on the tubes. Dynamic pressure sensors and heat flux gauges have been used to characterize the performance of the test bed. The PDE test bed is designed to utilize an existing calorimeter (for heat load measurement) and windowed (for optical access) combustor sections. It uses hydrogen as the fuel, and oxygen and nitrogen are mixed to simulate air. An electronic controller is used to open the hydrogen and air valves (or a continuous flow of air is used) and to fire the spark at the appropriate times. Scheduled tests on the test bed include an evaluation of the pumping ability of the train of detonation waves for use in an ejector and an evaluation of the pollutants formed in a PDE combustor. Glenn's Combustion Branch uses the National Combustor Code (NCC) to perform numerical analyses of PDE's as well as to evaluate alternative detonative combustion devices. Pulse Detonation Engine testbed.

Breisacher, Kevin J.↗

National Report on the NASA Sounding Rocket and Balloon Programs

The U. S. National Aeronautics and Space Administration (NASA) Sounding Rockets and Balloon Programs conduct a total of 30 to 40 missions per year in support of the NASA scientific community and other users. The NASA Sounding Rockets Program supports the science community by integrating their experiments into the sounding rocket payloads, and providing both the rocket vehicle and launch operations services. Activities since 2011 have included two flights from Andoya Rocket Range, more than eight flights from White Sands Missile Range, approximately sixteen flights from Wallops Flight Facility, two flights from Poker Flat Research Range, and four flights from Kwajalein Atoll. Other activities included the final developmental flight of the Terrier-Improved Malemute launch vehicle, a test flight of the Talos-Terrier-Oriole launch vehicle, and a host of smaller activities to improve program support capabilities. Several operational missions have utilized the new Terrier-Malemute vehicle. The NASA Sounding Rockets Program is currently engaged in the development of a new sustainer motor known as the Peregrine. The Peregrine development effort will involve one static firing and three flight tests with a target completion data of August 2014. The NASA Balloon Program supported numerous scientific and developmental missions since its last report. The program conducted flights from the U.S., Sweden, Australia, and Antarctica utilizing standard and experimental vehicles. Of particular note are the successful test flights of the Wallops Arc Second Pointer (WASP), the successful demonstration of a medium-size Super Pressure Balloon (SPB), and most recently, three simultaneous missions aloft over Antarctica. NASA continues its successful incremental design qualification program and will support a science mission aboard WASP in late 2013 and a science mission aboard the SPB in early 2015. NASA has also embarked on an intra-agency collaboration to launch a rocket from a balloon to conduct supersonic decelerator tests. An overview of NASA's Sounding Rockets and Balloon Operations, Technology Development and Science support activities will be presented.

Eberspeaker, Philip↗

Experimental and Semiclassical Stark Widths and Shifts for Spectral Lines of Neutral and Ionized Atoms (A Critical Review of Experimental and Semiclassical Data for the Period 2008 through 2020)

This Review Article compares new experimental and corresponding semiclassical Stark broadened data for non-hydrogenic spectral lines of neutral atoms and positive ions. This Review covers the period 2008 until the end of 2020 and presents the continuation of previous critical reviews from 1976 (two), 1984 (two), 1990, 2002, and one from 2009. This Review reports the Stark broadening parameters of 1665 spectral lines belonging to 35 elements with 61 different atomic and ion species. The data are taken from 60 papers. Experimental data are arranged by elements and spectra, and these are reported in tables in alphabetical and numerical order, respectively. Each experimental Stark broadening parameter is followed with estimated accuracy. The experimental Stark broadening and shift data presented in the tables are compared with the corresponding semiclassical results available in the literature. In addition, for comparison experiments versus semiclassical theoretical Stark broadening parameters, a numerically improved theoretical approach developed based on Griem and co-workers’ theory is used to evaluate Stark widths and shifts of all studied lines whenever required atomic energy levels data and transition probabilities are available. At the end of the text report, for each analyzed neutral or ion species, the information about the location of the same species data in preceding reviews is given.

Plasma line broadening↗

Thermal Actuator Identification and Control for Thermomechanical Real-Time Cyber–Physical Testing

Thermomechanical cyber–physical testing enables two-way thermal coupling between a numerical and an experimental subsystem. The interactions between the numerical model and the physical specimen occur through transfer systems, which enforce interface conditions. Thus, efficient control methodologies are necessary to achieve the desired interface interaction through thermal actuators with minimal error. This study introduces a novel thermal transfer system that imposes distributed cooling (or heating) thermal loads on a physical subsystem. First, the thermal actuator is identified considering switching-mode continuous dynamics for heating and cooling conditions. A switching-mode estimation algorithm is adopted to estimate the operating thermal cycle of the actuator in real-time. A control system is developed to experimentally impose the desired temperature and reduce tracking error (i.e., the error between the desired and actual temperature) under different thermal cycles. The identification and control of the thermal transfer system are then validated through a set of experiments considering different temperature rates of change. The developed control system is found to effectively minimize tracking errors in real-time cyber–physical experiments.

Herta Montoya↗

Planetesimals to planets - Numerical simulation of collisional evolution

In a simulation of collisional and gravitational interaction in the early solar system, planets of approximately 500 km diameter are generated from an initial swarm of kilometer-sized planetesimals. Collisions are treated in accordance with experimental and theoretical impact results (such as rebound, cratering, and catastrophic fragmentation) for a variety of materials whose parameters span plausible values for early solid objects. In this model, the small planets form in approximately 10,000 yr; during this time, most of the mass of the system continues to reside in particles near the original size. It is thought that the few 500-km planets may act as 'seeds' for the subsequent gradual accretional growth into full-sized planets.

Greenberg, R.↗

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ∼12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

97 MATHEMATICS AND COMPUTING↗

Implications of an improved water equation of state for water-rich planets

Water (H 2 O), in all forms, is an important constituent in planetary bodies, controlling habitability and influencing geological activity. Under conditions found in the interior of many planets, as the pressure increases, the H-bonds in water gradually weaken and are replaced by ionic bonds. Recent experimental measurements of the water equation of state (EOS) showed both a new phase of H-bonded water ice, ice-VIIt, and a relatively low transition pressure just above 30 GPa to ionic bonded ice-X, which has a bulk modulus 2.5 times larger. The higher bulk modulus of ice-X produces larger planets for a given mass, thereby either reducing the atmospheric contribution to the volume of many exoplanets or limiting their water content. We investigate the impact of the new EOS measurements on the planetary mass–radius relation and interior structure for water-rich planets. We find that the change in the planet mass–radius relation caused by the systematic differences between previous and new experimental EOS measurements is comparable to the observational uncertainties in some planet sizes – an issue that will become more important as observations continue to improve.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data-Driven Adaptive Damping Controller for Wind Power Plants with Doubly-Fed Induction Generators

This paper presents an adaptive damping controller for wind power plants in which the turbines are equipped with doubly-fed induction generators. The controller is designed to respond to an input control signal that is triggered according to the system operating conditions. A processing unit continuously estimates the electromechanical modes of oscillation based on real-time streaming data acquired from a phasor measurement unit that is strategically positioned on the grid. The decision to trigger (or not trigger) the control signal is automatic, based on the relative damping of the dominant mode. The modes are estimated using the dynamic mode decomposition algorithm with time-delay embedding. Numerical simulations performed on the two-area system demonstrate that the proposed controller enhances the rotor angle stability for both small-signal and large disturbances, and is adaptive to changing grid conditions.

damping controller↗