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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Additive Manufacturing of Biomechanically Tailored Meshes for Compliant Wearable and Implantable Devices

Abstract Additive manufacturing (AM) of medical devices such as orthopedic implants and hearing aids is highly attractive because of the potential of AM to match the complex form and mechanics of individual human bodies. Externally worn and implantable tissue‐support devices, such as ankle or knee braces, and hernia repair mesh, offer a new opportunity for AM to mimic tissue‐like mechanics and improve both patient outcomes and comfort. Here, it is demonstrated how explicit programming of the toolpath in an extrusion AM process can enable new, flexible mesh materials having digitally tailored mechanical properties and geometry. Meshes are fabricated by extrusion of thermoplastics, optionally with continuous fiber reinforcement, using a continuous toolpath that tailors the elasticity of unit cells of the mesh via incorporation of slack and modulation of filament–filament bonding. It is shown how the tensile mesh mechanics can be engineered to match the nonlinear response of muscle. An ankle brace with directionally specific inversion stiffness arising from embedded mesh is validated, and further concepts for 3D mesh devices are prototyped.

Pattinson, Sebastian W.↗
NASA NTRS1995

Supersonic civil airplane study and design: Performance and sonic boom

Since aircraft configuration plays an important role in aerodynamic performance and sonic boom shape, the configuration of the next generation supersonic civil transport has to be tailored to meet high aerodynamic performance and low sonic boom requirements. Computational fluid dynamics (CFD) can be used to design airplanes to meet these dual objectives. The work and results in this report are used to support NASA's High Speed Research Program (HSRP). CFD tools and techniques have been developed for general usages of sonic boom propagation study and aerodynamic design. Parallel to the research effort on sonic boom extrapolation, CFD flow solvers have been coupled with a numeric optimization tool to form a design package for aircraft configuration. This CFD optimization package has been applied to configuration design on a low-boom concept and an oblique all-wing concept. A nonlinear unconstrained optimizer for Parallel Virtual Machine has been developed for aerodynamic design and study.

Cheung, Samson↗
NASA NTRS2002

Unsteady Pressures in a Transonic Fan Cascade Due to a Single Oscillating Airfoil

An extensive set of unsteady pressure data was acquired along the midspan of a modern transonic fan blade for simulated flutter conditions. The data set was acquired in a nine-blade linear cascade with an oscillating middle blade to provide a database for the influence coefficient method to calculate instantaneous blade loadings. The cascade was set for an incidence of 10 dg. The data were acquired on three stationary blades on each side of the middle blade that was oscillated at an amplitude of 0.6 dg. The matrix of test conditions covered inlet Mach numbers of 0.5, 0.8, and 1.1 and the oscillation frequencies of 200, 300, 400, and 500 Hz. A simple quasiunsteady two-dimensional computer simulation was developed to aid in the running of the experimental program. For high Mach number subsonic inlet flows the blade pressures exhibit very strong, low-frequency, self-induced oscillations even without forced blade oscillations, while for low subsonic and supersonic inlet Mach numbers the blade pressure unsteadiness is quite low. The amplitude of forced pressure fluctuations on neighboring stationary blades strongly depends on the inlet Mach number and forcing frequency. The flowfield behavior is believed to be governed by strong nonlinear effects due to a combination of viscosity, compressibility, and unsteadiness. Therefore, the validity of the quasi-unsteady simplified computer simulation is limited to conditions when the flowfield is behaving in a linear, steady manner. Finally, an extensive set of unsteady pressure data was acquired to help development and verification of computer codes for blade flutter effects.

Lepicovsky, J.↗
DOE OSTI2015

Method for conducting nonlinear electrochemical impedance spectroscopy

A method for conducting nonlinear electrochemical impedance spectroscopy. The method includes quantifying the nonlinear response of an electrochemical system by measuring higher-order current or voltage harmonics generated by moderate-amplitude sinusoidal current or voltage perturbations. The method involves acquisition of the response signal followed by time apodization and fast Fourier transformation of the data into the frequency domain, where the magnitude and phase of each harmonic signal can be readily quantified. The method can be implemented on a computer as a software program.

25 ENERGY STORAGE↗
DOE OSTI2019

Designing Metabolic Division of Labor in Microbial Communities

ABSTRACT Microbes face a trade-off between being metabolically independent and relying on neighboring organisms for the supply of some essential metabolites. This balance of conflicting strategies affects microbial community structure and dynamics, with important implications for microbiome research and synthetic ecology. A “gedanken” (thought) experiment to investigate this trade-off would involve monitoring the rise of mutual dependence as the number of metabolic reactions allowed in an organism is increasingly constrained. The expectation is that below a certain number of reactions, no individual organism would be able to grow in isolation and cross-feeding partnerships and division of labor would emerge. We implemented this idealized experiment using in silico genome-scale models. In particular, we used mixed-integer linear programming to identify trade-off solutions in communities of Escherichia coli strains. The strategies that we found revealed a large space of opportunities in nuanced and nonintuitive metabolic division of labor, including, for example, splitting the tricarboxylic acid (TCA) cycle into two separate halves. The systematic computation of possible solutions in division of labor for 1-, 2-, and 3-strain consortia resulted in a rich and complex landscape. This landscape displayed a nonlinear boundary, indicating that the loss of an intracellular reaction was not necessarily compensated for by a single imported metabolite. Different regions in this landscape were associated with specific solutions and patterns of exchanged metabolites. Our approach also predicts the existence of regions in this landscape where independent bacteria are viable but are outcompeted by cross-feeding pairs, providing a possible incentive for the rise of division of labor. IMPORTANCE Understanding how microbes assemble into communities is a fundamental open issue in biology, relevant to human health, metabolic engineering, and environmental sustainability. A possible mechanism for interactions of microbes is through cross-feeding, i.e., the exchange of small molecules. These metabolic exchanges may allow different microbes to specialize in distinct tasks and evolve division of labor. To systematically explore the space of possible strategies for division of labor, we applied advanced optimization algorithms to computational models of cellular metabolism. Specifically, we searched for communities able to survive under constraints (such as a limited number of reactions) that would not be sustainable by individual species. We found that predicted consortia partition metabolic pathways in ways that would be difficult to identify manually, possibly providing a competitive advantage over individual organisms. In addition to helping understand diversity in natural microbial communities, our approach could assist in the design of synthetic consortia.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

NAS (Numerical Aerodynamic Simulation Program) technical summaries, March 1989 - February 1990

Given here are selected scientific results from the Numerical Aerodynamic Simulation (NAS) Program's third year of operation. During this year, the scientific community was given access to a Cray-2 and a Cray Y-MP supercomputer. Topics covered include flow field analysis of fighter wing configurations, large-scale ocean modeling, the Space Shuttle flow field, advanced computational fluid dynamics (CFD) codes for rotary-wing airloads and performance prediction, turbulence modeling of separated flows, airloads and acoustics of rotorcraft, vortex-induced nonlinearities on submarines, and standing oblique detonation waves.

Source record↗
NASA NTRS1993

Growth of new materials for solid state laser applications. The growth of ZnGeP2 by the vertical Bridgman method

This is the final technical report on the cooperative NASA program 'Growth of New Materials for Solid State Laser Applications,' covering the period from 9-1-86 through 3-31-91. The first two and one half years of the program, from 9-1-86 through 3-31-89, was devoted to the development of new eye-safe laser sources. Single crystal fibers of rare earth doped and co-doped YAG were grown by the laser-heated pedestal growth method, characterized for their structural properties and supplied to NASA/Langley for spectroscopic evaluation. From 4-1-89 though 3-31-91, research focused on the growth of zinc germanium phosphide for nonlinear optical applications in the IR. The vertical Bridgman growth process was evaluated as a means to achieve better control over the short wavelength optical absorption in this material that adversely affects 2 micron-pumped optical parametric oscillators.

Feigelson, R. S.↗
NASA NTRS1999

Crack Growth Simulation and Residual Strength Prediction in Airplane Fuselages

This is the final report for the NASA funded project entitled "Crack Growth Prediction Methodology for Multi-Site Damage." The primary objective of the project was to create a capability to simulate curvilinear fatigue crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage. The second objective was to validate the capability by way of comparisons to experimental results. Both objectives have been achieved and the results are detailed herein. In the first part of the report, the crack tip opening angle (CTOA) fracture criterion, obtained and correlated from coupon tests to predict fracture behavior and residual strength of built-up aircraft fuselages, is discussed. Geometrically nonlinear, elastic-plastic, thin shell finite element analyses are used to simulate stable crack growth and to predict residual strength. Both measured and predicted results of laboratory flat panel tests and full-scale fuselage panel tests show substantial reduction of residual strength due to the occurrence of multi-site damage (MSD). Detailed comparisons of n stable crack growth history, and residual strength between the predicted and experimental results are used to assess the validity of the analysis methodology. In the second part of the report, issues related to crack trajectory prediction in thin shells; an evolving methodology uses the crack turning phenomenon to improve the structural integrity of aircraft structures are discussed, A directional criterion is developed based on the maximum tangential stress theory, but taking into account the effect of T-stress and fracture toughness orthotropy. Possible extensions of the current crack growth directional criterion to handle geometrically and materially nonlinear problems are discussed. The path independent contour integral method for T-stress evaluation is derived and its accuracy is assessed using a p- and hp-version adaptive finite element method. Curvilinear crack growth is simulated in coupon tests and in full-scale fuselage panel tests. Both T-stress and fracture toughness orthotropy are found to be essential to predict the observed crack paths. The analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically while insuring continuous airworthiness. Consequently, it will improve the technology to support the safe operation of the current aircraft fleet as well as the design of more damage-tolerant aircraft for the next generation fleet.

Chen, Chuin-Shan↗
NASA NTRS1971

Nonlinear, nonlaminar-3D computation of electron motion through the output cavity of a klystron

The equations of motion used in the computation are discussed along with the space charge fields and the integration process. The following assumptions were used as a basis for the computation: (1) The beam is divided into N axisymmetric discs of equal charge and each disc into R rings of equal charge. (2) The velocity of each disc, its phase with respect to the gap voltage, and its radius at a specified position in the drift tunnel prior to the interaction gap is known from available large signal one dimensional programs. (3) The fringing rf fields are computed from exact analytical expressions derived from the wave equation assuming a known field shape between the tunnel tips at a radius a. (4) The beam is focused by an axisymmetric magnetic field. Both components of B, that is B sub z and B sub r, are taken into account. (5) Since this integration does not start at the cathode but rather further down the stream prior to entering the output cavity it is assumed that each electron moved along a laminar path from the cathode to the start of integration.

Albers, L. U.↗
DOE OSTI2026

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗
DOE OSTI2014

Understanding the many-body expansion for large systems. I. Precision considerations

Electronic structure methods based on low-order “n-body” expansions are an increasingly popular means to defeat the highly nonlinear scaling of ab initio quantum chemistry calculations, taking advantage of the inherently distributable nature of the numerous subsystem calculations. Here, we examine how the finite precision of these subsystem calculations manifests in applications to large systems, in this case, a sequence of water clusters ranging in size up to (H{sub 2}O){sub 47}. Using two different computer implementations of the n-body expansion, one fully integrated into a quantum chemistry program and the other written as a separate driver routine for the same program, we examine the reproducibility of total binding energies as a function of cluster size. The combinatorial nature of the n-body expansion amplifies subtle differences between the two implementations, especially for n ⩾ 4, leading to total energies that differ by as much as several kcal/mol between two implementations of what is ostensibly the same method. This behavior can be understood based on a propagation-of-errors analysis applied to a closed-form expression for the n-body expansion, which is derived here for the first time. Discrepancies between the two implementations arise primarily from the Coulomb self-energy correction that is required when electrostatic embedding charges are implemented by means of an external driver program. For reliable results in large systems, our analysis suggests that script- or driver-based implementations should read binary output files from an electronic structure program, in full double precision, or better yet be fully integrated in a way that avoids the need to compute the aforementioned self-energy. Moreover, four-body and higher-order expansions may be too sensitive to numerical thresholds to be of practical use in large systems.

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

Analysis of SNL/MSU/DOE Fatigue Database Trends for Wind Turbine Blade Materials 2010-2015.

Wind turbine blades are designed to several major structural conditions, including tip deflection, strength and b uckling during severe loading, as well as very high numbers of fatigue cycles and various service environments. The MSU Database Program has, since 1989, addressed the broad range of properties needed for current and potential blade materials through stati c and fatigue testing and test development in cooperation with Sandia National Laboratories and wind industry and supplier partners. This report is the latest in a series, giving test results and analysis for the period 2010 - 2015. Program data are compiled in a public database [1] and other reports and publications given in the cited references. The report begins with an executive summary and introductory material including background discussion of previous related studies. Section 3 describes experimental methods including processing, test methods, instrumentation and test development. Section 4 provides static tension, compression and shear stress - strain properties in three directions using coupons sectioned from a thick infused unidirectional glass/epoxy laminate. The nonlinear, shear dominated static properties were characterized with loading - u nloading - reloading (LUR) tests in tension and compression to increasing load levels, for +-45O laminates. Section 5 explores the origins of tensile fatigue sensitivity in glass fiber dominated laminates with variations in fabric architecture including speci ally prepared fabrics and aligned strand laminates. Several types of resins are considered, with variations in resin toughness and bonding to fibers, as well as cure cycle variations for an epoxy. Conclusions are drawn as to the limits of tensile fatigue r esistance and the effects of resin type and fabric architecture, including the behavior of a commercial aligned glass strand product. Interactions between cyclic fatigue response and creep are addressed for off - axis (+-45O) glass/epoxy laminates in Sectio n 6. The nonlinear fatigue and creep stress - strain and cumulative strain response are characterized in tension and compression as a function of stress level, cycles and cumulative time, using square and sinewave loading over a broad range of frequency. The results are analyzed in terms of the cycles and cumulative time under load. A cumulative strain failure criterion is established, and used to construct shear and tension constant life diagrams (CLD's) with data for nine R - values. The effects of a more duc tile urethne resin are also explored. A previous study of thick adhesives testing is extended to mixed mode fracture mechanics testing in Section 7. Mechanisms of static and fatigue crack extension near the laminate adherend interface are reported in deta il. Data are presented for mixed mode adhesive fracture, compared to mixed mode fracture in ply delamination. Fatigue crack growth exponents are also developed for a mixed mode cracked lap shear coupon. The data for fatigue trends and relative failure stra ins and exponents are compared for various blade component materials in Section 8. The effects of temperature and seawater saturation are considered for selected materials of interest for wind and hydrokinetic turbine blades in Section 9. Section 10 gives detailed conclusions for each section. A cknowledgements The research presented in this report was carried out under Sandia National Laboratories purchase orders 1325028 an d 1543945 between 2010 and 2015, with support from the DOE Wind and Water Technologies Office . In addition to the authors listed, significant contributions were made by Patrick Flaherty, Pancastya Agastra, Michael Schuster, and Michael Voth. Industry m aterials suppliers include Vectorply, Saertex, OCV, AGY, Bayer, Ashland, 3M and Nextel. Industry suppliers with significant contributions to the study were Hexion, PPG, Reichhold, Gurit and NEPTCO. Intentionally Left Blank

17 WIND ENERGY↗
NASA NTRS1973

Multiple regression technique for Pth degree polynominals with and without linear cross products

A multiple regression technique was developed by which the nonlinear behavior of specified independent variables can be related to a given dependent variable. The polynomial expression can be of Pth degree and can incorporate N independent variables. Two cases are treated such that mathematical models can be studied both with and without linear cross products. The resulting surface fits can be used to summarize trends for a given phenomenon and provide a mathematical relationship for subsequent analysis. To implement this technique, separate computer programs were developed for the case without linear cross products and for the case incorporating such cross products which evaluate the various constants in the model regression equation. In addition, the significance of the estimated regression equation is considered and the standard deviation, the F statistic, the maximum absolute percent error, and the average of the absolute values of the percent of error evaluated. The computer programs and their manner of utilization are described. Sample problems are included to illustrate the use and capability of the technique which show the output formats and typical plots comparing computer results to each set of input data.

Davis, J. W.↗
NASA NTRS1989

Vector algorithms for geometrically nonlinear 3D finite element analysis

Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.

Whitcomb, John D.↗
NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗
NASA NTRS1999

Delamination and Stitched Failure in Stitched Composite Joints

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and finite element study. The experimental program was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation of the debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The experimentally determined debond length vs. applied load was used as an input parameter in the finite element analysis of both configurations. The strain energy release rates at the debond from were calculated using plate finite elements. Nonlinear fastener elements were used to model the stitches and multipoint constraints were used to model the contact problem. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches were effective in reducing mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗
NASA NTRS2014

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel: Comparison of In- and Out-of-Autoclave Facesheet Configurations

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10-m diameter barrel section of the Heavy Lift Launch Vehicle, were manufactured and tested under the NASA Composites for Exploration and the NASA Constellation Ares V programs. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.0 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: in-autoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3 ft. by 5 ft. panels were cut from the 1/16th barrel sections and tested under compressive loading. Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of each 3 ft. by 5 ft. panel. To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yielded good qualitative and quantitative predictions. Additionally, it was correctly predicted that the panel would fail in buckling prior to failing in strength. Furthermore, several imperfection studies were performed to investigate the influence of geometric imperfections, fiber angle misalignments, and three-dimensional effects on the compressive response of the panel.

epoxy matrix composites↗