Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closed form solution”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Implications of Underground Nuclear Explosion Cavity Evolution for Radioxenon Isotopic Composition

Abstract Isotopic ratios of radioxenons sampled in the atmosphere or subsurface can be used to verify the occurrence of an underground nuclear explosion (UNE). Differences in the half-lives of radioactive xenon precursors and their decay-chain networks produce different time-dependent concentration profiles of xenon isotopes allowing isotopic ratios to be used for tracking UNE histories including estimating the time of detonation. In this study, we explore the potential effects of post-detonation cavity processes: precipitation of iodine precursors, gas seepage, and prompt venting on radioxenon isotopic evolution which influences UNE histories. Simplified analytical models and closed-form solutions yielding a potentially idealized radioactive decay/ingrowth chain in a closed and well-mixed system typically have limited application by not including the partitioning of the radionuclide inventory between a gas phase and rock melt created by the detonation and by ignoring gas transport from the cavity to host rock or ground surface. In reality, either subsurface transport or prompt release that is principally responsible for gas signatures violates the closed-system (or batch-mode) assumption. A closed-form solution representing time-dependent source-term activities is extended by considering the cavity partitioning process, slow seepage, and/or prompt release of gases from the cavity and applied to realistic systems.

Sun, Yunwei (ORCID:0000000194801439)↗
NASA NTRS2003

Spectral Characteristics of Wake Vortex Sound During Roll-Up

This report presents an analysis of the sound spectra generated by a trailing aircraft vortex during its rolling-up process. The study demonstrates that a rolling-up vortex could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB above threshold of human hearing) at a distance of 200 ft from the vortex core. The spectrum then drops o rapidly thereafter. A rigorous analytical approach has been adopted in this report to derive the spectrum of vortex sound. First, the sound pressure was solved from an alternative treatment of the Lighthill s acoustic analogy approach [1]. After the application of Green s function for free space, a tensor analysis was applied to permit the removal of the source term singularity of the wave equation in the far field. Consequently, the sound pressure is expressed in terms of the retarded time that indicates the time history and spacial distribution of the sound source. The Fourier transformation is then applied to the sound pressure to compute its spectrum. As a result, the Fourier transformation greatly simplifies the expression of the vortex sound pressure involving the retarded time, so that the numerical computation is applicable with ease for axisymmetric line vortices during the rolling-up process. The vortex model assumes that the vortex circulation is proportional to the time and the core radius is a constant. In addition, the velocity profile is assumed to be self-similar along the aircraft flight path, so that a benchmark vortex velocity profile can be devised to obtain a closed form solution, which is then used to validate the numerical calculations for other more realistic vortex profiles for which no closed form solutions are available. The study suggests that acoustic sensors operating at low frequency band could be profitably deployed for detecting the vortex sound during the rolling-up process.

Booth, Earl R., Jr.↗
NASA NTRS1988

Large-angle slewing maneuvers for flexible spacecraft

A new class of closed-form solutions for finite-time linear-quadratic optimal control problems is presented. The solutions involve Potter's solution for the differential matrix Riccati equation, which assumes the form of a steady-state plus transient term. Illustrative examples are presented which show that the new solutions are more computationally efficient than alternative solutions based on the state transition matrix. As an application of the closed-form solutions, the neighboring extremal path problem is presented for a spacecraft retargeting maneuver where a perturbed plant with off-nominal boundary conditions now follows a neighboring optimal trajectory. The perturbation feedback approach is further applied to three-dimensional slewing maneuvers of large flexible spacecraft. For this problem, the nominal solution is the optimal three-dimensional rigid body slew. The perturbation feedback then limits the deviations from this nominal solution due to the flexible body effects. The use of frequency shaping in both the nominal and perturbation feedback formulations reduces the excitation of high-frequency unmodeled modes. A modified Kalman filter is presented for estimating the plant states.

Chun, Hon M.↗
NASA NTRS1992

Collapse of composite tubes under end moments

The Reissner (1959) approximate closed-form solution of the collapse of isotropic tubes subject to end moments is presently extended to encompass specially orthotropic material. This closed-form solution is verified via extensive nonlinear FEM analysis of the collapse of long tubes under applied end-moments; the radius/thickness value and composite layup ranges are representative of space-station truss framework designs. The FEM analysis validates the assumption of an inextensional deformation of the cylindrical cross-section, as well as the approximation of the material as specially orthotropic.

Stockwell, Alan E.↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Collapse of composite tubes under end moments

Cylindrical tubes of moderate wall thickness such as those proposed for the original space station truss, may fail due to the gradual collapse of the tube cross section as it distorts under load. Sometimes referred to as the Brazier instability, it is a nonlinear phenomenon. This paper presents an extension of an approximate closed form solution of the collapse of isotropic tubes subject to end moments developed by Reissner in 1959 to include specially orthotropic material. The closed form solution was verified by an extensive nonlinear finite element analysis of the collapse of long tubes under applied end moments for radius to thickness ratios and composite layups in the range proposed for recent space station truss framework designs. The finite element analysis validated the assumption of inextensional deformation of the cylindrical cross section and the approximation of the material as specially orthotropic.

Stockwell, Alan E.↗
NASA NTRS2004

Automatic Generation of Algorithms for the Statistical Analysis of Planetary Nebulae Images

Analyzing data sets collected in experiments or by observations is a Core scientific activity. Typically, experimentd and observational data are &aught with uncertainty, and the analysis is based on a statistical model of the conjectured underlying processes, The large data volumes collected by modern instruments make computer support indispensible for this. Consequently, scientists spend significant amounts of their time with the development and refinement of the data analysis programs. AutoBayes [GF+02, FS03] is a fully automatic synthesis system for generating statistical data analysis programs. Externally, it looks like a compiler: it takes an abstract problem specification and translates it into executable code. Its input is a concise description of a data analysis problem in the form of a statistical model as shown in Figure 1; its output is optimized and fully documented C/C++ code which can be linked dynamically into the Matlab and Octave environments. Internally, however, it is quite different: AutoBayes derives a customized algorithm implementing the given model using a schema-based process, and then further refines and optimizes the algorithm into code. A schema is a parameterized code template with associated semantic constraints which define and restrict the template s applicability. The schema parameters are instantiated in a problem-specific way during synthesis as AutoBayes checks the constraints against the original model or, recursively, against emerging sub-problems. AutoBayes schema library contains problem decomposition operators (which are justified by theorems in a formal logic in the domain of Bayesian networks) as well as machine learning algorithms (e.g., EM, k-Means) and nu- meric optimization methods (e.g., Nelder-Mead simplex, conjugate gradient). AutoBayes augments this schema-based approach by symbolic computation to derive closed-form solutions whenever possible. This is a major advantage over other statistical data analysis systems which use numerical approximations even in cases where closed-form solutions exist. AutoBayes is implemented in Prolog and comprises approximately 75.000 lines of code. In this paper, we take one typical scientific data analysis problem-analyzing planetary nebulae images taken by the Hubble Space Telescope-and show how AutoBayes can be used to automate the implementation of the necessary anal- ysis programs. We initially follow the analysis described by Knuth and Hajian [KHO2] and use AutoBayes to derive code for the published models. We show the details of the code derivation process, including the symbolic computations and automatic integration of library procedures, and compare the results of the automatically generated and manually implemented code. We then go beyond the original analysis and use AutoBayes to derive code for a simple image segmentation procedure based on a mixture model which can be used to automate a manual preproceesing step. Finally, we combine the original approach with the simple segmentation which yields a more detailed analysis. This also demonstrates that AutoBayes makes it easy to combine different aspects of data analysis.

Fischer, Bernd↗
NASA NTRS1960

An Approximate Analytical Method for Studying Atmosphere Entry of Vehicles with Modulated Aerodynamic Forces

The dimensionless, transformed, nonlinear differential equation developed in NASA TR R-11 for describing the approximate motion and heating during entry into planetary atmospheres for constant aerodynamic coefficients and vehicle shape has been modified to include entries during which the aerodynamic coefficients and the vehicle shape are varied. The generality of the application of the original equation to vehicles of arbitrary weight, size, and shape and to arbitrary atmospheres is retained. A closed-form solution for the motion, heating, and the variation of drag loading parameter m/C(D)A has been obtained for the case of constant maximum resultant deceleration during nonlifting entries. This solution requires certain simplifying assumptions which do not compromise the accuracy of the results. The closed-form solution has been used to determine the variation of m/C(D)A required to reduce peak decelerations and to broaden the corridor for nonlifting entry into the earth's atmosphere at escape velocity. The attendant heating penalty is also studied.

Levy, Lionel L., Jr.↗
NASA NTRS1991

Vibration suppression using a constrained rate-feedback-threshold control strategy

A finite time, minimum force rate-feedback-threshold controller is developed to bring a system with or without known external disturbances back into an 'allowable' state bound in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The optimal control is defined by a two-point boundary value problem coupled to a set of definite integral constraints. Quasi-closed form solutions are derived which replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples are provided which demonstrate the threshold control technique and compare the quasi-closed form solutions with numerical and MACSYMA generated exact solutions.

Zimmerman, D. C.↗
NASA NTRS2013

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology. Appendix C -- Finite Element Models Solution Database File, Appendix D -- Benchmark Finite Element Models Solution Database File

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

Allen, Phillip A.↗
NASA NTRS1988

Composite transport wing technology development

The design, fabrication, testing, and analysis of stiffened wing cover panels to assess damage tolerance criteria are discussed. The damage tolerance improvements were demonstrated in a test program using full-sized cover panel subcomponents. The panels utilized a hard skin concept with identical laminates of 44-percent 0-degree, 44-percent plus or minus 45-degree, and 12-percent 90-degree plies in the skins and stiffeners. The panel skins were impacted at midbay between the stiffeners, directly over the stiffener, and over the stiffener flange edge. The stiffener blades were impacted laterally. Impact energy levels of 100 ft-lb and 200 ft-lb were used. NASTRAN finite-element analyses were performed to simulate the nonvisible damage that was detected in the panels by nondestructive inspection. A closed-form solution for generalized loading was developed to evaluate the peel stresses in the bonded structure. Two-dimensional delamination growth analysis was developed using the principle of minimum potential energy in terms of closed-form solution for critical strain. An analysis was conducted to determine the residual compressive stress in the panels after impact damage, and the analytical predictions were verified by compression testing of the damaged panels.

Madan, Ram C.↗
DOE OSTI2021

Verification Testing For Solid-Element Material Models 11-19 in DYNA3D/ParaDyn

This technical report documents the creation and implementation of verification tests for solid-element material models 11 through 19 available in DYNA3D/ParaDyn. The verification tests covered all aspects of each material model, except for the Weibull distribution functionality in material models 15 and 19. General test cases were created to verify the elastic and plastic behavior of the material models. Other additional tests were developed to examine the intricacies of each material model. Each test involved the use a kinematic load case and specification of material parameters necessary to activate corresponding features of the material model. The load cases prescribed the full time history of the kinematic motion for the solid elements, and these loads are independent of the material model or element formulation. When possible, closed form solutions were then derived for each verification test in a continuum setting. The DYNA3D simulations for each test were carried out over a unit time interval, t ϵ [0, 1], and the as implemented DYNA3D response was compared to the closed form solutions evaluated at discrete points in time. A relative error measure was determined for each test to justify the proper implementation of the material model. The relative errors comparing the DYNA3D solution to the analytical solution were, in general, on the order of machine precision except where noted. This signifies the proper implementation of solid-element material models 11-19. In the development of these verification tests, six bugs were found and fixed in the source code. Additionally, this work generated eighteen DYNA3D input decks and answer extraction scripts in the DYNA3D/ParaDyn Software Quality Assurance test suite, which are comprised of a total of 285 solid-element tests. Testing for each material model utilizes two input decks and answer extraction scripts, where one focuses on the linear elastic response and the other examines the inelastic and remaining functionalities of the material model. In total, this work added 285 individual verification test problems in the DYNA3D/ParaDyn test suite.

42 ENGINEERING↗
NASA NTRS1982

Stability of flow over axisymmetric bodies with porous suction strips

Linear triple deck, closed form solutions for mean-flow quantities are developed for axisymmetric incompressible flow past a body with porous strips. The solutions account for upstream influence and are linear superpositions of the flow past the body without suction plus the perturbations due to the suction strips. Flow past the suctionless body is calculated using the Transition Analysis Program System, and a simple linear optimization scheme to determine number, spacing, and mass flow rate through the strips on an axisymmetric body is developed using the linear, triple-deck, closed-form solutions. The theory is demonstrated by predicting optimal strip distributions, and the effect of various adverse pressure-gradient situations on stability is studied.

Nayfeh, A. H.↗
NASA NTRS1998

Influence of Material Distribution on Impact Resistance of Hybrid Composites

Impact events occur in a wide variety of circumstances. A typical example is a bullet impacting a target made of composite material. These impact events produce time-varying loads on a structure that can result in damage. As a first step to understanding the damage resistance issue in composite laminates, an accurate prediction of the transient response during an impact event is necessary. The analysis of dynamic loadings on laminated composite plates has undergone considerable development recently. Rayleigh-Ritz energy method was used to determine the impact response of laminated plates. The impact response of composite plates using shear deformation plate theory was analyzed. In recent work a closed-form solution was obtained for a rectangular plate with four edges simply supported subjected to a center impact load using classical plate theory. The problem was further investigated and the analysis results compared of both classical plate theory and shear deformation theory, and found that classical plate theory predicts very accurate results for the range of small deformations considered. In this study, the influence of cross sectional material distribution on the comparative impact responses of hybrid metal laminates subjected to low and medium velocity impacts is investigated. A simple linear model to evaluate the magnitude of the impact load is proposed first, and it establishes a relation between the impact velocity and the impact force. Then a closed-form solution for impact problem is presented. The results were compared with the finite element analysis results. For an 11 layer-hybrid laminate, the impact response as a function of material distribution in cross-section is presented. With equal areal weight, the effect of the number of laminate layers on the impact resistance is also investigated. Finally, the significance of the presented results is discussed.

Abatan, Ayu↗
NASA NTRS1991

Combined structures-controls optimization of lattice trusses

The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.

Balakrishnan, A. V.↗
NASA NTRS1971

Approximate optimal atmospheric entry trajectories.

For a reentry glider, approximate solutions are found in closed form for the problem of maximizing a function of the terminal velocity, altitude, flight-path angle and heading angle subject to, at most, three terminal nonlinear constraints. The results given here extend the previous results in two important ways: (1) the second-order approximation of the entry dynamics of Loh is used instead of the first-order approximation of Allen and Eggers. This approximation is found to compare extremely well with an exact numerical optimal path. (2) A three-dimensional optimization problem is solved which includes Loh's second-order approximation where the roll angle as well as the lift coefficient are determined subject to constraints which include the terminal heading angle. Furthermore, closed-form solutions can also be obtained subject to in-flight constraints on the lift coefficient and roll angle.

Speyer, J. L.↗