Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

Here, in this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The influence of Al2O3 on the structural properties of MgSiO3 akimotoite

Abstract Akimotoite, a MgSiO3 polymorph present in the lower transition zone within ultramafic portions of subducting slabs and potentially also in the ambient mantle, will partition some amount of Al, raising the question of how this will affect its crystal structure and properties. In this study, a series of samples along the MgSiO3-Al2O3 (akimotoite-corundum) solid solution have been investigated by means of single-crystal X-ray diffraction to examine their crystal chemistry. Results show a strong nonlinear behavior of the a- and c-axes as a function of Al content, which arises from fundamentally different accommodation mechanisms in the akimotoite and corundum structures. Furthermore, two Al2O3-bearing akimotoite samples were investigated at high pressure to determine the different compression mechanisms associated with Al substitution. Al2O3-bearing akimotoite becomes more compressible at least up to 20 mol% Al2O3, due likely to an increase in compressibility as the Al cation is incorporated into the SiO6 octahedron. This observation is in strong contrast to the stiffer corundum end-member having a KT = 250 GPa, which is larger than that of the akimotoite end-member [KT = 205(1) GPa]. These findings have implications for mineral physics models of elastic properties, which have in the past assumed linear mixing behavior between the MgSiO3 akimotoite and Al2O3 corundum end-members to calculate sound wave velocities for Al-bearing akimotoite at high pressure and temperature.

Geochemistry & Geophysics↗
DOE OSTI2021

HIRM v1.0: a hybrid impulse response model for climate modeling and uncertainty analyses

Simple climate models (SCMs) are frequently used in research and decision-making communities because of their flexibility, tractability, and low computational cost. SCMs can be idealized, flexibly representing major climate dynamics as impulse response functions, or process-based, using explicit equations to model possibly nonlinear climate and Earth system dynamics. Each of these approaches has strengths and limitations. Here we present and test a hybrid impulse response modeling framework (HIRM) that combines the strengths of process-based SCMs in an idealized impulse response model, with HIRM's input derived from the output of a process-based model. This structure enables the model to capture some of the major nonlinear dynamics that occur in complex climate models as greenhouse gas emissions transform to atmospheric concentration to radiative forcing to climate change. As a test, the HIRM framework was configured to emulate the total temperature of the simple climate model Hector 2.0 under the four Representative Concentration Pathways and the temperature response of an abrupt 4 times CO 2 concentration step. HIRM was able to reproduce near-term and long-term Hector global temperature with a high degree of fidelity. Additionally, we conducted two case studies to demonstrate potential applications for this hybrid model: examining the effect of aerosol forcing uncertainty on global temperature and incorporating more process-based representations of black carbon into a SCM. The open-source HIRM framework has a range of applications including complex climate model emulation, uncertainty analyses of radiative forcing, attribution studies, and climate model development.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1971

An optimal system design process for a Mars roving vehicle

The problem of determining the optimal design for a Mars roving vehicle is considered. A system model is generated by consideration of the physical constraints on the design parameters and the requirement that the system be deliverable to the Mars surface. An expression which evaluates system performance relative to mission goals as a function of the design parameters only is developed. The use of nonlinear programming techniques to optimize the design is proposed and an example considering only two of the vehicle subsystems is formulated and solved.

Pavarini, C.↗
NASA NTRS1973

Quasi-static extension of a tensile crack contained in a viscoelastic-plastic solid

Final stretch criterion of failure is applied to the problem of quasi-static extension of a crack embedded in an elastic-plastic or viscoelastic-plastic matrix. The slow growth under subcritical conditions in a rate-sensitive Tresca solid is shown to be a superposition of creep rupture and McClintock's ductile growth. This type of growth occurs at subcritical magnitude of the imposed K-factor and can be accounted for only through a recognition of inelastic properties of solids. In the subcritical range there is no unique value for K sub c independent of geometrical configuration and flaw size. Not only the produced states of stress and strain are dependent on the loading path, but also the material resistance to fracture turns out to be a function of the history of loading that precedes catastrophic failure. A nonlinear integro-differential equation of motion is derived for a crack progressing through a viscoelastic medium with some limited ability to plastic flow. Examples of numerical integration are given incorporating both monotonic and cyclic loading programs.

Wnuk, M. P.↗
NASA NTRS1973

An all digital phase locked loop for synchronization of a sinusoidal signal embedded in white Gaussian noise

An all digital phase locked loop which tracks the phase of the incoming sinusoidal signal once per carrier cycle is proposed. The different elements and their functions and the phase lock operation are explained in detail. The nonlinear difference equations which govern the operation of the digital loop when the incoming signal is embedded in white Gaussian noise are derived, and a suitable model is specified. The performance of the digital loop is considered for the synchronization of a sinusoidal signal. For this, the noise term is suitably modelled which allows specification of the output probabilities for the two level quantizer in the loop at any given phase error. The loop filter considered increases the probability of proper phase correction. The phase error states in modulo two-pi forms a finite state Markov chain which enables the calculation of steady state probabilities, RMS phase error, transient response and mean time for cycle skipping.

Reddy, C. P.↗
NASA NTRS1974

Nonlinear initial-boundary value solutions by the finite element method

The finite-element method is used to establish a numerical solution algorithm for the Navier-Stokes equations for two-dimensional flows of a viscous compressible fluid. Numerical experiments confirm the advection property for the finite-element equivalent of the nonlinear convection term for both unidirectional and recirculating flowfields. For linear functionals, the algorithm demonstrates good accuracy using coarse discretizations and h squared convergence with discretization refinement.

Baker, A. J.↗
NASA NTRS1976

Polyatomic molecule vibrations

Polyatomic molecule vibrations are analyzed as harmonic vibrations along normal coordinates. The energy eigenvalues are found for linear and nonlinear symmetric triatomic molecules for valence bond models of the potential function with arbitrary coupling coefficients; such models can usually be fitted to observed energy levels with reasonably good accuracy. Approximate normal coordinates for the H2O molecule are discussed. Degenerate vibrational modes such as occur in CO2 are analyzed and expressions for Fermi resonance between close-lying states of the same symmetry are developed. The bending modes of linear triatomic molecules are expressed in terms of Laguerre polynomials in cylindrical coordinates as well as in terms of Hermite polynomials in Cartesian coordinates. The effects of large-amplitude bending such as occur in the C3 molecule are analyzed, along with anharmonic effects, which split the usually degenerate bending mode energy levels. Finally, the vibrational frequencies, degeneracies, and symmetry properties of XY3, X2Y2, and XY4 type molecules are discussed.

Source record↗
NASA NTRS1977

Impact cratering and regolith dynamics

The most recent models concerning mechanical aspects of lunar regolith dynamics related to impact cratering use probabilistic approaches to account for the randomness of the meteorite environment in both space and time. Accordingly the absolute regolith thickness is strictly a function of total bombardment intensity and absolute regolith growth rate in nonlinear through geologic time. Regoliths of increasing median thickness will have larger and larger proportions of more and more deep seated materials. An especially active zone of reworking on the lunar surface of about 1 mm depth has been established. With increasing depth, the probability of excavation and regolith turnover decreases very rapidly. Thus small scale stratigraphy - observable in lunar core materials - is perfectly compatible with regolith gardening, though it is also demonstrated that any such stratigraphy does not necessarily present a complete record of the regolith's depositional history. At present, the lifetimes of exposed lunar rocks against comminution by impact processes can be modeled; it appears that catastrophic rupture dominates over single particle abrasion.

Hoerz, F.↗
NASA NTRS1981

Aerodynamic mathematical modeling - basic concepts

This lecture is a review of the basic concepts involved In the mathematical modeling of the aerodynamic response of an aircraft to arbitrary maneuvers. The review leads from Bryan's original formulation through the concepts of linear aerodynamic indicial functions and superposition to the extension of these concepts into the nonlinear regime. The nonlinear generalization yields a form for the aerodynamic response that can be built up from the responses to a limited number of well-defined characteristic motions, reproducible in principle either in wind-tunnel experiments or flow-field computations. A further generalization leads to a torn, accommodating the discontinuous and double-valued behavior characteristic of hysteresis in the steady-state aerodynamic response

Murray Tobak↗
NASA NTRS1988

Application of fuzzy theories to formulation of multi-objective design problems

Much of the decision making in real world takes place in an environment in which the goals, the constraints, and the consequences of possible actions are not known precisely. In order to deal with imprecision quantitatively, the tools of fuzzy set theory can by used. This paper demonstrates the effectiveness of fuzzy theories in the formulation and solution of two types of helicopter design problems involving multiple objectives. The first problem deals with the determination of optimal flight parameters to accomplish a specified mission in the presence of three competing objectives. The second problem addresses the optimal design of the main rotor of a helicopter involving eight objective functions. A method of solving these multi-objective problems using nonlinear programming techniques is presented. Results obtained using fuzzy formulation are compared with those obtained using crisp optimization techniques. The outlined procedures are expected to be useful in situations where doubt arises about the exactness of permissible values, degree of credibility, and correctness of statements and judgements.

Dhingra, A. K.↗
NASA NTRS1993

Transition to turbulence in an elliptic vortex

We study the three dimensional instability and nonlinear growth of the two dimensional flow described by the stream function Psi = (A sin b1 x sin b2 y)/(b1(exp 2) + b2(exp 2)) where b1 = pi/L1, b2 = pi/L2. This is a swirling flow in a box which is bounded by 0 less than x less than L1, 0 less than y less than L2 and is infinite in the z direction. This flow is a solution of the Navier-Stokes equation with A = exp(-v(b1(exp 2) + b2(exp 2))t) which slowly decays. We seek a viscous solution which starts near this one and slips along but does not penetrate the bounding walls. The vorticity of the basic flow is w(sub z) = A sin b1 x sin b2 y which has maximum value A at the center of the box and drops to zero at the boundaries. We can think of the resulting flow as that of a captive vortex.

Lundgren, Thomas S.↗
NASA NTRS1996

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗
NASA NTRS2003

Flaw Tolerance In Lap Shear Brazed Joints

This paper presents results of the second part of an on-going effort to gain better understanding of defect tolerance in braze joints. In the first part of this three-part series, we mechanically tested and modeled the strength of the lap joints as a function of the overlap distance. A failure criterion was established based on the zone damage theory, which predicts the dependence of the lap joint shear strength on the overlap distance, based on the critical size of a finite damage zone or an overloaded region in the joint. In this second part of the study, we experimentally verified the applicability of the damage zone criterion on prediction of the shear strength of the lap joint and introduced controlled flaws into the lap joints. The purpose of the study was to evaluate the lap joint strength as a function of flaw size and its location through mechanical testing and nonlinear finite element analysis (FEA) employing damage zone criterion for definition of failure. The results obtained from the second part of the investigation confirmed that the failure of the ductile lap shear brazed joints occurs when the damage zone reaches approximately 10% of the overlap width. The same failure criterion was applicable to the lap joints containing flaws.

Wang, Len↗
NASA NTRS2008

Acoustic Rectification in Dispersive Media

It is shown that the shapes of acoustic radiation-induced static strain and displacement pulses (rectified acoustic pulses) are defined locally by the energy density of the generating waveform. Dispersive properties are introduced analytically by assuming that the rectified pulses are functionally dependent on a phase factor that includes both dispersive and nonlinear terms. The dispersion causes an evolutionary change in the shape of the energy density profile that leads to the generation of solitons experimentally observed in fused silica.

Cantrell, John H.↗
DOE OSTI2022

McMillan map and nonlinear Twiss parameters

In this article we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings originally introduced by Edwin McMillan; the second one is also known as canonical McMillan map. Both of them are simplest symmetric McMillan maps with only one intrinsic parameter, the trace of the Jacobian at the fixed point. While these dynamical systems have numerous of applications and are used in many areas of math and physics, some of their dynamical properties have not been described yet. We fulfill the gap and provide complete description of all stable trajectories including parametrization of invariant curves, Pioncaré rotation numbers and canonical action-angle variables. In the second part we relate these maps with general chaotic map in McMillan-Turaev form. We show that McMillan sextupole and octupole mappings are first order approximations of dynamics around the fixed point, in a similar way as linear map and quadratic invariant (Courant-Snyder invariant in accelerator physics) is the zeroth order approximation (known as linearization). Finally we suggest the new formalism of nonlinear Twiss parameters which incorporate dependence of rotation number as a function of amplitude, in contrast to e.g. betatron phase advance used in accelerator physics which is independent of amplitude. Specifically in application to accelerator physics this new formalism is capable of predicting dynamical aperture around 1-st, 2-nd, 3-rd and 4-th order resonances for flat beams, which is critical for beam injection/extraction.

43 PARTICLE ACCELERATORS↗