Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear elimination”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Stability, accuracy, and efficiency of some underintegrated methods in finite element computations

In an attempt to increase computational efficiency in the numerical solution of highly nonlinear problems in solid and fluid mechanics, underintegrated finite element methods have been employed by many analysts. Underintegration refers to the use of a rule of an order lower than that required to integrate polynomial integrands exactly. The main drawback of this technique is related to the production of rank-deficient stiffness matrices, or equivalently an expanded kernel of the governing linear momentum operators. Such a development can introduce numerical instabilities. In order to overcome this difficulty, artificial stiffness or viscosity methods, or other stabilization methods have been proposed. One approach involves the elimination of spurious modes in a postprocessing operation. The present study is concerned with this a posteriori elimination method, taking into account the results which can be expected from it, and some of its possible extensions.

Jacquotte, O.-P.↗
DOE OSTI2024

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Image charge interaction correction in charged-defect calculations

Charged-defect calculation using a periodic supercell is a significant class of problems in solid state physics. However, the finite supercell size induces an undesirable long-range image charge Coulomb interaction. Although a variety of methods have been proposed to eliminate such image Coulomb interaction, most of the previous schemes are based on a rough approximation of the defect charge screening. In this work, we present a rigorous derivation of the image charge interaction with a defect screening model where the use of a bulk macroscopic dielectric constant can be avoided. Here, we have verified this approach in comparison with a widely used approach for 12 different defects. Our correction scheme offers a much faster convergence concerning the supercell size for cases with considerable image charge interactions. In those cases, we also found that the nonlinear dielectric screening might play an important role. Our proposed defect screening model will also shed new light on understanding the defect screening properties and can be applied to other defect systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

The use of laterally vectored thrust to counter thrust asymmetry in a tactical jet aircraft

A nonlinear, six degree-of-freedom flight simulator for a twin engine tactical jet was built on a hybrid computer to investigate lateral vectoring of the remaining thrust component for the case of a single engine failure at low dynamic pressures. Aircraft control was provided by an automatic controller rather than a pilot, and thrust vector control was provided by an open-loop controller that deflected a vane (located on the periphery of each exhaust jet and normally streamlined for noninterference with the flow). Lateral thrust vectoring decreased peak values of lateral control deflections, eliminated the requirement for steady-state lateral aerodynamic control deflections, and decreased the amount of altitude lost for a single engine failure.

Source record↗
NASA NTRS2002

Observer-Based Magnetic Bearing Controller Developed for Aerospace Flywheels

A prototype of a versatile, observer-based magnetic bearing controller for aerospace flywheels was successfully developed and demonstrated on a magnetic bearing test rig (see the photograph) and an actual flywheel module. The objective of this development included a fast, yet low risk, control development process, and a robust, high-performance controller for a large variety of flywheels. This required a good system model, an efficient development procedure, and a model-based controller that addressed the key problems associated with flywheel and bearing imbalance, sensor error, and vibration. The model used in this control development and tuning procedure included the flexible rotor dynamics and motor-induced vibrations. Such a model was essential for low-risk scheduling of speed-dependent control parameters and for reliable evaluation of novel control strategies. The successfully tested control prototype utilized an extended Kalman filter to estimate the true rotor principal-axis motion from the raw sensor position feedback. For control refinement, the extended Kalman filter also estimated and eliminated the combined effects of mass-imbalance and sensor runouts from the input data. A key advantage of the design based on the extended Kalman filter is its ability to accurately estimate both the rotor's principal-axis position and gyroscopic rates with the least amount of phase lag. This is important for control parameter scheduling to dampen the gyroscopic motions. Because of large uncertainties in the magnetic bearing and imbalance characteristics, this state-estimation scheme alone is insufficient for containing the rotor motion within the desired 1-mil excursion radius. A nonlinear gain adjustment based on an estimation of the principal-axis orbit size was needed to provide a coarse (nonoptimal), but robust, control of the orbit growth. Control current minimization was achieved with a (steepest gradient) search of synchronous errors in the principal-axis position input data. Actual flywheel tests of this observer-based controller (developed entirely in-house) at the NASA Glenn Research Center showed that the model correctly predicted the rotor orbit growth as a function of rotational speed, and it demonstrated the capability of gain adjustments to arrest this growth. Data from these tests on an actual flywheel module spun to 26,000 rpm proved that the controller was able to contain the shaft motion to within much less than 0.5 mils of radial excursion with axis currents less than 300 mA in root-mean-square estimate. The test speed range was limited because of thermal expansion concerns for this particular flywheel unit, not because of any deficiency in the controller. Simulations for this unit indicated that the controller should be robust up to its top operating speed of 60,000 rpm. Aside from these important achievements, and most significantly, it took less than 1 week to adapt this controller from the simple test rig to the actual flywheel and to demonstrate full five-axis levitation and control. This demonstration showed that both the controller and the model-based development and tuning framework are easily adaptable to a wide range of rotors and bearing configurations and, hence, are capable of reducing design risks and costs for many future flywheel technology developments.

Le, Dzu K.↗
DOE OSTI2023

Semiclassical theory and the Koopman-van Hove equation *

Abstract The phase space Koopman-van Hove (KvH) equation can be derived from the asymptotic semiclassical analysis of partial differential equations. Semiclassical theory yields the Hamilton–Jacobi equation for the complex phase factor and the transport equation for the amplitude. These two equations can be combined to form a nonlinear semiclassical version of the KvH equation in configuration space. There is a natural injection of configuration space solutions into phase space and a natural projection of phase space solutions onto configuration space. Hence, every solution of the configuration space KvH equation satisfies both the semiclassical phase space KvH equation and the Hamilton–Jacobi constraint. For configuration space solutions, this constraint resolves the paradox that there are two different conserved densities in phase space. For integrable systems, the KvH spectrum is the Cartesian product of a classical and a semiclassical spectrum. If the classical spectrum is eliminated, then, with the correct choice of Jeffreys–Wentzel–Kramers–Brillouin (JWKB) matching conditions, the semiclassical spectrum satisfies the Einstein–Brillouin–Keller quantization conditions which include the correction due to the Maslov index. However, semiclassical analysis uses different choices for boundary conditions, continuity requirements, and the domain of definition. For example, use of the complex JWKB method allows for the treatment of tunneling through the complexification of phase space. Finally, although KvH wavefunctions include the possibility of interference effects, interference is not observable when all observables are approximated as local operators on phase space. Observing interference effects requires consideration of nonlocal operations, e.g. through higher orders in the asymptotic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1997

Design of a Minimum Surface-Effect Tendon-Based Microactuator for Micromanipulation

A piezoelectric (PZT) stack-based actuator was developed to provide a means of actuation with dynamic characteristics appropriate for small-scale manipulation. In particular, the design incorporates a highly nonlinear, large-ratio transmission that provides approximately two orders of magnitude motion amplification from the PZT stack. In addition to motion amplification, the nonlinear transmission was designed via optimization methods to distort the highly non-uniform properties of a piezoelectric actuator so that the achievable actuation force is nearly constant throughout the actuator workspace. The package also includes sensors that independently measure actuator output force and displacement, so that a manipulator structure need not incorporate sensors nor the associated wires. Specifically, the actuator was designed to output a maximum force of at least one Newton through a stroke of at least one millimeter. For purposes of small-scale precision position and/or force control, the actuator/sensor package was designed to eliminate stick-slip friction and backlash. The overall dimensions of the actuator/sensor package are approximately 40 x 65 x 25 mm.

Goldfarb, Michael↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Accelerating an iterative process by explicit annihilation

A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.

Jespersen, D. C.↗
NASA NTRS1985

Accelerating an iterative process by explicit annihilation

A slowly convergent stationary iterative process can be accelerated by explicitly annihilating (i.e., eliminating) the dominant eigenvector component of the error. The dominant eigenvalue or complex pair of eigenvalues can be estimated from the solution during the iteration. The corresponding eigenvector or complex pair of eigenvectors can then be annihilated by applying an explicit Richardson process over the basic iterative method. This can be done entirely in real arithmetic by analytically combining the complex conjugate annihilation steps. The technique is applied to an implicit algorithm for the calculation of two dimensional steady transonic flow over a circular cylinder using the equations of compressible inviscid gas dynamics. This demonstrates the use of explicit annihilation on a nonlinear problem.

Jespersen, D. C.↗
NASA NTRS1990

Flight test of a trajectory controller using linearizing transformations with measurement feedback

The design of nonlinear controllers has relied on the use of detailed aerodynamic and engine models that must be associated with the control law in the flight system implementation. Many of these controllers have been applied to vehicle flightpath control problems and have attempted to combine both inner- and outer-loop control functions in a single controller. This paper presents an alternate approach to the design of outer-loop controllers. The approach simplifies the outer-loop design problem by separating the inner-loop (stabilization and control) from the outer-loop (guidance and navigation) functions. Linearizing transformations are applied using measurement feedback to eliminate the need for detailed aircraft models in outer-loop control applications. Also discussed is an implementation of the controller. This implementation was tested on a six-degree-of-freedom F-15 simulation and in flight on an F-15 aircraft. Proof of the concept is provided by flight test data which is presented and discussed.

Antoniewicz, Robert F.↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Toward Malaria Risk Prediction in Afghanistan Using Remote Sensing

Malaria causes more than one million deaths every year worldwide, with most of the mortality in Sub-Saharan Africa. It is also a significant public health concern in Afghanistan, with approximately 60% of the population, or nearly 14 million people, living in a malaria-endemic area. Malaria transmission has been shown to be dependent on a number of environmental and meteorological variables. For countries in the tropics and the subtropics, rainfall is normally the most important variable, except for regions with high altitude where temperature may also be important. Afghanistan s diverse landscape contributes to the heterogeneous malaria distribution. Understanding the environmental effects on malaria transmission is essential to the effective control of malaria in Afghanistan. Provincial malaria data gathered by Health Posts in 23 provinces during 2004-2007 are used in this study. Remotely sensed geophysical parameters, including precipitation from TRMM, and surface temperature and vegetation index from MODIS are used to derive the empirical relationship between malaria cases and these geophysical parameters. Both neural network methods and regression analyses are used to examine the environmental dependency of malaria transmission. And the trained models are used for predicting future transmission. While neural network methods are intrinsically more adaptive for nonlinear relationship, the regression approach lends itself in providing statistical significance measures. Our results indicate that NDVI is the strongest predictor. This reflects the role of irrigation, instead of precipitation, in Afghanistan for agricultural production. The second strongest prediction is surface temperature. Precipitation is not shown as a significant predictor, contrary to other malarious countries in the tropics or subtropics. With the regression approach, the malaria time series are modelled well, with average R2 of 0.845. For cumulative 6-month prediction of malaria cases, the average provincial accuracy reaches 91%. The developed predictive and early warning capabilities support the Third Strategic Approach of the WHO EMRO Malaria Control and Elimination Plan.

Safi, N.↗
DOE OSTI2024

Parasitic conduction loss of lithium niobate on insulator platform

Thin-film lithium niobate on insulator (LNOI) has recently emerged as an exciting platform for nonlinear integrated photonic devices, drawing parallels with the impact of silicon on insulator (SOI) in the electronic device realm. While the influence of interfacial states of SOI on electronic device performance, particularly their contribution to parasitic radio frequency (RF) losses, is well-documented, this aspect remains relatively unexplored in LNOI, despite their shared silicon/oxide interface structure. Here, we investigate the RF loss arising from mobile charges at the silicon/oxide interface of the LNOI platform. We show that this parasitic conduction can be exacerbated, especially during plasma etching, and leads to substantial RF losses. Through strategic incorporation of a trap-rich poly-silicon layer at this interface, we effectively eliminate this parasitic conduction, resulting in the suppression of RF losses. This work establishes crucial guidelines for harnessing LNOI in applications such as high-frequency electro-mechanical devices, high-speed electro-optic modulators, and high-efficiency microwave-to-optics converters.

Physics↗
NASA NTRS1991

The effects of pressure sensor acoustics on airdata derived from a High-angle-of-attack Flush Airdata Sensing (HI-FADS) system

The accuracy of a nonintrusive high angle-of-attack flush airdata sensing (HI-FADS) system was verified for quasi-steady flight conditions up to 55 deg angle of attack during the F-18 High Alpha Research Vehicle (HARV) Program. The system is a matrix of nine pressure ports arranged in annular rings on the aircraft nose. The complete airdata set is estimated using nonlinear regression. Satisfactory frequency response was verified to the system Nyquist frequency (12.5 Hz). The effects of acoustical distortions within the individual pressure sensors of the nonintrusive pressure matrix on overall system performance are addressed. To quantify these effects, a frequency-response model describing the dynamics of acoustical distortion is developed and simple design criteria are derived. The model adjusts measured HI-FADS pressure data for the acoustical distortion and quantifies the effects of internal sensor geometries on system performance. Analysis results indicate that sensor frequency response characteristics very greatly with altitude, thus it is difficult to select satisfactory sensor geometry for all altitudes. The solution used presample filtering to eliminate resonance effects, and short pneumatic tubing sections to reduce lag effects. Without presample signal conditioning the system designer must use the pneumatic transmission line to attenuate the resonances and accept the resulting altitude variability.

Whitmore, Stephen A.↗
NASA NTRS1996

Optimization of Air-Breathing Engine Concept

The design optimization of air-breathing propulsion engine concepts has been accomplished by soft-coupling the NASA Engine Performance Program (NEPP) analyzer with the NASA Lewis multidisciplinary optimization tool COMETBOARDS. Engine problems, with their associated design variables and constraints, were cast as nonlinear optimization problems with thrust as the merit function. Because of the large number of mission points in the flight envelope, the diversity of constraint types, and the overall distortion of the design space; the most reliable optimization algorithm available in COMETBOARDS, when used by itself, could not produce satisfactory, feasible, optimum solutions. However, COMETBOARDS' unique features-which include a cascade strategy, variable and constraint formulations, and scaling devised especially for difficult multidisciplinary applications-successfully optimized the performance of subsonic and supersonic engine concepts. Even when started from different design points, the combined COMETBOARDS and NEPP results converged to the same global optimum solution. This reliable and robust design tool eliminates manual intervention in the design of air-breathing propulsion engines and eases the cycle analysis procedures. It is also much easier to use than other codes, which is an added benefit. This paper describes COMETBOARDS and its cascade strategy and illustrates the capabilities of the combined design tool through the optimization of a high-bypass- turbofan wave-rotor-topped subsonic engine and a mixed-flow-turbofan supersonic engine.

Patnaik, Surya N.↗
NASA NTRS1976

Analytical and experimental study of structurally efficient composite hat-stiffened panels loaded in axial compression

Structural efficiency studies were made to determine the weight saving potential of graphite/epoxy composite structures for compression panel applications. Minimum weight hat-stiffened and open corrugation configurations were synthesized using a nonlinear mathematical programming technique. Selected configurations were built and tested to study local and Euler buckling characteristics. Test results for 23 panels critical in local buckling and six panels critical in Euler buckling are compared with analytical results obtained using the BUCLASP-2 branched plate buckling program. A weight efficiency comparison is made between composite and aluminum compression panels using metal test data generated by the NACA. Theoretical studies indicate that potential weight savings of up to 50% are possible for composite hat-stiffened panels when compared with similar aluminum designs. Weight savings of 32% to 42% were experimentally achieved. Experience suggests that most of the theoretical weight saving potential is available if design deficiencies are eliminated and strict fabrication control is exercised.

Williams, J. G.↗