Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NONLINEAR PROGRAMMING”

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,423 records · Page 79

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

An Overview of the Semi-Span Super-Sonic Transport (S4T) Wind-Tunnel Model Program

A summary of computational and experimental aeroelastic (AE) and aeroservoelastic (ASE) results for the Semi-Span Super-Sonic Transport (S4T) wind-tunnel model is presented. A broad range of analyses and multiple AE and ASE wind-tunnel tests of the S4T wind-tunnel model have been performed in support of the ASE element in the Supersonics Program, part of the NASA Fundamental Aeronautics Program. This paper is intended to be an overview of multiple papers that comprise a special S4T technical session. Along those lines, a brief description of the design and hardware of the S4T wind-tunnel model will be presented. Computational results presented include linear and nonlinear aeroelastic analyses, and rapid aeroelastic analyses using CFD-based reduced-order models (ROMs). A brief survey of some of the experimental results from two open-loop and two closed-loop wind-tunnel tests performed at the NASA Langley Transonic Dynamics Tunnel (TDT) will be presented as well.

Silva, Walter A.↗

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Composite materials research and education program: The NASA-Virginia Tech composites program

Major areas of study include: (1) edge effects in finite width laminated composites subjected to mechanical, thermal and hygroscopic loading with temperature dependent material properties and the influence of edge effects on the initiation of failure; (2) shear and compression testing of composite materials at room and elevated temperatures; (3) optical techniques for precise measurement of coefficients of thermal expansion of composites; (4) models for the nonlinear behavior of composites including material nonlinearity and damage accumulation and verification of the models under biaxial loading; (5) compressive failure of graphite/epoxy plates with circular holes and the buckling of composite cylinders under combined compression and torsion; (6) nonlinear mechanical properties of borsic/aluminum, graphite/polyimide and boron/aluminum; (7) the strength characteristics of spliced sandwich panels; and (8) curved graphite/epoxy panels subjected to internal pressure.

Herakovich, C. T.↗

Domain-decomposition nonlinear manifold reduced order model

This software combines nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD) techniques. NM-ROMs, which utilize a shallow, sparse autoencoder trained with full order model (FOM) snapshot data, approximate the FOM state on a nonlinear manifold. These models offer advantages over linear-subspace ROMs (LS-ROMs) particularly in scenarios with slowly decaying Kolmogorov n-width. However, the training of NM-ROMs involves a number of parameters that scale with the size of the FOM, and storing high-dimensional FOM snapshots can significantly increase the cost of ROM training for extreme-scale problems. To mitigate these costs, the software employs DD to partition the FOM into smaller subdomains, computes NM-ROMs for each, and then integrates these to form a global NM-ROM. This strategy offers multiple benefits: it enables parallel training of subdomain NM-ROMs, reduces the number of parameters needed, decreases the dimensional requirements of subdomain FOM training data, and allows for customization to the unique characteristics of each FOM subdomain. The use of a shallow, sparse autoencoder architecture in each subdomain NM-ROM facilitates the application of hyper-reduction (HR), simplifying the nonlinear complexities and enhancing computational speed. This software marks the inaugural application of NM-ROM combined with HR to a DD problem. It features an algebraic DD reformulation of the FOM, training of NM-ROMs with HR for each subdomain, and employs a sequential quadratic programming (SQP) solver for the evaluation of the coupled global NMROM. The effectiveness of the DD NM-ROM with HR is numerically demonstrated on the 2D steady-state Burgers' equation, showing an order of magnitude improvement in accuracy over the DD LS-ROM with HR.

Diaz, AlejandroN↗

Unsteady aerodynamic flow field analysis of the space shuttle configuration. Part 2: Launch vehicle aeroelastic analysis

An exploratory analysis has been made of the aeroelastic stability of the Space Shuttle Launch Configuration, with the objective of defining critical flow phenomena with adverse aeroelastic effects and developing simple analytic means of describing the time-dependent flow-interference effects so that they can be incorporated into a computer program to predict the aeroelastic stability of all free-free modes of the shuttle launch configuration. Three critical flow phenomana have been identified: (1) discontinuous jump of orbiter wing shock, (2) inlet flow between orbiter and booster, and (3) H.O. tank base flow. All involve highly nonlinear and often discontinuous aerodynamics which cause limit cycle oscillations of certain critical modes. Given the appropriate static data, the dynamic effects of the wing shock jump and the HO tank bulbous base effect can be analyzed using the developed quasi-steady techniques. However, further analytic and experimental efforts are required before the dynamic effects of the inlet flow phenomenon can be predicted for the shuttle launch configuration.

Reding, J. P.↗

Evaluation of inelastic constitutive models for nonlinear structural analysis

The influence of inelastic material models on computed stress-strain states, and therefore predicted lives, was studied for thermomechanically loaded structures. Nonlinear structural analyses were performed on a fatigue specimen which had been subjected to thermal cycling in fluidized beds and on a mechanically load cycled benchmark notch specimen. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, combined plus transient creep) were exercised using the MARC program. Of the plasticity models, kinematic hardening gave results most consistent with experimental observations. Life predictions using the computed strain histories at the critical location with a strainrange partitioning approach considerably overpredicted the crack initiation life of the thermal fatigue specimen.

Kaufman, A.↗

Performance of a quantitative study of instability-related delamination growth

An experimental program was conducted to quantify instability-induced imbedded delamination growth. Static tests on double cantilever beam (DCB) specimens yielded the critical mode 1 strain energy release rate (G sub 1 C) for T300/5208 graphite/epoxy. Static tests on mixed mode cracked lap shear (CLS) specimens, and a nonlinear finite element analysis (NFEA) of the CLS specimen to separate mode 1 and mode 2 effects, yielded the critical mode 2 strain energy release rate (G sub 2 C) for T300/5208. Constant amplitude fatigue tests on DCB and CLS specimens, along with the NFEA results on CLS specimens, quantified mode 1 and mode 2 contributions to delamination growth rate. Fatigue tests were conducted at a frequency (omega) of 10 Hertz, maintaining the minimum to maximum cyclic load ratio (R) at 0.05. Static compression and constant amplitude compression fatigue tests were also conducted on specimens with imbedded through the width (ITTW) delaminations. Kapton imbeddments were located below 3, 4 or 6 plies in a 64-ply laminate, during layup, to simulate ITTW delaminations.

Ramkumar, R. L.↗

GaussFit - A system for least squares and robust estimation

GaussFit is a new computer program for solving least-squares and robust estimation problems. It has a number of unique features, including a complete programming language designed especially to formulate estimation problems, a built-in compiler and interpreter to support the programming language, and a built-in algebraic manipulator for calculating the required partial derivatives analytically. These features make GaussFit very easy to use, so that even complex problems can be set up and solved with minimal effort. GaussFit can correctly handle many cases of practical interest: nonlinear models, exact constraints, correlated observations, and models where the equations of condition contain more than one observed quantity. An experimental robust estimation capability is built into GaussFit so that data sets contaminated by outliers can be handled simply and efficiently.

Jefferys, W. H.↗

Component specific modeling

The overall objective of this program is to develop and verify a series of interdisciplinary modeling and analysis techniques that were specialized to address three specific hot section components. These techniques incorporate data as well as theoretical methods from many diverse areas including cycle and performance analysis, heat transfer analysis, linear and nonlinear stress analysis, and mission analysis. Building on the proven techniques already available in these fields, the new methods developed are integrated to provide an accurate, efficient, and unified approach to analyzing combustor burner liners, hollow air-cooled turbine blades, and air-cooled turbine vanes. For these components, the methods developed predict temperature, deformation, stress, and strain histories throughout a complete flight mission.

Maffeo, R. J.↗

Comparison of Rolling Moment Characteristics During Roll Oscillations for a Low and a High Aspect Ratio Configuration

Improvements in testing and modeling of nonlinear and unsteady aerodynamic effects for flight dynamics predictions of vehicle performance is critical to enable the design and implementation of new, innovative vehicle concepts. Any configuration which exhibits significant flow separation, nonlinear aerodynamics, control interactions or attempts maneuvering through one or more conditions such as these is, at present, a challenge to test, model or predict flight dynamic responses prior to flight. Even in flight test experiments, adequate models are not available to study and characterize the complex nonlinear and time-dependent flow effects occurring during portions of the maneuvering envelope. Traditionally, airplane designs have been conducted to avoid these areas of the flight envelope. Better understanding and characterization of these flight regimes may not only reduce risk and cost of flight test development programs, but also may pave the way for exploitation of those characteristics that increase airplane capabilities. One of the hurdles is that the nonlinear/unsteady effects appear to be configuration dependent. This paper compares some of the dynamic aerodynamic stability characteristics of two very different configurations - representative of a fighter and a transport airplane - during dynamic body-axis roll wind tunnel tests. The fighter model shows significant effects of oscillation frequency which are not as apparent for the transport configuration.

Brandon, Jay M.↗

Thermal barrier coating life prediction model development, phase 2

The objective of this program was to generate a life prediction model for electron-beam-physical vapor deposited (EB-PVD) zirconia thermal barrier coating (TBC) on gas turbine engine components. Specific activities involved in development of the EB-PVD life prediction model included measurement of EB-PVD ceramic physical and mechanical properties and adherence strength, measurement of the thermally grown oxide (TGO) growth kinetics, generation of quantitative cyclic thermal spallation life data, and development of a spallation life prediction model. Life data useful for model development was obtained by exposing instrumented, EB-PVD ceramic coated cylindrical specimens in a jet fueled burner rig. Monotonic compression and tensile mechanical tests and physical property tests were conducted to obtain the EB-PVD ceramic behavior required for burner rig specimen analysis. As part of that effort, a nonlinear constitutive model was developed for the EB-PVD ceramic. Spallation failure of the EB-PVD TBC system consistently occurred at the TGO-metal interface. Calculated out-of-plane stresses were a small fraction of that required to statically fail the TGO. Thus, EB-PVD spallation was attributed to the interfacial cracking caused by in-plane TGO strains. Since TGO mechanical properties were not measured in this program, calculation of the burner rig specimen TGO in-plane strains was performed by using alumina properties. A life model based on maximum in-plane TGO tensile mechanical strain and TGO thickness correlated the burner rig specimen EB-PVD ceramic spallation lives within a factor of about plus or minus 2X.

Meier, Susan Manning↗

User's guide for ENSAERO: A multidisciplinary program for fluid/structural/control interaction studies of aircraft (release 1)

Strong interactions can occur between the flow about an aerospace vehicle and its structural components resulting in several important aeroelastic phenomena. These aeroelastic phenomena can significantly influence the performance of the vehicle. At present, closed-form solutions are available for aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for aeroelasticity involving complex nonlinear flows with shock waves, vortices, flow separations, and aerodynamic heating, computational methods are still under development. These complex aeroelastic interactions can be dangerous and limit the performance of aircraft. Examples of these detrimental effects are aircraft with highly swept wings experiencing vortex-induced aeroelastic oscillations, transonic regime at which the flutter speed is low, aerothermoelastic loads that play a critical role in the design of high-speed vehicles, and flow separations that often lead to buffeting with undesirable structural oscillations. The simulation of these complex aeroelastic phenomena requires an integrated analysis of fluids and structures. This report presents a summary of the development, applications, and procedures to use the multidisciplinary computer code ENSAERO. This code is based on the Euler/Navier-Stokes flow equations and modal/finite-element structural equations.

Guruswamy, Guru P.↗

The Physics of Gain Mechanisms in a Self-Amplified Spontaneous Emission Free-Electron Laser

This is the final report for the DoE BES contract DE-FG02-07ER46272. This grant began under a title that gives its history in a few words: The Physics of Gain Mechanisms in a Self-Amplified Spontaneous Emission Free-Electron Laser. At the time of adopting this name, the UCLA Particle Beam Physics Laboratory (PBPL) was at the forefront of the first demonstrations of SASE FEL, giving the experimental basis for the 4th generation X-ray light source now exemplified by the LCLS, a remarkable generation of machines that have revolutionize nanoscale and ultrafast science. The grant title indicates that physical aspects of the FEL are of interest, and implies that progress requires comprehensive investigations. These aspects include: the physics and technology of high brightness electron beam generation, phase space manipulations, advanced beam diagnostics, state-of-the-art undulators, studies of FEL mode characteristics, as well as the introduction of new theoretical models and start-to-end simulations, all yielding a deep understanding of the SASE FEL. The latest period (2010-13) of BES grant support has been the most productive yet for the UCLA PBPL, and the extensive list of accomplishments is illustrative of our central contributions to advanced FEL research. We give a concise overview here and review our research results in further detail below. In FEL and ultrafast beam physics, we have experimentally investigated: production of helical beams and FEL modes with orbital angular momentum; chirp-taper single-spike FEL gain; new FROG techniques optimized for FEL, introduced coherent imaging of beam profiles using transition radiation; and demonstrated a next generation cryoundulator. We have performed deep theoretical studies of these phenomena and other advanced FEL physics concepts, such as the space-charge amplifier. In high brightness beams, we have: introduced a new class of hybrid photoinjector for <100 fs beam production; developed new approaches to self-shaped ellipsoidal beams; used nonlinear plasma oscillations to create beam pulse trains; tested new plasmonic, high quantum efficiency cathodes, and developed new methods of beam measurement with um spatial/sub-fs temporal resolution. In ultrafast imaging, we have: developed advanced electro-optic timing systems, performed the 1st pump-probe study using single shot UED, demonstrated continuously time-resolved UED with RF deflection; made initial use of inverse Compton scattering (ICS) X-ray sources in application through single-shot inline phase contrast imaging, established single-shot X-ray diffraction, and observed ICS high harmonic generation and redshifting due to nonlinear effects. In 5th generation light source research, we have developed the techniques of narrowband, very high power coherent Cerenkov radiation sources based on dielectric wakefield devices. This work is represented by numerous publications in leading journals: e.g. Nature Physics, Physical Review Letters, Applied Physics Letters. The PBPL program’s influence on the field also extends to organization of the electron beam applications community; we have organized the first international UED workshops and the initial two meetings dedicated to 5th generation light source.

43 PARTICLE ACCELERATORS↗

Structural model optimization using statistical evaluation

The results of research in applying statistical methods to the problem of structural dynamic system identification are presented. The study is in three parts: a review of previous approaches by other researchers, a development of various linear estimators which might find application, and the design and development of a computer program which uses a Bayesian estimator. The method is tried on two models and is successful where the predicted stiffness matrix is a proper model, e.g., a bending beam is represented by a bending model. Difficulties are encountered when the model concept varies. There is also evidence that nonlinearity must be handled properly to speed the convergence.

Collins, J. D.↗

Design considerations for flight test of a fault inferring nonlinear detection system algorithm for avionics sensors

The modifications to the design of a fault inferring nonlinear detection system (FINDS) algorithm to accommodate flight computer constraints and the resulting impact on the algorithm performance are summarized. An overview of the flight data-driven FINDS algorithm is presented. This is followed by a brief analysis of the effects of modifications to the algorithm on program size and execution speed. Significant improvements in estimation performance for the aircraft states and normal operating sensor biases, which have resulted from improved noise design parameters and a new steady-state wind model, are documented. The aircraft state and sensor bias estimation performances of the algorithm's extended Kalman filter are presented as a function of update frequency of the piecewise constant filter gains. The results of a new detection system strategy and failure detection performance, as a function of gain update frequency, are also presented.

Caglayan, A. K.↗

The Krigifier: A Procedure for Generating Pseudorandom Nonlinear Objective Functions for Computational Experimentation

Comprehensive computational experiments to assess the performance of algorithms for numerical optimization require (among other things) a practical procedure for generating pseudorandom nonlinear objective functions. We propose a procedure that is based on the convenient fiction that objective functions are realizations of stochastic processes. This report details the calculations necessary to implement our procedure for the case of certain stationary Gaussian processes and presents a specific implementation in the statistical programming language S-PLUS.

Trosset, Michael W.↗