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 289 records · Page 16

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 programing 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 to date 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.↗

The NCOREL computer program for 3D nonlinear supersonic potential flow computations

An innovative computational technique (NCOREL) was established for the treatment of three dimensional supersonic flows. The method is nonlinear in that it solves the nonconservative finite difference analog of the full potential equation and can predict the formation of supercritical cross flow regions, embedded and bow shocks. The method implicitly computes a conical flow at the apex (R = 0) of a spherical coordinate system and uses a fully implicit marching technique to obtain three dimensional cross flow solutions. This implies that the radial Mach number must remain supersonic. The cross flow solutions are obtained by using type dependent transonic relaxation techniques with the type dependency linked to the character of the cross flow velocity (i.e., subsonic/supersonic). The spherical coordinate system and marching on spherical surfaces is ideally suited to the computation of wing flows at low supersonic Mach numbers due to the elimination of the subsonic axial Mach number problems that exist in other marching codes that utilize Cartesian transverse marching planes.

Siclari, M. J.↗

Verification and Validation Challenges for Adaptive Flight Control of Complex Autonomous Systems

Autonomy of aerospace systems requires the ability for flight control systems to be able to adapt to complex uncertain dynamic environment. In spite of the five decades of research in adaptive control, the fact still remains that currently no adaptive control system has ever been deployed on any safety-critical or human-rated production systems such as passenger transport aircraft. The problem lies in the difficulty with the certification of adaptive control systems since existing certification methods cannot readily be used for nonlinear adaptive control systems. Research to address the notion of metrics for adaptive control began to appear in the recent years. These metrics, if accepted, could pave a path towards certification that would potentially lead to the adoption of adaptive control as a future control technology for safety-critical and human-rated production systems. Development of certifiable adaptive control systems represents a major challenge to overcome. Adaptive control systems with learning algorithms will never become part of the future unless it can be proven that they are highly safe and reliable. Rigorous methods for adaptive control software verification and validation must therefore be developed to ensure that adaptive control system software failures will not occur, to verify that the adaptive control system functions as required, to eliminate unintended functionality, and to demonstrate that certification requirements imposed by regulatory bodies such as the Federal Aviation Administration (FAA) can be satisfied. This presentation will discuss some of the technical issues with adaptive flight control and related V&V challenges.

Nguyen, Nhan T.↗

Natural convection in steady solidification - Finite element analysis of a two-phase Rayleigh-Benard problem

Galerkin finite-element approximations and Newton's method for solving free boundary problems are combined with computer-implemented techniques from nonlinear perturbation analysis to study solidification problems with natural convection in the melt. The Newton method gives rapid convergence to steady state velocity, temperature and pressure fields and melt-solid interface shapes, and forms the basis for algebraic methods for detecting multiple steady flows and assessing their stability. The power of this combination is demonstrated for a two-phase Rayleigh-Benard problem composed of melt and solid in a veritical cylinder with the thermal boundary conditions arranged so that a static melt with a flat melt-solid interface is always a solution. Multiple cellular flows bifurcating from the static state are detected and followed as Rayleigh number is varied. Changing the boundary conditions to approach those appropriate for the vertical Bridgman solidification system causes imperfections that eliminate the static state. The flow structure in the Bridgman system is related to those for the Rayleigh-Benard system by a continuous evolution of the boundary conditions.

Chang, C. J.↗

A finite micro-rotation material point method for micropolar solid and fluid dynamics with three-dimensional evolving contacts and free surfaces

This paper introduces an explicit material point method designed specifically for simulating the micropolar continuum dynamics in the finite deformation and finite microrotation regime. The material point method enables us to simulate large deformation problems while circumventing the potential mesh distortion without remeshing. To eliminate rotational motion damping and loss of angular momentum during the projection, we introduce the mapping for microinertia and angular momentum between particles and grids through the affine particle-in-cell approach. The microrotation and the curvature at each particle are updated through zero-order forward integration of the microgyration and its spatial gradient. We show that the microinertia and the angular momentum are conserved during the projections between particles and grids in our formulation. We verify the formulation and implementation by comparing with the analytical dispersion relation of micropolar waves under the small strain and small microrotation, as well as the analytical soliton solution for solids undergoing large deformation and large microrotation. Additionally we also demonstrate the capacity of the proposed computational framework to handle a wide spectrum of simulations that exhibit size effects in the geometrical nonlinear regime through three representative numerical examples, i.e., a cantilever beam torsion problem, a fragment-impact penetration problem, and a micropolar fluid discharging problem.

42 ENGINEERING↗

An analysis of rotor blade twist variables associated with different Euler sequences and pretwist treatments

A nonlinear analysis which is necessary to adequately model elastic helicopter rotor blades experiencing moderately large deformations was examined. The analysis must be based on an appropriate description of the blade's deformation geometry including elastic bending and twist. Built-in pretwist angles complicate the deformation process ant its definition. Relationships between the twist variables associated with different rotation sequences and corresponding forms of the transformation matrix are lasted. Relationships between the twist variables associated with first, the pretwist combined with the deformation twist are included. Many of the corresponding forms of the transformation matrix for the two cases are listed. It is shown that twist variables connected with the combined twist treatment are related to those where the pretwist is applied initially. A method to determine the relationships and some results are outlined. A procedure to evaluate the transformation matrix that eliminates the Eulerlike sequence altogether is demonstrated. The resulting form of the transformation matrix is unaffected by rotation sequence or pretwist treatment.

Alkire, K.↗

Design and Analysis of High Heat Flux Plasma-Facing Components for NSTX Upgrade

The National Spherical Tokamak Experiment (NSTX) upgrade (NSTX-U) requirements lead to enhanced heat loads on plasma-facing components (PFCs) especially in the divertor regions, where normal heat flux density can reach up to 8 MW/m 2 . For these high heat flux (HHF) regions, the PFCs were redesigned, to use castellations which relieved the thermal stresses associated with high incident heat flux. Another design feature of HHF PFCs is the absence of front face mounting holes which create localized areas of high temperature and thermal stress concentrations. Optimized shaping of the front face of the HHF PFCs eliminates regions where the front face is perpendicular to the heat flux direction during normal helicity operation and spreads the heat load. A unique mechanism was designed to connect PFC tiles to the NSTX-U center stack casing using locking features accessible from the low heat flux regions. Isotropic graphite was selected as the HHF PFC tile material. Initial tile design parameters were accessed using analytical expressions for pulsed heat flux loading. A working prototype of the locking mechanism was created during the initial stages of the design to prove the concept performance. This article presents an overview of the divertor HHF PFC design and the results of the thermal and structural analyses performed using ANSYS software. The results of the analyses cover normal operating conditions and disruptions which impose electromagnetic (EM) loads from eddy and halo currents. The 3-D, transient, nonlinear analyses took into account the temperature-dependent properties of the materials, friction interfaces between the parts, and variable electric properties of the parts and interfaces. The results confirm that the HHF PFC tiles remain within the allowable limits for the loads defined by the NSTX-U Recovery Project. In addition, tolerance stack up analyses were performed to ensure tile performance in the worst possible assembly configuration. Finally, the design process was completed successfully, and the NSTX HHF PFC tiles are currently in the production phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Titanium Niobium Oxide: From Discovery to Application in Fast-Charging Lithium-Ion Batteries

Lithium-ion batteries are essential for portable technology and are now poised to disrupt a century of combustion-based transportation. The electrification revolution could eliminate our reliance on fossil fuels and enable a clean energy future; advanced batteries would facilitate this transition. However, owing to the demanding performance, cost, and safety requirements, it is challenging to translate new materials from laboratory prototypes to industrial-scale products. This Perspective describes that journey for a new lithium-ion battery anode material, TiNb 2 O 7 (TNO). TNO is intended as an alternative to graphite or Li 4 Ti 5 O 12 with better rate and safety characteristics than the former and higher energy density than the latter. The high capacity of TNO stems from the multielectron redox of Nb 5+ to Nb 3+ , its operating voltage window well above the Li + /Li reduction potential prevents lithium dendrite formation, and its open crystal structure leads to high-power performance. Nevertheless, the creation of a practical TNO anode was nonlinear and nontrivial. Its history is built on 30 years of fundamental science that preceded its application as a battery anode, and its battery development included a nearly 30-year gap. The insights and lessons contained in this Perspective, many of them acquired firsthand, serve two purposes: (i) to unite the disparate studies of TiNb 2 O 7 into a coherent modern understanding relevant to its application as a battery material and (ii) to highlight briefly some of the challenges faced when scaling up a new material that affect TiNb 2 O 7 as well as new electrode candidates more generally.

25 ENERGY STORAGE↗

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗

An evaluation of corrosion protection by two epoxy primers on 2219-T87 and 7075-T73 aluminum

A comparison of the corrosion protection provided by two amine epoxy primers was made using salt fog, alternate immersion, and total immersion as exposure media. The study is the result of a request to use an unqualified low volatile organic carbon (VOC) primer (AKZO 463-6-78) in place of the current primer (AKZO 463-6-3) because environmental regulations have eliminated use of the current primer in many states. Primed, scribed samples of 2219-T87 and 7075-T73 aluminum were exposed to 5-percent NaCl salt fog and 3.5-percent NaCl alternate immersion for a period of 90 days. In addition, electrode samples immersed in 3.5-percent NaCl were tested using electrochemical impedance spectroscopy (EIS). The EG&G model 368 ac impedance measurement system was used to monitor changing properties of AKZO 463-6-78 and AKZO 463-6-3 primed 2219-T87 aluminum for a period of 30 days. The response of the corroding system of a frequency scan can be modeled in terms of an equivalent circuit consisting of resistors and capacitors in a specific arrangement. Each resistor/capacitor combination represents physical processes taking place within the electrolyte, at the electrolyte/primer surface, within the coating, and at the coating/substrate surface. Values for the resistors and capacitors are assigned following a nonlinear least squares fit of the data to the equivalent circuit. Changes in the values of equivalent circuit parameters during the 30-day exposure allow assessment of the time to and mechanism of coating breakdown.

Mendrek, M. J.↗

Apollo 9 multiband photography experiment S065

Fourier analysis was applied to microdensitometer scans of a selected region of one SO65 frame in each of the three black-and-white bands. The approach was unique because a somewhat arbitrary section of the image was used and not limited to available targets as in edge analysis. Comparison of duplicates and calculation of absolute SO65 MTF were done by applying linear systems theory to the spatial spectra of the image scans. It was found that the duplication process was nonlinear and resulted in general amplification of spatial frequency modulation. However, the increase in modulation was offset by a corresponding increase in the granularity of the copies. The amount of increase seemed to be related to the initial granularity, but a direct relationship was not verified. Band-to-band comparison of image quality was achieved in the form of signal-to-noise ratio curves as a function of spatial frequency for each band. From this standpoint the DD band was an order of magnitude better than the other two. These were several factors that restricted the analysis of the Apollo 9 imagery. Among these were the lack of precise sensitometric and optical system data on the high altitude photography. In addition it was determined that there was only one simultaneous pair of high altitude and SO65 frames. Finally, the original SO65 photography was not available for scanning (for obvious reasons), thus eliminating a reference base for granularity and SO65 MTF determination.

Schowengerdt, R. A.↗

Investigating Thaw and Plant Productivity Constraints on Old Soil Carbon Respiration From Permafrost

Isotopic radiocarbon (Δ 14 C) signatures of ecosystem respiration (Reco) can identify old soil carbon (C) loss and serve as an early indicator of permafrost destabilization in a warming climate. Warming also stimulates plant productivity causing plant respiration to dominate Reco Δ 14 C signatures and potentially obscuring old soil C loss. Here, we investigate how a wide spatio-temporal gradient of permafrost thaw and plant productivity affects Reco Δ 14 C patterns and isotopic partitioning. Spatial gradients came from a warming experiment with doubling thaw depth and variable biomass, and a vegetation removal manipulation to eliminate plant contributions. We sampled in August and September to capture transitions from high to low plant productivity, decreased surface soil temperature, and relatively small seasonal thaw extensions. We found that surface processes dominate spatial variation in old soil C loss and a process-based partitioning approach was crucial for constraining old soil C loss. Resampling the same plots in different times of the year revealed that old soil C losses tripled with cooling surface temperature, and the largest old soil C losses were detected when the organic-to-mineral soil horizons thawed (~50–60 cm). We suggest that the measured increase in old soil respiration over the season and when the organic-to-mineral horizon thawed, may be explained by mobilization of nitrogen that stimulates microbial decomposition at depth. Our results suggest that soil C in the organic to mineral horizon may be an important source of soil C loss as the entire Arctic region warms and could lead to nonlinearities in projected permafrost climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

The coherent backscattering opposition effect

We have measured the opposition effect, the nonlinear surge in reflectance seen in particulate materials as phase angle approaches zero degrees, in a suite of materials of varying particle size and reflectance. These samples were illuminated by linearly and circularly polarized monochromatic radiation at two wavelengths, 0.442 and 0.633 microns. By measuring the linear and circular polarization ratios for each sample, we have found that in highly reflective materials the behavior of the reflected radiation is consistent with the coherent backscattering process which has recently been proposed to explain the opposition surge that is seen in such media. The size and width of the coherent backscattering opposition peak vary as a function of reflectance of the sample. The opposition effect has been observed in particulate materials studied in the laboratory and it is also observed in the radiation reflected from solar system bodies which present a regolith to the earth based observer. The traditional explanation of the opposition effect, the shadow-hiding hypothesis, is that it was caused by the elimination of mutual shadows cast between the regolith grains as the phase angle of the observation became smaller. This shadow-hiding hypothesis, however, is unable to explain the opposition effect seen in highly reflective materials such as magnesium oxide and barium sulfate powders. This is because highly reflective media will multiply scatter the incident radiation between the regolith grains. This causes the shadows to be eliminated. We have measured the angular scattering properties of a suite of materials of different reflectivity. We have observed polarization ratios in reflective particulates that are consistent with coherent backscattering as the principal process which causes the opposition surge.

Nelson, R. M.↗

A Compact, Tunable Near-UV Source for Quantitative Microgravity Combustion Diagnostics

There is a need for improved optical diagnostic methods for use in microgravity combustion research. Spectroscopic methods with fast time response that can provide absolute concentrations and concentration profiles of important chemical species in flames are needed to facilitate the understanding of combustion kinetics in microgravity. Although a variety of sophisticated laser-based diagnostics (such as planar laser induced fluorescence, degenerate four wave mixing and coherent Raman methods) have been applied to the study of combustion in laboratory flames, the instrumentation associated with these methods is not well suited to microgravity drop tower or space station platforms. Important attributes of diagnostic systems for such applications include compact size, low power consumption, ruggedness, and reliability. We describe a diode laser-based near-UV source designed with the constraints of microgravity research in mind. Coherent light near 420 nm is generated by frequency doubling in a nonlinear crystal. This light source is single mode with a very narrow bandwidth suitable for gas phase diagnostics, can be tuned over several 1/cm and can be wavelength modulated at up to MHz frequencies. We demonstrate the usefulness of this source for combustion diagnostics by measuring CH radical concentration profiles in an atmospheric pressure laboratory flame. The radical concentrations are measured using wavelength modulation spectroscopy (WMS) to obtain the line-of-sight integrated absorption for different paths through the flame. Laser induced fluorescence (LIF) measurements are also demonstrated with this instrument, showing the feasibility of simultaneous WMS absorption and LIF measurements with the same light source. LIF detection perpendicular to the laser beam can be used to map relative species densities along the line-of-sight while the integrated absorption available through WMS provides a mathematical constraint on the extraction of quantitative information from the LIF data. Combining absorption with LIF - especially if the measurements are made simultaneously with the same excitation beam - may allow elimination of geometrical factors and effects of intensity fluctuations (common difficulties with the analysis of LIF data) from the analysis.

Peterson, K. A.↗

Development of a near-wall Reynolds-stress closure based on the SSG model for the pressure strain

In this research, a near-wall second-order closure based on the Speziable et al.(1991) or SSG model for the pressure-strain term is proposed. Unlike the LRR model, the SSG model is quasi-nonlinear and yields better results when applied to calculate rotating homogeneous turbulent flows. An asymptotic analysis near the wall is applied to both the exact and modeled, equations so that appropriate near-wall corrections to the SSG model and the modeled dissipation-rate equation can be derived to satisfy the physical wall boundary conditions as well as the asymptotic near-wall behavior of the exact equations. Two additional model constants are introduced and they are determined by calibrating against one set of near-wall channel flow data. Once determined, their values are found to remain constant irrespective of the type of flow examined. The resultant model is used to calculate simple turbulent flows, near separating turbulent flows, complex turbulent flows and compressible turbulent flows with a freestream Mach number as high as 10. In all the flow cases investigated, the calculated results are in good agreement with data. This new near-wall model is less ad hoc, physically and mathematically more sound and eliminates the empiricism introduced by Zhang. Therefore, it is quite general, as demonstrated by the good agreement achieved with measurements covering a wide range of Reynolds numbers and Mach numbers.

So, R. M. C.↗

Effect of Scattering on Instabilities in Line-driven Stellar Winds

Line driven O-B stellar winds are unstable to perturbations of short spatial wavelength and the growth rates for such instabilities are very rapid. The nature of the nonlinear development of this instability is unknown, but might possibly be one of blobs of gas driven through ambient gas or a quasiregular train of outward moving shocks. In either case the resulting dissipation of mechanical energy might explain the observed anomalous heating n O-B stars as evidenced by their X-ray emission and high ionization state. It also might explain the observed fine structure of the absorption lines and their time variability. The driving due to the absorption of the stellar continuum flux was considered and the effects of the diffuse, scattered radiation field were neglected. It was shown that under certain special conditions the effect of scattering could reduce the instability growth rate to zero. A stability analysis that includes scattering, but that uses the more physically realistic assumption of complete redistribution instead of coherent scattering, and that includes the effects of transverse velocity gradients, which become important as the flow moves away from the stellar surface is presented. It is found that the instability is eliminated right at the base of the wind, but that as the flow moves outward the instability rate rapidly becomes equal to a substantial fraction of the calculated value, the fraction asymptotically reaching 80% at large radii. Since this still implies many e folds in a characteristic outflow time, the primary conclusion that these winds are highly unstable is unchanged.

Owocki, S. P.↗

Acceleration of electrons and ions by strong lower-hybrid turbulence in solar flares

One of the outstanding problems in solar flare theory is how to explain the 10-20 keV and greater hard x-ray emissions by a thick target bremsstrahlung model. The model requires the acceleration mechanism to accelerate approximately 10(exp 35) electrons sec(exp -l) with comparable energies, without producing a large return current which persists for long time scales after the beam ceases to exist due to Lenz's law, thereby, producing a self-magnetic field of order a few mega-Gauss. In this paper, we investigate particle acceleration resulting from the relaxation of unstable ion ring distributions, producing strong wave activity at the lower hybrid frequency. It is shown that strong lower hybrid wave turbulence collapses in configuration space producing density cavities containing intense electrostatic lower hybrid wave activity. The collapse of these intense nonlinear wave packets saturate by particle acceleration producing energetic electron and ion tails. There are several mechanisms whereby unstable ion distributions could be formed in the solar atmosphere, including reflection at perpendicular shocks, tearing modes, and loss cone depletion. Numerical simulations of ion ring relaxation processes, obtained using a 2 1/2-D fully electromagnetic, relativistic particle in cell code are discussed. We apply the results to the problem of explaining energetic particle production in solar flares. The results show the simultaneous acceleration of both electrons and ions to very high energies: electrons are accelerated to energies in the range 10-500 keV, while ions are accelerated to energies of the order of MeVs, giving rise to x-ray emission and gamma-ray emission respectively. Our simulations also show wave generation at the electron cyclotron frequency. We suggest that these waves are the solar millisecond radio spikes. The strong turbulence collapse process leads to a highly filamented plasma producing many localized regions for particle acceleration and resulting in approximately 10(exp 17) electron 'beamlets' of width approximately equal to 10 lambda sub De which eliminates the production of large magnetic fields. In this paper, we demonstrate that the model produces an energetic electron spectrum with the right flux to account for the hard x-ray observations.

Spicer, D. S.↗

Cascade Optimization Strategy for Aircraft and Air-Breathing Propulsion System Concepts

Design optimization for subsonic and supersonic aircraft and for air-breathing propulsion engine concepts has been accomplished by soft-coupling the Flight Optimization System (FLOPS) and the NASA Engine Performance Program analyzer (NEPP), to the NASA Lewis multidisciplinary optimization tool COMETBOARDS. Aircraft and engine design problems, with their associated constraints and design variables, were cast as nonlinear optimization problems with aircraft weight and engine thrust as the respective merit functions. Because of the diversity of constraint types and the overall distortion of the design space, the most reliable single optimization algorithm available in COMETBOARDS could not produce a satisfactory feasible optimum solution. Some of 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 both aircraft and engines. The cascade method has two principal steps: In the first, the solution initiates from a user-specified design and optimizer, in the second, the optimum design obtained in the first step with some random perturbation is used to begin the next specified optimizer. The second step is repeated for a specified sequence of optimizers or until a successful solution of the problem is achieved. A successful solution should satisfy the specified convergence criteria and have several active constraints but no violated constraints. The cascade strategy available in the combined COMETBOARDS, FLOPS, and NEPP design tool converges to the same global optimum solution even when it starts from different design points. This reliable and robust design tool eliminates manual intervention in the design of aircraft and of air-breathing propulsion engines where it eases the cycle analysis procedures. The combined code is also much easier to use, which is an added benefit. This paper describes COMETBOARDS and its cascade strategy and illustrates the capability of the combined design tool through the optimization of a subsonic aircraft and a high-bypass-turbofan wave-rotor-topped engine.

Patnaik, Surya N.↗