SEARCH · Search NASA
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.
A particle model of the magnetotail neutral sheet self-consistent structure was developed by assuming that the thickness of the neutral sheet is much less than the thickness of the plasma sheet. This made it possible to approximate the neutral sheet as an x-independent structure. The role of nonlinear particle motions in determining the conductivity and the global current sheet structure is investigated.
We present Stokes I Zeeman splitting measurements of sunspots using the highly sensitive (g = 3) Fe I line at a wavelength of 1.5649 micron. The splittings are compared with simultaneous intensity measurements in the adjacent continuum. The relation between magnetic field strength and temperature has a characteristic, nonlinear shape in all the spots studied. In the umbra, there is an approximately linear relation between B super 2 and T sub b, consistent with magnetohydrostatic equilibrium in a nearly vertical field. A distinct flattening of the B super 2 vs T sub b relationship in the inner penumbra may be due to changes in the lateral pressure balance as the magnetic field becomes more horizontal; spatially unresolved intensity inhomogeneities may also influence the observed relation.
Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfven speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis on the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere, and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most- commonly -quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating TRACE satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.
Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfven speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis an the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most commonly quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating Transition Region and Coronal Explorer (TRACE) satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.
Various color histogram equalization (CHE) methods have been proposed to extend grayscale histogram equalization (GHE) for color images. In this paper a new method called histogram diffusion that extends the GHE method to arbitrary dimensions is proposed. Ranges in a histogram are specified as overlapping bars of uniform heights and variable widths which are proportional to their frequencies. This diagram is called the vistogram. As an alternative approach to GHE, the squared error of the vistogram from the uniform distribution is minimized. Each bar in the vistogram is approximated by a Gaussian function. Gaussian particles in the vistoram diffuse as a nonlinear autonomous system of ordinary differential equations. CHE results of color images showed that the approach is effective.
The second and third-order Brugger elastic constants are obtained for liquids and ideal gases having an initial hydrostatic pressure p(sub 1). For liquids the second-order elastic constants are C(sub 11) = A + p(sub 1), C(sub 12) = A -- p(sub 1), and the third-order constants are C(sub 111) = --(B + 5A + 3p(sub 1)), C(sub 112) = --(B + A -- p(sub 1)), and C(sub 123) = A -- B -- p1, where A and B are the Beyer expansion coefficients in the liquid equation of state. For ideal gases the second order constants are C(sub 11) = p(sub 1)gamma + p9sub 1), C(sub 12) = p(sub 1)gamma -- p(sub 1), and the third-order constants are C(sub 111) = p(sub 1)(gamma(2) + 4gamma + 3), C(sub 112) = --p(sub 1)(gamma(2) -- 1), and C(sub 123) = --p(sub 1) (gamma(2) -- 2gamma + 1), where gamma is the ratio of specific heats. The inequality of C(sub 11) and C(sub 12) results in a nonzero shear constant C(sub 44) = (1/2)(C(sub 11) C(sub 12)) = p(sub 1) for both liquids and gases. For water at standard temperature and pressure the ratio of terms p1/A contributing to the second-order constants is approximately 4.3 x 10(-5). For atmospheric gases the ratio of corresponding terms is approximately 0.7. Analytical expressions that include initial stresses are derived for the material 'nonlinearity parameters' associated with harmonic generation and acoustoelasticity for fluids and solids of arbitrary crystal symmetry. The expressions are used to validate the relationships for the elastic constants of fluids.
This paper covers the application of a model-based engine control (MBEC) methodology featuring a self tuning on-board model for an aircraft turbofan engine simulation. The nonlinear engine model is capable of modeling realistic engine performance, allowing for a verification of the advanced control methodology over a wide range of operating points and life cycle conditions. The on-board model is a piece-wise linear model derived from the nonlinear engine model and updated using an optimal tuner Kalman Filter estimation routine, which enables the on-board model to self-tune to account for engine performance variations. MBEC is used here to show how advanced control architectures can improve efficiency during the design phase of a turbofan engine by reducing conservative operability margins. The operability margins that can be reduced, such as stall margin, can expand the engine design space and offer potential for efficiency improvements. Application of MBEC architecture to a nonlinear engine simulation is shown to reduce the thrust specific fuel consumption by approximately 1% over the baseline design, while maintaining safe operation of the engine across the flight envelope.
A mathematical procedure is presented for the calculation of exact cumulative distribution statistics for a viscosity-free variant of Burgers nonlinear partial differential equation (PDE) in one space dimension and time subject to sinusoidal initial data with uncertain (random variable) amplitude or phase shift. Analytical solutions of nonlinear PDEs with uncertain initial and/or boundary data are invaluable benchmarks in assessing approximate uncertainty quantification techniques. The Burgers equation solution with uncertain initial data results in nonsmooth solution behavior in both physical and random variable dimensions which provides a severe test for approximate uncertainty quantification techniques. Mathematical proofs are provided to verify that exact cumulative distribution statistics can be systematically and robustly obtained for all forward time.
Accurate and quick identification of high-impedance faults (HIFs) is critical for the reliable operation of distribution systems. Unlike other faults in power grids, HIFs are very difficult to detect by conventional overcurrent relays due to the low fault current. Although HIFs can be affected by various factors, the voltage-current characteristics can substantially imply how the system responds to the disturbance and thus provides opportunities to effectively localize HIFs. In this work, we propose a data-driven approach for the identification of HIF events. To tackle the nonlinearity of the voltage-current trajectory, first, we formulate optimization problems to approximate the trajectory with piecewise functions. Then we collect the function features of all segments as inputs and use the support vector machine approach to efficiently identify HIFs at different locations. Numerical studies on the IEEE 123-node test feeder demonstrate the validity and accuracy of the proposed approach for real-time HIF identification.
Experiments were performed to determine the plasma propagation velocities resulting from KrF laser irradiation of copper wire target (75 microns diameter) and water droplets (75 microns diameter) at irradiance levels ranging from 25 to 150 GW/sq cm. Plasma propagation velocities were measured using a streak camera system oriented orthogonally to the high-energy laser propagation axis. Plasma velocities were studied as a function of position in the focused beam. Results show that both the shape of the plasma formation and material removal from the copper wire are different and depend on whether the targets are focused or slightly defocused (approximately = 0.5 mm movement in the beam axis). Plasma formation and its position relative to the target is an important factor in determining the practical focal point during high-energy laser interaction with materials. At irradiance of 100 GW/sq cm, the air plasma has two weak-velocity components which propagate toward and away from the incident laser while a strong-velocity component propagates away from the laser beam as a detonation wave. Comparison of the measured breakdown velocities (in the range of 2.22-2.27 x 10(exp 5) m/s) for air and the value calculated by the nonlinear breakdown wave theory at irradiance of 100 GW/sq cm showed a quantitative agreement within approximately 50% while the linear theory and Gaussian pulse theory failed. The detonation wave velocities of plasma generated from water droplets and copper wire targets for different focused cases were measured and analyzed theoretically. The propagation velocities of laser-induced plasma liquid droplets obtained by previous research are compared with current work.
An optical heterodyne interferometer that can be used to measure linear displacements with an error <=20 pm has been developed. The remarkable accuracy of this interferometer is achieved through a design that includes (1) a wavefront split that reduces (relative to amplitude splits used in other interferometers) self interference and (2) a common-optical-path configuration that affords common-mode cancellation of the interference effects of thermal-expansion changes in optical-path lengths. The most popular method of displacement- measuring interferometry involves two beams, the polarizations of which are meant to be kept orthogonal upstream of the final interference location, where the difference between the phases of the two beams is measured. Polarization leakages (deviations from the desired perfect orthogonality) contaminate the phase measurement with periodic nonlinear errors. In commercial interferometers, these phase-measurement errors result in displacement errors in the approximate range of 1 to 10 nm. Moreover, because prior interferometers lack compensation for thermal-expansion changes in optical-path lengths, they are subject to additional displacement errors characterized by a temperature sensitivity of about 100 nm/K. Because the present interferometer does not utilize polarization in the separation and combination of the two interfering beams and because of the common-mode cancellation of thermal-expansion effects, the periodic nonlinear errors and the sensitivity to temperature changes are much smaller than in other interferometers
The development of relaxation methods to calculate transonic flows is discussed. Rather accurate predictions can be made for a number of flows of interest using the transonic potential flow equation, which may be derived from the Euler equations for inviscid compressible flow by introducing the assumption that the flow is irrotational. Once the choice of a mathematical model has been settled, the numerical procedure for actually computing a solution contains two main elements: the construction of a discrete approximation which converges to the solution of the continuous problem in the limit as the mesh width is reduced to zero, and the solution of the resulting set of nonlinear difference equations by a convergent iterative scheme. The choice of an appropriate coordinate system and its influence on the accuracy of the discrete approximation is discussed. Several applications of the general method are described.
A finite element study is presented of double-diffusive convection driven by g-jitter in a microgravity environment. Mathematical formulations are presented and extensive simulations are carried out for g-jitter induced fluid flow, temperature distribution, and solutal transport in an alloy system under consideration for space flights. Computations include the use of idealized single-frequency and multi-frequency g-jitter as well as the real g-jitter data taken during an actual Space Shuttle fight. Little correlation is seen between these velocity components for the g-jitter components studied. The temperature field is basically undisturbed by convection because of a small Pr number for the fluid. The disturbance of the concentration field, however, is pronounced, and the local variation of the concentration follows the velocity oscillation in time. It is found that although the concentration field varies in both position and time, the local concentration gradient remains approximately constant in time. Numerical study further indicates that with an increase in g-jitter force (or amplitude), the nonlinear convective effects become much more obvious, which in turn drastically change the concentration fields. The simulated results computed using the g-jitter data taken during space flights show that both the velocity and concentration become random, following approximately the same pattern as the g-jitter perturbations.
Several laws are designed for simultaneous control of the orientation of an Earth-pointing spacecraft, the energy stored by counter-rotating flywheels, and the angular momentum of the flywheels and control moment gyroscopes used together as an integrated set of actuators for attitude control. General, nonlinear equations of motion are presented in vector-dyadic form, and used to obtain approximate expressions which are then linearized in preparation for design of control laws that include feedback of flywheel kinetic energy error as a means of compensating for damping exerted by rotor bearings. Two flywheel steering laws are developed such that torque commanded by an attitude control law is achieved while energy is stored or discharged at the required rate. Using the International Space Station as an example, numerical simulations are performed to demonstrate control about a torque equilibrium attitude, and illustrate the benefits of kinetic energy error feedback. Control laws for attitude hold are also developed, and used to show the amount of propellant that can be saved when flywheels assist the CMGs. Nonlinear control laws for large-angle slew maneuvers perform well, but excessive momentum is required to reorient a vehicle like the International Space Station.
Implicit, noniterative, finite difference schemes were recently developed by several authors for multidimensional systems of nonlinear hyperbolic partial differential equations. When applied to linear model equations with periodic boundary conditions those schemes are unconditionally stable (A-stable). As applied in practice the algorithms often face a severe time step restriction. A major source of the difficulty is the treatment of the numerical boundary conditions. One conjecture was that unconditional stability requires implicit numerical boundary conditions. An apparent counter example was the space time extrapolation considered by Gustafsson, Kreiss, and Sunstrom. Spatial (implicit) and space time (explicit) extrapolation using normal mode analysis for a finite and infinite number of spatial mesh intervals are examined. The results indicate that for unconditional stability with a finite number of spatial mesh intervals, the numerical boundary conditions must be implicit.
We approach the problem of parallelizing state-of-the-art nonlinear programming optimization algorithms.Specifically, we focus on parallelizing quasi-Newton interior-point methods that use limited-memory secant Hessian approximations. Such interior-point methods are known to have better convergence properties and to be more effective on large-scale problems than gradient-based and derivative-free optimization algorithms. We target nonlinear and potentially nonconvex optimization problems with an arbitrary number of bound constraints and a small number of general equality and inequality constraints on the optimization variables. These problems occur for example in the form of optimal control, optimal design, and inverse problems governed by ordinary or partial differential equations, whenever they are expressed in a “reduced-space” optimization approach. We introduce and analyze the time and space complexity of a decomposition method for solving the quasi-Newton linear systems that leverages the fact that the quasi-Newton Hessian matrix has a small number of dense blocks that border a low-rank update of a diagonal matrix. This enables an efficient parallelization on memory-distributed computers of the iterations of the optimization algorithm, a state-of-the-art filter line-search interior-point algorithm by Wächter et. al. We illustrate the efficiency of the proposed method by solving structural topology optimization problems on up to 4608 cores on a parallel machine.
Our study details the derivation of the nonlinear equations of motion for the axial, biaxial bending and torsional vibrations of an aeroelastic cantilever undergoing rigid body (pitch) rotation at the base. The primary attenstion is focussed on the geometric nonlinearities of the system, whereby the aeroelastic load is modeled by the theory of linear quasisteady aerodynamics. This modelling effort is intended to mimic the wind-tunnel experimental setup at the Royal Military College of Canada. While the derivation closely follows the work of Hodges and Dowell [1] for rotor blades, this aeroelastic system contains new inertial terms which stem from the fundamentally different kinematics than those exhibited by helicopter or wind turbine blades. Using the Hamilton’s principle, a set of coupled nonlinear partial differential equations (PDEs) and an ordinary differential equation (ODE) are derived which describes the coupled axial-bending-bending-torsion-pitch motion of the aeroelastic cantilever with the pitch rotation. The finite dimensional approximation of the coupled system of PDEs are obtained using the Galerkin projection, leading to a coupled system of ODEs. Subsequently, these nonlinear ODEs are solved numerically using the built-in MATLAB implicit ODE solver and the associated numerical results are compared with those obtained using Houbolt’s method. It is demonstrated that the system undergoes coalescence flutter, leading to a limit cycle oscillation (LCO) due to coupling between the rigid body pitching mode and teh flexible mode arising from the flapwise bending motion.
In this work, the multimodal nonlinear elastic behavior of concrete, which is representative of a consolidated granular material, is modeled numerically. Starting from a local three-dimensional softening law, the initial stiffness properties are re-estimated according to the local strain field. The experiments deal with samples of thermally damaged concrete blocks successively excited around their first three modes of vibration. The geometry of these samples cannot be described by a one-dimensional approximation in these experiments where compressional and shear motions are strongly coupled. Despite this added complexity, the nonlinear behavior for the three modes of vibration of the samples is well captured by the simulations using a single scalar nonlinear parameter appropriately integrated into the elasticity equations. Finally, it is shown that without sufficient attention paid to the latter, the conclusions would have brought erroneous statements such as nonlinearity dispersion or strain type dependence.