SEARCH · Search NASA

Results for “MATHEMATICAL EVOLUTION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2013

On the Theory and Numerical Simulation of Cohesive Crack Propagation with Application to Fiber-Reinforced Composites

The phenomenon of crack propagation is among the predominant modes of failure in many natural and engineering structures, often leading to severe loss of structural integrity and catastrophic failure. Thus, the ability to understand and a priori simulate the evolution of this failure mode has been one of the cornerstones of applied mechanics and structural engineering and is broadly referred to as "fracture mechanics." The work reported herein focuses on extending this understanding, in the context of through-thickness crack propagation in cohesive materials, through the development of a continuum-level multiscale numerical framework, which represents cracks as displacement discontinuities across a surface of zero measure. This report presents the relevant theory, mathematical framework, numerical modeling, and experimental investigations of through-thickness crack propagation in fiber-reinforced composites using the Variational Multiscale Cohesive Method (VMCM) developed by the authors.

Rudraraju, Siva Shankar↗

NASA NTRS1982

Theoretical studies in isoelectric focusing

A mathematical model has been developed which describes the steady state in an isoelectric focusing (IEF) system with ampholytes or monovalent buffers. The model is based on the fundamental equations describing the component dissociation equilibria, mass transport due to diffusion and electromigration, electroneutrality, and the conservation of charge. The validity and usefulness of the model has been confirmed by using it to formulate buffer systems in actual laboratory experiments. The model has been recently extended to include the evolution of transient states not only in IEF but also in other modes of electrophoresis.

Mosher, R. A.↗

NASA NTRS1979

The mean age of mantle and crustal reservoirs

An analysis of two transport models for trace elements in mantle and crust evolution and related to the observed abundance patterns. In model I, continents are derived by melt extraction over the earth history from undepleted mantle, and the residue forms a depleted mantle which is the current source of mid-ocean ridge basalts. In model II, new additions to continents are derived from a mantle reservoir 2, which becomes more depleted by repeated extraction of melts. The isotopic composition and concentrations of trace elements are shown to reduce to simple mathematical expressions which permit calculations of basic evolutionary parameters. The mean age of the crust mass and isotopic data for the continental crust and the mantle are discussed, concluding with a consideration of their difference in compositions of newly derived crust.

Jacobsen, S. B.↗

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗

NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗

NASA NTRS1986

SRB/SLEEC (Solid Rocket Booster/Shingle Lap Extendible Exit Cone) feasibility study, volume 2. Appendix A: Design study for a SLEEC actuation system

The results are presented of a design feasibility study of a self-contained (powered) actuation system for a Shingle Lap Extendible Exit Cone (SLEEC) for Transportation System (STS). The evolution of the SLEEC actuation system design is reviewed, the final design concept is summarized, and the results of the detailed study of the final concept of the actuation system are treated. A conservative design using proven mechanical components was established as a major program priority. The final mechanical design has a very low development risk since the components, which consist of ballscrews, gearing, flexible shaft drives, and aircraft cables, have extensive aerospace applications and a history of proven reliability. The mathematical model studies have shown that little or no power is required to deploy the SLEEC actuation system because acceleration forces and internal pressure from the rocket plume provide the required energies. A speed control brake is incorporated in the design in order to control the rate of deployment.

Thompson, D. S.↗

NASA NTRS2005

Coevolutionary Free Lunches

Recent work on the mathematical foundations of optimization has begun to uncover its rich structure. In particular, the "No Free Lunch" (NFL) theorems state that any two algorithms are equivalent when their performance is averaged across all possible problems. This highlights the need for exploiting problem-specific knowledge to achieve better than random performance. In this paper we present a general framework covering more search scenarios. In addition to the optimization scenarios addressed in the NFL results, this framework covers multi-armed bandit problems and evolution of multiple co-evolving players. As a particular instance of the latter, it covers "self-play" problems. In these problems the set of players work together to produce a champion, who then engages one or more antagonists in a subsequent multi-player game. In contrast to the traditional optimization case where the NFL results hold, we show that in self-play there are free lunches: in coevolution some algorithms have better performance than other algorithms, averaged across all possible problems. We consider the implications of these results to biology where there is no champion.

Wolpert, David H.↗

NASA NTRS2006

Stochastic Representation of Chaos Using Terminal Attractors

A nonlinear version of the Liouville equation based on terminal attractors is part of a mathematical formalism for describing postinstability motions of dynamical systems characterized by exponential divergences of trajectories leading to chaos (including turbulence as a form of chaos). The formalism can be applied to both conservative systems (e.g., multibody systems in celestial mechanics) and dissipative systems (e.g., viscous fluids). The development of the present formalism was undertaken in an effort to remove positive Lyapunov exponents. The means chosen to accomplish this is coupling of the governing dynamical equations with the corresponding Liouville equation that describes the evolution of the flow of error probability. The underlying idea is to suppress the divergences of different trajectories that correspond to different initial conditions, without affecting a target trajectory, which is one that starts with prescribed initial conditions.

Zak, Michail↗

NASA NTRS1993

Kalman filtering applied to real-time monitoring of apogee maneuvers

Part of the Space Mathematics Division in CNES, the Flight Dynamics Center provides attitude and orbit determinations and maneuvers during the Launch and Early Operation Phase (LEOP) of geostationary satellites. Orbit determination is based on a Kalman filter method; when the 2 GHz CNES/NASA network is used, Doppler measurements are available and allow orbit determination during the apogee maneuvers. This method was used for TELE-X and TDF 2 LEOP (3-axis controlled satellites) and also for TELECOM 2 and HISPASAT (spun satellites): it enables us to follow the evolution of the maneuver and gives out a quite accurate estimation of the reached orbit. In this paper, we briefly describe the dynamic models of the orbit evolution in both cases, '3-axis' and 'inertial' thrust. Then, we present the results obtained for each case. Afterwards, we present some cases to show the robustness of the filter.

Deboer, Frederic↗

NASA NTRS2003

Extension of Liouville Formalism to Postinstability Dynamics

A mathematical formalism has been developed for predicting the postinstability motions of a dynamic system governed by a system of nonlinear equations and subject to initial conditions. Previously, there was no general method for prediction and mathematical modeling of postinstability behaviors (e.g., chaos and turbulence) in such a system. The formalism of nonlinear dynamics does not afford means to discriminate between stable and unstable motions: an additional stability analysis is necessary for such discrimination. However, an additional stability analysis does not suggest any modifications of a mathematical model that would enable the model to describe postinstability motions efficiently. The most important type of instability that necessitates a postinstability description is associated with positive Lyapunov exponents. Such an instability leads to exponential growth of small errors in initial conditions or, equivalently, exponential divergence of neighboring trajectories. The development of the present formalism was undertaken in an effort to remove positive Lyapunov exponents. The means chosen to accomplish this is coupling of the governing dynamical equations with the corresponding Liouville equation that describes the evolution of the flow of error probability. The underlying idea is to suppress the divergences of different trajectories that correspond to different initial conditions, without affecting a target trajectory, which is one that starts with prescribed initial conditions.

Zak, Michail↗

NASA NTRS1990

A simple analytic model for the evolution of captured galactic disks

The general analytic solution developed by Steinman-Cameron and Durisen and published in 1988 for the evolution of dissipative nonplanar disks is applied to captured galactic disks in model galaxies with nonspherical, scale-free, logarithmic gravitational potentials. Such potentials produce flat rotation curves, similar to those seen in real galaxies. In this case, the analytic solution yields a self-similar structure for the warps and twists that develop while the disks is settling. Being scale-free in a simple, analytic form, this solution is completely defined by only a few dimensionless fitting parameters. As a result, it can be utilized as a mathematical tool to fit settling disks in real galaxies. The minimum time it takes for a disk to settle into a steady state orientation is also a scale-free quantity when expressed in units of the precession period or the orbit period. For realistic parameters, settling times are on the order of one-half to two periods. The use of the time-dependent structure of settling disks as a probe of the three-dimensional mass distribution of the host galaxies, including dark halos, is discussed.

Steiman-Cameron, Thomas Y.↗

NASA NTRS1990

Networks for image acquisition, processing and display

The human visual system comprises layers of networks which sample, process, and code images. Understanding these networks is a valuable means of understanding human vision and of designing autonomous vision systems based on network processing. Ames Research Center has an ongoing program to develop computational models of such networks. The models predict human performance in detection of targets and in discrimination of displayed information. In addition, the models are artificial vision systems sharing properties with biological vision that has been tuned by evolution for high performance. Properties include variable density sampling, noise immunity, multi-resolution coding, and fault-tolerance. The research stresses analysis of noise in visual networks, including sampling, photon, and processing unit noises. Specific accomplishments include: models of sampling array growth with variable density and irregularity comparable to that of the retinal cone mosaic; noise models of networks with signal-dependent and independent noise; models of network connection development for preserving spatial registration and interpolation; multi-resolution encoding models based on hexagonal arrays (HOP transform); and mathematical procedures for simplifying analysis of large networks.

Ahumada, Albert J., Jr.↗

NASA NTRS1999

A Quantitative Study of Oxygen as a Metabolic Regulator

An acute reduction in oxygen (O2) delivery to a tissue is generally associated with a decrease in phosphocreatine, increases in ADP, NADH/NAD, and inorganic phosphate, increased rates of glycolysis and lactate production, and reduced rates of pyruvate and fatty acid oxidation. However, given the complexity of the human bioenergetic system and its components, it is difficult to determine quantitatively how cellular metabolic processes interact to maintain ATP homeostasis during stress (e.g., hypoxia, ischemia, and exercise). Of special interest is the determination of mechanisms relating tissue oxygenation to observed metabolic responses at the tissue, organ, and whole body levels and the quantification of how changes in tissue O2 availability affect the pathways of ATP synthesis and the metabolites that control these pathways. In this study, we extend a previously developed mathematical model of human bioenergetics to provide a physicochemical framework that permits quantitative understanding of O2 as a metabolic regulator. Specifically, the enhancement permits studying the effects of variations in tissue oxygenation and in parameters controlling the rate of cellular respiration on glycolysis, lactate production, and pyruvate oxidation. The whole body is described as a bioenergetic system consisting of metabolically distinct tissue/organ subsystems that exchange materials with the blood. In order to study the dynamic response of each subsystem to stimuli, we solve the ordinary differential equations describing the temporal evolution of metabolite levels, given the initial concentrations. The solver used in the present study is the packaged code LSODE, as implemented in the NASA Lewis kinetics and sensitivity analysis code, LSENS. A major advantage of LSENS is the efficient procedures supporting systematic sensitivity analysis, which provides the basic methods for studying parameter sensitivities (i.e., changes in model behavior due to parameter variation). Sensitivity analysis establishes relationships between model predictions and problem parameters (i.e., initial concentrations, rate coefficients, etc). It helps determine the effects of uncertainties or changes in these input parameters on the predictions, which ultimately are compared with experimental observations in order to validate the model. Sensitivity analysis can identify parameters that must be determined accurately because of their large effect on the model predictions and parameters that need not be known with great precision because they have little or no effect on the solution. This capability may prove to be important in optimizing the design of experiments, thereby reducing the use of animals. This approach can be applied to study the metabolic effects of reduced oxygen delivery to cardiac muscle due to local myocardial ischemia and the effects of acute hypoxia on brain metabolism. Other important applications of sensitivity analysis include identification of quantitatively relevant pathways and biochemical species within an overall mechanism, when examining the effects of a genetic anomaly or pathological state on energetic system components and whole system behavior.

Radhakrishnan, Krishnan↗

NASA NTRS1972

Apollo experience report: S-band system signal design and analysis

A description is given of the Apollo communications-system engineering-analysis effort that ensured the adequacy, performance, and interface compatibility of the unified S-band system elements for a successful lunar-landing mission. The evolution and conceptual design of the unified S-band system are briefly reviewed from a historical viewpoint. A comprehensive discussion of the unified S-band elements includes the salient design features of the system and serves as a basis for a better understanding of the design decisions and analyses. The significant design decisions concerning the Apollo communications-system signal design are discussed providing an insight into the role of systems analysis in arriving at the current configuration of the Apollo communications system. Analyses are presented concerning performance estimation (mathematical-model development through real-time mission support) and system deficiencies, modifications, and improvements.

Rosenberg, H. R.↗

NASA NTRS1994

The art of mapping the magnetosphere

A comprehensive review is presented of the mathematical models used to represent magnetic fields in the Earth's magnetosphere, of the way existing data-based models use these methods and of the associated problems and concepts. The magnetic field has five main components: the internal field, the magnetopause, the ring current, the tail and Birkeland currents. Methods of representing separately each of these are discussed, as is the deformation of magnetic fields; Appendix B traces the connection between deformations and the Cauchy integral. A summary section lists the uses of data-based models and their likely future evolution, and Appendix A supplements the text with a set of problems.

Stern, David P.↗

NASA NTRS1992

A comparison of the reduced and approximate systems for the time dependent computation of the polar wind and multiconstituent stellar winds

The paper derives the 'reduced' system of equations commonly used to describe the time evolution of the polar wind and multiconstituent stellar winds from the equations for a multispecies plasma with known temperature profiles by assuming that the electron thermal speed approaches infinity. The reduced system is proved to have unbounded growth near the sonic point of the protons for many of the standard parameter cases. For the same parameter cases, the unmodified system exhibits growth in some of the Fourier modes, but this growth is bounded. An alternate system (the 'approximate' system) in which the electron thermal speed is slowed down is introduced. The approximate system retains the mathematical behavior of the unmodified system and can be shown to accurately describe the smooth solutions of the unmodified system. Other advantages of the approximate system over the reduced system are discussed.

Browning, G. L.↗

NASA NTRS1975

Interactive computer aided technology, evolution in the design/manufacturing process

A powerful computer-operated three dimensional graphic system and associated auxiliary computer equipment used in advanced design, production design, and manufacturing was described. This system has made these activities more productive than when using older and more conventional methods to design and build aerospace vehicles. With the use of this graphic system, designers are now able to define parts using a wide variety of geometric entities, define parts as fully surface 3-dimensional models as well as "wire-frame" models. Once geometrically defined, the designer is able to take section cuts of the surfaced model and automatically determine all of the section properties of the planar cut, lightpen detect all of the surface patches and automatically determine the volume and weight of the part. Further, his designs are defined mathematically at a degree of accuracy never before achievable.

English, C. H.↗

NASA NTRS1987

A global time-dependent model of thunderstorm electricity. I - Mathematical properties of the physical and numerical models

A time-dependent model is introduced that can be used to simulate the interaction of a thunderstorm with its global electrical environment. The model solves the continuity equation of the Maxwell current, which is assumed to be composed of the conduction, displacement, and source currents. Boundary conditions which can be used in conjunction with the continuity equation to form a well-posed initial-boundary value problem are determined. Properties of various components of solutions of the initial-boundary value problem are analytically determined. The results indicate that the problem has two time scales, one determined by the background electrical conductivity and the other by the time variation of the source function. A numerical method for obtaining quantitative results is introduced, and its properties are studied. Some simulation results on the evolution of the displacement and conduction currents during the electrification of a storm are presented.

Browning, G. L.↗