Search NASA⌕ Search

SEARCH · Search NASA

Results for “first principles simulations”

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 181 records · Page 10

Refined Zigzag Theory for Laminated Composite and Sandwich Plates

A refined zigzag theory is presented for laminated-composite and sandwich plates that includes the kinematics of first-order shear deformation theory as its baseline. The theory is variationally consistent and is derived from the virtual work principle. Novel piecewise-linear zigzag functions that provide a more realistic representation of the deformation states of transverse-shear-flexible plates than other similar theories are used. The formulation does not enforce full continuity of the transverse shear stresses across the plate s thickness, yet is robust. Transverse-shear correction factors are not required to yield accurate results. The theory is devoid of the shortcomings inherent in the previous zigzag theories including shear-force inconsistency and difficulties in simulating clamped boundary conditions, which have greatly limited the accuracy of these theories. This new theory requires only C(sup 0)-continuous kinematic approximations and is perfectly suited for developing computationally efficient finite elements. The theory should be useful for obtaining relatively efficient, accurate estimates of structural response needed to design high-performance load-bearing aerospace structures.

Tessler, Alexander↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Aerospace Medicine and Biology: A continuing bibliography with indexes, supplement 143

This supplement to Aerospace Medicine and Biology (NASA SP-7011) lists 251 reports, articles and other documents announced during June 1975 in Scientific and Technical Aerospace Reports (STAR) or in International Aerospace Abstracts (IAA). The first issue of the bibliography was published in July 1964; since that time, monthly supplements have been issued. In its subject coverage, Aerospace Medicine and Biology concentrates on the biological, physiological, and environmental effects to which man is subjected during and following simulated or actual flight in the earth's atmosphere or in interplanetary space. References describing similar effects of biological organisms of lower order are also included. Such related topics as sanitary problems, pharmacology, toxicology, safety and survival, life support systems, exobiology, and personnel factors receive appropriate attention. In general, emphasis is placed on applied research, but references to fundamental studies and theoretical principles related to experimental development also qualify for inclusion.

Source record↗

Use of Soft Computing Technologies For Rocket Engine Control

The problem to be addressed in this paper is to explore how the use of Soft Computing Technologies (SCT) could be employed to further improve overall engine system reliability and performance. Specifically, this will be presented by enhancing rocket engine control and engine health management (EHM) using SCT coupled with conventional control technologies, and sound software engineering practices used in Marshall s Flight Software Group. The principle goals are to improve software management, software development time and maintenance, processor execution, fault tolerance and mitigation, and nonlinear control in power level transitions. The intent is not to discuss any shortcomings of existing engine control and EHM methodologies, but to provide alternative design choices for control, EHM, implementation, performance, and sustaining engineering. The approaches outlined in this paper will require knowledge in the fields of rocket engine propulsion, software engineering for embedded systems, and soft computing technologies (i.e., neural networks, fuzzy logic, and Bayesian belief networks), much of which is presented in this paper. The first targeted demonstration rocket engine platform is the MC-1 (formerly FASTRAC Engine) which is simulated with hardware and software in the Marshall Avionics & Software Testbed laboratory that

Trevino, Luis C.↗

Task-Level Control for a Full Semi-Autonomous Mission: Test Platform Development and Demonstration

Pseudolites can extend the availability of GPS-type positioning systems to a wide range of applications not possible with satellite-only GPS, including indoor and deep-space applications. Conventional GPS pseudolite arrays require that the devices be pre-calibrated through a survey of their locations, typically to sub-centimeter accuracy. This can sometimes be a difficult task, especially in remote or hazardous environments. By using the GPS signals that the pseudolites broadcast, however, it is possible to have the array self-survey its own relative locations, creating a Self-Calibrating Pseudolite Array (SCPA). In order to provide the bi-directional ranging signals between devices necessary for array self-calibration, pseudolite transceivers must be used. The basic principles behind the use of transceivers to create an SCPA were first presented in paper presented to the Institute of Navigation GPS-98 Conference. This paper begins with a brief review of the transceiver architecture and the fundamental direct-ranging algorithm presented in that paper. This is followed by a description of a prototype self-differencing transceiver system that has been constructed, and a presentation of experimental code- and carrier-phase ranging data obtained using that system. A second algorithm is then described which uses these fundamental range measurements between transceiver pairs to self-calibrate a larger stationary array and to provide positioning information for a vehicle moving within that array. Simulation results validating the accuracy and effective convergence of this algorithm are also presented.

Rock, Stephen M.↗

Fitting Vibrational Energy Parameters to Falloff Experiments: CH3 + CH3 + He (0.6-2 Torr 200-298 K)

There are many reported experiments and theoretical calculations for the combination of methyl free radicals diluted in a non-reactive 'deactivator. At high pressure the ethane (C2H6) combination product is the sole product; with decreasing pressure the chemically activated C2H6 will decompose. The falloff in the observed rate coefficient is the result of the competition between collisional stabilization of the chemically activated CA and unimolecular decomposition. The dependence of the rate coefficient on pressure, temperature and collision properties is complex and can not be calculated from first principles. The understanding of this system is not only of fundamental importance but is relevant to the recent detection of methyl free radicals in the atmospheres of Saturn - and Neptune. The temperatures of these outer planet atmospheres are in the 140-200 K region with total pressures (predominately H2 and He) less than 0.2 Torr. Experimentally determined rate coefficients have been reported for this reaction at T = 296-906 K and T = 200-408 K mostly with argon as the deactivator. At T = 200 K only the high pressure rate coefficient has been determined. Complete falloff curves over a wide temperature range (200-1600 K) with a variety of weak collider models used to simulate argon as the deactivator have also been reported by Klippenstein and Harding (KH). More recently we have reported the experimental rate coefficients in the falloff region with helium as the deactivator at 200 and 298 K. In this paper we have used the calculated falloff curves reported by KH for argon to determine the average energy transferred per collision for helium in our recently reported experiments. Collision rates were converted using Lennard Jones parameters; the temperature dependence of this conversion factor is noted. The helium experiments were consistent with a down of approximately 100 cm (exp-1); the temperature dependence was slight. The magnitude of down and its temperature dependence will be compared to other systems with similar substrate complexity and reaction energetics.

Tardy, Dwight C.↗

Application of a First-Principles Anomalous Transport Model for Electrons to Multiple Hall Thrusters and Operating Conditions

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Chaplin, Vernon H.↗

Fault tolerance of artificial neural networks with applications in critical systems

This paper investigates the fault tolerance characteristics of time continuous recurrent artificial neural networks (ANN) that can be used to solve optimization problems. The principle of operations and performance of these networks are first illustrated by using well-known model problems like the traveling salesman problem and the assignment problem. The ANNs are then subjected to 13 simultaneous 'stuck at 1' or 'stuck at 0' faults for network sizes of up to 900 'neurons'. The effects of these faults is demonstrated and the cause for the observed fault tolerance is discussed. An application is presented in which a network performs a critical task for a real-time distributed processing system by generating new task allocations during the reconfiguration of the system. The performance degradation of the ANN under the presence of faults is investigated by large-scale simulations, and the potential benefits of delegating a critical task to a fault tolerant network are discussed.

Protzel, Peter W.↗

The rainout parameterization in a photochemical model

Photochemical models of tropospheric trace species typically simulate the effects of wet removal of soluble species from the atmosphere through a first-order rainout parameterization, which is independent of the species solubility and constant in time. In this work a modified first-order rainout parameterization is proposed that is a function of the species' solubility and takes into account the intermittent nature of atmospheric wet removal. In the case of a highly soluble species, such as HNO3, the inclusion of the effects of intermittent rainout in the parameterization leads to averaged rainout lifetimes considerably higher than ones used in most photochemical models. For illustrative purposes this parameterization is applied to a one-dimensional steady state model with parameterized eddy diffusion, although in principle it can be used in more sophisticated multidimensional models. The model results are compared to analagous calculations obtained with a 3-D tracer transport model.

Giorgi, F.↗

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗

Characterization of Space Environmental Effects on Candidate Solar Sail Material

The National Aeronautics and Space Administration's (NASA) Marshall Space Flight Center (MSFC) is concentrating research into the utilization of photonic materials for spacecraft propulsion. Spacecraft propulsion, using photonic materials, will be achieved using a solar sail. A solar sail operates on the principle that photons, originating from the sun, impart pressure to the sail and therefore provide a source for spacecraft propulsion. The pressure imparted to a solar sail can be increased, up to a factor of two if the sunfacing surface is perfectly reflective. Therefore, these solar sails are generally composed of a highly reflective metallic sun-facing layer, a thin polymeric substrate and occasionally a highly emissive back surface. The Space Environmental Effects Team, at MSFC, is actively characterizing candidate solar sail material to evaluate the thermo-optical and mechanical properties after exposure to radiation environments simulating orbital environments. This paper describes the results of three candidate materials after exposure to a simulated Geosynchronous Transfer Orbit (GTO). This is the first known characterization of solar sail material exposed to space simulated radiation environments. The technique of radiation dose versus material depth profiling was used to determine the orbital equivalent exposure doses. The solar sail exposure procedures and results of the material characterization will be discussed.

Edwards, David↗

Minimizing selective availability error on Topex GPS measurements

GPS measurements made at Topex/Poseidon and the accompanying ground tracking sites will be affected by the selective availability. Although in principle the effects may be removed by differencing between receivers observing the same GPS satellites, this requires accurate synchronization of all receiver clocks. In the case of Topex/Poseidon application, there are two sources of imperfect clock synchronization. The first and larger is due to the constantly drifting clock onboard Topex, which may cause a residual effect as large as 10 cm on Topex carrier phase and 1 m on Topex pseudorange. The second is due to light-time differences between receivers observing the same GPS satellites, which may amount to a few mm error. In this paper a data reduction scheme which incorporates a low-order polynomial interpolation and carrier phase smoothing on pseudorange acquired at Topex and ground receivers is described; a simulation analysis is given demonstrating the effectiveness of the scheme for reducing the GPS S/A effects; and comparison with other schemes is discussed.

Wu, S. C.↗

Investigation of television transmission using adaptive delta modulation principles

The results are presented of a study on the use of the delta modulator as a digital encoder of television signals. The computer simulation of different delta modulators was studied in order to find a satisfactory delta modulator. After finding a suitable delta modulator algorithm via computer simulation, the results were analyzed and then implemented in hardware to study its ability to encode real time motion pictures from an NTSC format television camera. The effects of channel errors on the delta modulated video signal were tested along with several error correction algorithms via computer simulation. A very high speed delta modulator was built (out of ECL logic), incorporating the most promising of the correction schemes, so that it could be tested on real time motion pictures. Delta modulators were investigated which could achieve significant bandwidth reduction without regard to complexity or speed. The first scheme investigated was a real time frame to frame encoding scheme which required the assembly of fourteen, 131,000 bit long shift registers as well as a high speed delta modulator. The other schemes involved the computer simulation of two dimensional delta modulator algorithms.

Schilling, D. L.↗

Kinetics of template-directed pyrophosphate-linked dideoxyguanylate synthesis as a function of 2-MeImpdG and poly(C) concentration: insights into the mechanism

Aqueous solutions of deoxyguanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpdG, yield primarily deoxyguanosine 5'-monophosphate, 5'dGMP, and pyrophosphate-linked dideoxyguanylate, dG5'ppdG, abbreviated G2p (see Chart 1). The initial rate of G2p formation, d[G2p]/dt in M h-1, determined at 23 degrees C, pH 7.8, 1.0 M NaCl and 0.2 M Mg2+ by timed high-performance liquid chromatography (HPLC) analysis, exhibits a second-order dependence on 2-MeImpdG concentration, [G]o, indicating a bimolecular mechanism of dimerization in the range 0.02 M < or = [G]o < or = 0.09 M. In the presence of polycytidylate, poly(C), G2p synthesis is accelerated and oligodeoxyguanylate products are formed by incorporation of 2-MeImpdG molecules. The kinetics of G2p formation as a function of both monomer and polymer concentration, expressed in C equivalents, were also determined under the above conditions and exhibited a complex behavior. Specifically, at a constant [poly(C)], values of d[G2p]/dt typically increased with [G]o with a parabolic upward curvature. At a constant [G]o, values of d[G2p]/dt increase with [poly(C)], but level off at the higher poly(C) concentrations. As [G]o increases this saturation occurs at a higher poly(C) concentration, a result opposite to expectation for a simple complexation of two reacting monomers with the catalyst prior to reaction. Nevertheless, these results are shown to be quantitatively consistent with a template-directed (TD) mechanism of dimerization where poly(C) acts as the template to bind 2-MeImpdG in a cooperative manner and lead, for the first time, to the formulation of principles that govern template-directed chemistry. Analysis of the kinetic data via a proposed TD cooperative model provides association constants for the affinity between polymer and monomer and the intrinsic reactivity of 2-MeImpdG toward pyrophosphate synthesis. To the best of our knowledge, poly(C)/2-MeImpdG is the first system that could serve as a textbook example of a TD reaction under conditions such that the template is fully saturated by monomers and under conditions that it is not.

NASA Discipline Exobiology↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗

Recent Advances in Interstellar Astrochemistry

This presentation will address the possibility that molecules created in the interstellar medium may play a role in the origin and evolution of life on planetary surfaces. The argument will be presented in three main parts. First, I will review the various organic and volatile compounds that are now known or suspected to exist in a variety of circumstellar and interstellar environments. This information is largely the result of the combined applications of observational infrared and radio spectroscopy, laboratory astrophysical simulations, and theoretical astrochemistry. Second, I will devote additional attention to those materials that exhibit chemical properties that are of biological interest. Finally, I will discuss the evidence, largely gathered from the laboratory study of extraterrestrial materials (meteorites and cosmic dust), that interstellar materials, including organics, can survive the transition from the dense cloud phase into forming stellar systems. Once there, some of this material can be delivered largely unaltered to planetary surfaces where it can play key roles in the origin and subsequent evolution of life. Thus, the guiding principle of this presentation will be to attempt to 'put the astro into astrobiology.'

Sandford, Scott A.↗

Lattice Boltzmann and Navier-Stokes Cartesian CFD Approaches for Airframe Noise Predictions

Lattice Boltzmann (LB) and compressible Navier-Stokes (NS) equations based computational fluid dynamics (CFD) approaches are compared for simulating airframe noise. Both LB and NS CFD approaches are implemented within the Launch Ascent and Vehicle Aerodynamics (LAVA) framework. Both schemes utilize the same underlying Cartesian structured mesh paradigm with provision for local adaptive grid refinement and sub-cycling in time. We choose a prototypical massively separated, wake-dominated flow ideally suited for Cartesian-grid based approaches in this study - The partially-dressed, cavity-closed nose landing gear (PDCC-NLG) noise problem from AIAA's Benchmark problems for Airframe Noise Computations (BANC) series of workshops. The relative accuracy and computational efficiency of the two approaches are systematically compared. Detailed comments are made on the potential held by LB to significantly reduce time-to-solution for a desired level of accuracy within the context of modeling airframes noise from first principles.

Airframe Noise Prediction↗