Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

On the Assessment of Acoustic Scattering and Shielding by Time Domain Boundary Integral Equation Solutions

Based on the time domain boundary integral equation formulation of the linear convective wave equation, a computational tool dubbed Time Domain Fast Acoustic Scattering Toolkit (TD-FAST) has recently been under development. The time domain approach has a distinct advantage that the solutions at all frequencies are obtained in a single computation. In this paper, the formulation of the integral equation, as well as its stabilization by the Burton-Miller type reformulation, is extended to cases of a constant mean flow in an arbitrary direction. In addition, a "Source Surface" is also introduced in the formulation that can be employed to encapsulate regions of noise sources and to facilitate coupling with CFD simulations. This is particularly useful for applications where the noise sources are not easily described by analytical source terms. Numerical examples are presented to assess the accuracy of the formulation, including a computation of noise shielding by a thin barrier motivated by recent Historical Baseline F31A31 open rotor noise shielding experiments. Furthermore, spatial resolution requirements of the time domain boundary element method are also assessed using point per wavelength metrics. It is found that, using only constant basis functions and high-order quadrature for surface integration, relative errors of less than 2% may be obtained when the surface spatial resolution is 5 points-per-wavelength (PPW) or 25 points-per-wavelength squared (PPW2).

Hu, Fang Q.↗
NASA NTRS1994

Near-surface ice on Mercury and the Moon: A topographic thermal model

A thermal model that can be easily adapted to craters of arbitrary shape is developed and applied to high-latitude impact craters on Mercury and the Moon, Chao Meng Fu crater at -87.5 deg L on Mercury, an unnamed bowl-shaped crater at 86.7 deg L on Mercury, and Peary crater at 88.6 deg L on the Moon. For an assumed input topography and grid of surface elements, the model computes for each element the irradiation from direct insolation and reflected and emitted radiation from other elements, taking into account shadowing by walls of the crater, partial obscuration of the solar disk near the poles and the diurnal, orbital, and seasonal cycles. Temperatures are computed over the surface grid as functions of depth and time from the surface to a specified depth and over the pertinent astronomical cycles, including the effects of direct and indirect surface irradiation, infrared radiation, heat conduction, and interior heating. Vapor fluxes and ice recession times are computed as functions of ice depth over the surface grid. Temperatures profiles, vapor fluxes, and ice recession times were computed for flat surfaces not associated with craters near the poles of Mercury and the Moon. It was found that water ice could have existed throughout geologic time within the maximum radar detection depth of recent observation of Mercury (J. K. Harmon and M. A. Slade, 1992, Science 258, 640-643) poleward of approximately 87 - 88 deg L on Mercury and poleward of approximately 73 deg L on the Moon. For Chao Meng Fu crater it was found that approximately 40% of the crater floor is permanently shadowed from direct solar insolation, while the remainder of the crater floor is periodically illuminated by a partially obscured Sun. Temperatures at the upper levels of the south wall can slightly exceed 550 K. Surface temperatures in the permanently shadowed region of the crater floor are under approximately 130 K, which could have allowed water ice to exist throughout geologic time within the radar detection depth of recent observation of Mercury. For small bowl-shaped crater on Mercury, it was found that most of the crater is permanently shadowed from direct solar radiation, except for a narrow semicircular band bordering the north rim. However, temperatures in the permanently shadowed region periodically reach a maximum near approximately 315 K due to efficient heating of the small crater by thermal emission and reflection from the small sunlit region, which periodically reaches temperatures exceeding 630 K. Water ice could not have existed throughout geologic time anywhere in this crater within the radar detection depth. For Peary crater on the Moon, the entire crater floor is permanently shadowed from direct solar insolation with maximum temperature under 120 K. The upper level of the north wall periodically reaches a maximum temperature near 310 K. The low temperatures on the crater floor would have allowed water ice to exist near the surface throughout geologic time, provided that the Moon's obliquity was always as low as it is at present.

Salvail, James R.↗
NASA NTRS1993

A local dynamic model for large eddy simulation

The dynamic model is a method for computing the coefficient C in Smagorinsky's model for the subgrid-scale stress tensor as a function of position from the information already contained in the resolved velocity field rather than treating it as an adjustable parameter. A variational formulation of the dynamic model is described that removes the inconsistency associated with taking C out of the filtering operation. This model, however, is still unstable due to the negative eddy-viscosity. Next, three models are presented that are mathematically consistent as well as numerically stable. The first two are applicable to homogeneous flows and flows with at least one homogeneous direction, respectively, and are, in fact, a rigorous derivation of the ad hoc expressions used by previous authors. The third model in this set can be applied to arbitrary flows, and it is stable because the C it predicts is always positive. Finally, a model involving the subgrid-scale kinetic energy is presented which attempts to model backscatter. This last model has some desirable theoretical features. However, even though it gives results in LES that are qualitatively correct, it is outperformed by the simpler constrained variational models. It is suggested that one of the constrained variational models should be used for actual LES while theoretical investigation of the kinetic energy approach should be continued in an effort to improve its predictive power and to understand more about backscatter.

Ghosal, Sandip↗
NASA NTRS1991

An area model for on-chip memories and its application

An area model suitable for comparing data buffers of different organizations and arbitrary sizes is described. The area model considers the supplied bandwidth of a memory cell and includes such buffer overhead as control logic, driver logic, and tag storage. The model gave less than 10 percent error when verified against real caches and register files. It is shown that, comparing caches and register files in terms of area for the same storage capacity, caches generally occupy more area per bit than register files for small caches because the overhead dominates the cache area at these sizes. For larger caches, the smaller storage cells in the cache provide a smaller total cache area per bit than the register set. Studying cache performance (traffic ratio) as a function of area, it is shown that, for small caches, direct-mapped caches perform significantly better than four-way set-associative caches and, for caches of medium areas, both direct-mapped and set-associative caches perform better than fully associative caches.

Mulder, Johannes M.↗
NASA NTRS1994

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗
NASA NTRS1989

The Hole-Count Test Revisited: Effects of Test Specimen Thickness

For historical reasons the hole count, an important performance test for the Analytical Electron Microscope (AEM), is somewhat arbitrary yielding different numbers for different investigators. This was not a problem a decade ago when AEM specimens were often bathed with large fluxes of stray electrons and hard x rays. At that time the presence or absence of a thick Pt second condenser (C2) aperture could be detected by a simple comparison of the x-ray spectrum taken 'somewhere in the hole' with a spectrum collected on a 'typical thickness' of Mo or Ag foil. A high hole count of about 10-20% indicated that the electron column needed modifications; whereas a hole count of 1-2% was accepted for most AEM work. The absolute level of the hole count is a function of test specimen atomic number, overall specimen shape, and thin-foil thickness. In order that equivalent results may be obtained for any AEM in any laboratory in the world, this test must become standardized. The hole-count test we seek must be as simpl and as nonsubjective as the graphite 0.344nm lattice-line-resolution test. This lattice-resolution test spurred manufacturers to improve the image resolution of the TEM significantly in the 1970s and led to the even more stringent resolution tests of today. A similar phenomenon for AEM instruments would be welcome. The hole-count test can also indicate whether the spurious x-ray signal is generated by high-energy continuum x rays (bremsstrahlung) generated in the electron column (high K-line to L-line ratio) or uncollimated electrons passing through or around the C2 aperture (low K/L ratio).

Lyman, C. E.↗
NASA NTRS1972

Apollo 9 multiband photography experiment S065

Fourier analysis was applied to microdensitometer scans of a selected region of one SO65 frame in each of the three black-and-white bands. The approach was unique because a somewhat arbitrary section of the image was used and not limited to available targets as in edge analysis. Comparison of duplicates and calculation of absolute SO65 MTF were done by applying linear systems theory to the spatial spectra of the image scans. It was found that the duplication process was nonlinear and resulted in general amplification of spatial frequency modulation. However, the increase in modulation was offset by a corresponding increase in the granularity of the copies. The amount of increase seemed to be related to the initial granularity, but a direct relationship was not verified. Band-to-band comparison of image quality was achieved in the form of signal-to-noise ratio curves as a function of spatial frequency for each band. From this standpoint the DD band was an order of magnitude better than the other two. These were several factors that restricted the analysis of the Apollo 9 imagery. Among these were the lack of precise sensitometric and optical system data on the high altitude photography. In addition it was determined that there was only one simultaneous pair of high altitude and SO65 frames. Finally, the original SO65 photography was not available for scanning (for obvious reasons), thus eliminating a reference base for granularity and SO65 MTF determination.

Schowengerdt, R. A.↗
NASA NTRS2011

A Self-Stabilizing Synchronization Protocol for Arbitrary Digraphs

This paper presents a self-stabilizing distributed clock synchronization protocol in the absence of faults in the system. It is focused on the distributed clock synchronization of an arbitrary, non-partitioned digraph ranging from fully connected to 1-connected networks of nodes while allowing for differences in the network elements. This protocol does not rely on assumptions about the initial state of the system, other than the presence of at least one node, and no central clock or a centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. There is no theoretical limit on the maximum number of participating nodes. The only constraint on the behavior of the node is that the interactions with other nodes are restricted to defined links and interfaces. This protocol deterministically converges within a time bound that is a linear function of the self-stabilization period. We present an outline of a deductive proof of the correctness of the protocol. A bounded model of the protocol was mechanically verified for a variety of topologies. Results of the mechanical proof of the correctness of the protocol are provided. The model checking results have verified the correctness of the protocol as they apply to the networks with unidirectional and bidirectional links. In addition, the results confirm the claims of determinism and linear convergence. As a result, we conjecture that the protocol solves the general case of this problem. We also present several variations of the protocol and discuss that this synchronization protocol is indeed an emergent system.

Malekpour, Mahyar R.↗
NASA NTRS1999

NASA Lewis Launch Collision Probability Model Developed and Analyzed

There are nearly 10,000 tracked objects orbiting the earth. These objects encompass manned objects, active and decommissioned satellites, spent rocket bodies, and debris. They range from a few centimeters across to the size of the MIR space station. Anytime a new satellite is launched, the launch vehicle with its payload attached passes through an area of space in which these objects orbit. Although the population density of these objects is low, there always is a small but finite probability of collision between the launch vehicle and one or more of these space objects. Even though the probability of collision is very low, for some payloads even this small risk is unacceptable. To mitigate the small risk of collision associated with launching at an arbitrary time within the daily launch window, NASA performs a prelaunch mission assurance Collision Avoidance Analysis (or COLA). For the COLA of the Cassini spacecraft, the NASA Lewis Research Center conducted an in-house development and analysis of a model for launch collision probability. The model allows a minimum clearance criteria to be used with the COLA analysis to ensure an acceptably low probability of collision. If, for any given liftoff time, the nominal launch vehicle trajectory would pass a space object with less than the minimum required clearance, launch would not be attempted at that time. The model assumes that the nominal positions of the orbiting objects and of the launch vehicle can be predicted as a function of time, and therefore, that any tracked object that comes within close proximity of the launch vehicle can be identified. For any such pair, these nominal positions can be used to calculate a nominal miss distance. The actual miss distances may differ substantially from the nominal miss distance, due, in part, to the statistical uncertainty of the knowledge of the objects positions. The model further assumes that these position uncertainties can be described with position covariance matrices. With these, and some additional simplifying assumptions, a closed-form solution for the probability of collision is obtained. This solution provides clear insights into how each of the independent parameters affects the probability of collision. It shows that for a given maximum probability of collision and prior knowledge of the objects position uncertainties and sizes, only knowledge of the nominal closest approach distance is required to make the launch/no launch decision. The model was completed and used for the mission assurance COLA analysis for the Cassini spacecraft, which was launched on a Titan IV/Centaur rocket on October 15, 1997. Although the model was specifically developed for the Cassini mission, it is clearly applicable for other launches as well. The effect of COLA closures on the launch window is shown. The bar represents the entire 140-min launch window on October 15, 1997; the blackened areas represent the loss of launch opportunities due to the potential for a collision with an orbiting object.

Bollenbacher, Gary↗
NASA NTRS2005

Training Software in Artificial-Intelligence Computing Techniques

The Artificial Intelligence (AI) Toolkit is a computer program for training scientists, engineers, and university students in three soft-computing techniques (fuzzy logic, neural networks, and genetic algorithms) used in artificial-intelligence applications. The program promotes an easily understandable tutorial interface, including an interactive graphical component through which the user can gain hands-on experience in soft-computing techniques applied to realistic example problems. The tutorial provides step-by-step instructions on the workings of soft-computing technology, whereas the hands-on examples allow interaction and reinforcement of the techniques explained throughout the tutorial. In the fuzzy-logic example, a user can interact with a robot and an obstacle course to verify how fuzzy logic is used to command a rover traverse from an arbitrary start to the goal location. For the genetic algorithm example, the problem is to determine the minimum-length path for visiting a user-chosen set of planets in the solar system. For the neural-network example, the problem is to decide, on the basis of input data on physical characteristics, whether a person is a man, woman, or child. The AI Toolkit is compatible with the Windows 95,98, ME, NT 4.0, 2000, and XP operating systems. A computer having a processor speed of at least 300 MHz, and random-access memory of at least 56MB is recommended for optimal performance. The program can be run on a slower computer having less memory, but some functions may not be executed properly.

Howard, Ayanna↗
NASA NTRS2004

A Novel, Real-Valued Genetic Algorithm for Optimizing Radar Absorbing Materials

A novel, real-valued Genetic Algorithm (GA) was designed and implemented to minimize the reflectivity and/or transmissivity of an arbitrary number of homogeneous, lossy dielectric or magnetic layers of arbitrary thickness positioned at either the center of an infinitely long rectangular waveguide, or adjacent to the perfectly conducting backplate of a semi-infinite, shorted-out rectangular waveguide. Evolutionary processes extract the optimal physioelectric constants falling within specified constraints which minimize reflection and/or transmission over the frequency band of interest. This GA extracted the unphysical dielectric and magnetic constants of three layers of fictitious material placed adjacent to the conducting backplate of a shorted-out waveguide such that the reflectivity of the configuration was 55 dB or less over the entire X-band. Examples of the optimization of realistic multi-layer absorbers are also presented. Although typical Genetic Algorithms require populations of many thousands in order to function properly and obtain correct results, verified correct results were obtained for all test cases using this GA with a population of only four.

Hall, John Michael↗
NASA NTRS2010

Reactive Collision Avoidance Algorithm

The reactive collision avoidance (RCA) algorithm allows a spacecraft to find a fuel-optimal trajectory for avoiding an arbitrary number of colliding spacecraft in real time while accounting for acceleration limits. In addition to spacecraft, the technology can be used for vehicles that can accelerate in any direction, such as helicopters and submersibles. In contrast to existing, passive algorithms that simultaneously design trajectories for a cluster of vehicles working to achieve a common goal, RCA is implemented onboard spacecraft only when an imminent collision is detected, and then plans a collision avoidance maneuver for only that host vehicle, thus preventing a collision in an off-nominal situation for which passive algorithms cannot. An example scenario for such a situation might be when a spacecraft in the cluster is approaching another one, but enters safe mode and begins to drift. Functionally, the RCA detects colliding spacecraft, plans an evasion trajectory by solving the Evasion Trajectory Problem (ETP), and then recovers after the collision is avoided. A direct optimization approach was used to develop the algorithm so it can run in real time. In this innovation, a parameterized class of avoidance trajectories is specified, and then the optimal trajectory is found by searching over the parameters. The class of trajectories is selected as bang-off-bang as motivated by optimal control theory. That is, an avoiding spacecraft first applies full acceleration in a constant direction, then coasts, and finally applies full acceleration to stop. The parameter optimization problem can be solved offline and stored as a look-up table of values. Using a look-up table allows the algorithm to run in real time. Given a colliding spacecraft, the properties of the collision geometry serve as indices of the look-up table that gives the optimal trajectory. For multiple colliding spacecraft, the set of trajectories that avoid all spacecraft is rapidly searched on-line. The optimal avoidance trajectory is implemented as a receding-horizon model predictive control law. Therefore, at each time step, the optimal avoidance trajectory is found and the first time step of its acceleration is applied. At the next time step of the control computer, the problem is re-solved and the new first time step is again applied. This continual updating allows the RCA algorithm to adapt to a colliding spacecraft that is making erratic course changes.

Scharf, Daniel↗
NASA NTRS1991

Estimation of time averages from irregularly spaced observations - With application to coastal zone color scanner estimates of chlorophyll concentration

The sampling error of an arbitrary linear estimate of a time-averaged quantity constructed from a time series of irregularly spaced observations at a fixed located is quantified through a formalism. The method is applied to satellite observations of chlorophyll from the coastal zone color scanner. The two specific linear estimates under consideration are the composite average formed from the simple average of all observations within the averaging period and the optimal estimate formed by minimizing the mean squared error of the temporal average based on all the observations in the time series. The resulting suboptimal estimates are shown to be more accurate than composite averages. Suboptimal estimates are also found to be nearly as accurate as optimal estimates using the correct signal and measurement error variances and correlation functions for realistic ranges of these parameters, which makes it a viable practical alternative to the composite average method generally employed at present.

Chelton, Dudley B.↗
NASA NTRS2010

Ka-Band Waveguide Two-Way Hybrid Combiner for MMIC Amplifiers

The design, simulation, and characterization of a novel Ka-band (32.05 0.25 GHz) rectangular waveguide two-way branch-line hybrid unequal power combiner (with port impedances matched to that of a standard WR-28 waveguide) has been created to combine input signals, which are in phase and with an amplitude ratio of two. The measured return loss and isolation of the branch-line hybrid are better than 22 and 27 dB, respectively. The measured combining efficiency is 92.9 percent at the center frequency of 32.05 GHz. This circuit is efficacious in combining the unequal output power from two Ka-band GaAs pseudomorphic high electron mobility transistor (pHEMT) monolithic microwave integrated circuit (MMIC) power amplifiers (PAs) with high efficiency. The component parts include the branch-line hybrid-based power combiner and the MMIC-based PAs. A two-way branch-line hybrid is a four-port device with all ports matched; power entering port 1 is divided in phase, and into the ratio 2:1 between ports 3 and 4. No power is coupled to port 2. MMICs are a type of integrated circuit fabricated on GaAs that operates at microwave frequencies, and performs the function of signal amplification. The power combiner is designed to operate over the frequency band of 31.8 to 32.3 GHz, which is NASA's deep space frequency band. The power combiner would have an output return loss better than 20 dB. Isolation between the output port and the isolated port is greater than 25 dB. Isolation between the two input ports is greater than 25 dB. The combining efficiency would be greater than 90 percent when the ratio of the two input power levels is two. The power combiner is machined from aluminum with E-plane split-block arrangement, and has excellent reliability. The flexibility of this design allows the combiner to be customized for combining the power from MMIC PAs with an arbitrary power output ratio. In addition, it allows combining a low-power GaAs MMIC with a high-power GaN MMIC. The arbitrary port impedance allows matching the output impedance of the MMIC PA directly to the waveguide impedance without transitioning first into a transmission line with characteristic impedance of 50 ohms. Thus, by eliminating the losses associated with a transition, the overall SSPA efficiency is enhanced. For reducing the cost and weight when required in very large quantities, such as in the beam-forming networks of phased-array antenna systems, the combiner can be manufactured using metal-plated plastic. Two hybrid unequal power combiners can be cascaded to realize a non-binary combiner (for e.g., a three-way) and can be synergistically optimized for low VSWR (voltage standing wave ratio), low insertion loss, high isolation, and wide bandwidth using commercial off-the-shelf electromagnetic software design tools.

Simons, Rainee N.↗
NASA NTRS2013

Tool for Automated Retrieval of Generic Event Tracks (TARGET)

Methods have been developed to identify and track tornado-producing mesoscale convective systems (MCSs) automatically over the continental United States, in order to facilitate systematic studies of these powerful and often destructive events. Several data sources were combined to ensure event identification accuracy. Records of watches and warnings issued by National Weather Service (NWS), and tornado locations and tracks from the Tornado History Project (THP) were used to locate MCSs in high-resolution precipitation observations and GOES infrared (11-micron) Rapid Scan Operation (RSO) imagery. Thresholds are then applied to the latter two data sets to define MCS events and track their developments. MCSs produce a broad range of severe convective weather events that are significantly affecting the living conditions of the populations exposed to them. Understanding how MCSs grow and develop could help scientists improve their weather prediction models, and also provide tools to decision-makers whose goals are to protect populations and their property. Associating storm cells across frames of remotely sensed images poses a difficult problem because storms evolve, split, and merge. Any storm-tracking method should include the following processes: storm identification, storm tracking, and quantification of storm intensity and activity. The spatiotemporal coordinates of the tracks will enable researchers to obtain other coincident observations to conduct more thorough studies of these events. In addition to their tracked locations, their areal extents, precipitation intensities, and accumulations all as functions of their evolutions in time were also obtained and recorded for these events. All parameters so derived can be catalogued into a moving object database (MODB) for custom queries. The purpose of this software is to provide a generalized, cross-platform, pluggable tool for identifying events within a set of scientific data based upon specified criteria with the possibility of storing identified events into a searchable database. The core of the application uses an implementation of the connected component labeling (CCL) algorithm to identify areas of interest, then uses a set of criteria to establish spatial and temporal relationships between identified components. The CCL algorithm is used for identifying objects within images for computer vision. This application applies it to scientific data sets using arbitrary criteria. The most novel concept was applying a generalized CCL implementation to scientific data sets for establishing events both spatially and temporally. The combination of several existing concepts (pluggable components, generalized CCL algorithm, etc.) into one application is also novel. In addition, how the system is designed, i.e., its extensibility with pluggable components, and its configurability with a simple configuration file, is innovative. This allows the system to be applied to new scenarios with ease.

Clune, Thomas↗
NASA NTRS1971

Analysis of solar ultraviolet lines

The formation of the strongest ultra-violet emission lines of Mg II, O I, C II, and C III in the solar atmosphere is studied in detail. The equations of statistical equilibrium and radiative transfer for each ion are solved using a general computer program that is capable of solving non-LTE line-formation problems for arbitrary atmospheric and atomic models. Interpreting the results in terms of the structure of the solar atmosphere, it is concluded that the HSRA atmosphere has a temperature too low by about 500 K near h = 1100 km and that a temperature plateau with T sub e approximately = 18,000 K and width close to 60 km exists in the upper chromosphere. The structure of the solar atmosphere in the range 20,000 to 100,000 K and the effects of microturbulence on the formation of lines are also investigated. Approximate analytic line-formation problems are solved, and more exact solutions are derived later. An attempt is made to make the best possible fit to the Ca II K line center-to-limb profiles with a one-component atmosphere, with an assumed source function and microturbulent velocity.

Chipman, E.↗
NASA NTRS2008

Scripting Module for the Satellite Orbit Analysis Program (SOAP)

This add-on module to the SOAP software can perform changes to simulation objects based on the occurrence of specific conditions. This allows the software to encompass simulation response of scheduled or physical events. Users can manipulate objects in the simulation environment under programmatic control. Inputs to the scripting module are Actions, Conditions, and the Script. Actions are arbitrary modifications to constructs such as Platform Objects (i.e. satellites), Sensor Objects (representing instruments or communication links), or Analysis Objects (user-defined logical or numeric variables). Examples of actions include changes to a satellite orbit ( v), changing a sensor-pointing direction, and the manipulation of a numerical expression. Conditions represent the circumstances under which Actions are performed and can be couched in If-Then-Else logic, like performing v at specific times or adding to the spacecraft power only when it is being illuminated by the Sun. The SOAP script represents the entire set of conditions being considered over a specific time interval. The output of the scripting module is a series of events, which are changes to objects at specific times. As the SOAP simulation clock runs forward, the scheduled events are performed. If the user sets the clock back in time, the events within that interval are automatically undone. This script offers an interface for defining scripts where the user does not have to remember the vocabulary of various keywords. Actions can be captured by employing the same user interface that is used to define the objects themselves. Conditions can be set to invoke Actions by selecting them from pull-down lists. Users define the script by selecting from the pool of defined conditions. Many space systems have to react to arbitrary events that can occur from scheduling or from the environment. For example, an instrument may cease to draw power when the area that it is tasked to observe is not in view. The contingency of the planetary body blocking the line of sight is a condition upon which the power being drawn is set to zero. It remains at zero until the observation objective is again in view. Computing the total power drawn by the instrument over a period of days or weeks can now take such factors into consideration. What makes the architecture especially powerful is that the scripting module can look ahead and behind in simulation time, and this temporal versatility can be leveraged in displays such as x-y plots. For example, a plot of a satellite s altitude as a function of time can take changes to the orbit into account.

Carnright, Robert↗
NASA NTRS2017

A Compilation of MATLAB Scripts and Functions for MACGMC Analyses

The primary aim of the current effort is to provide scripts that automate many of the repetitive pre- and post-processing tasks associated with composite materials analyses using the Micromechanics Analysis Code with the Generalized Method of Cells. This document consists of a compilation of hundreds of scripts that were developed in MATLAB (The Mathworks, Inc., Natick, MA) programming language and consolidated into 16 MATLAB functions. (MACGMC). MACGMC is a composite material and laminate analysis software code developed at NASA Glenn Research Center. The software package has been built around the generalized method of cells (GMC) family of micromechanics theories. The computer code is developed with a user-friendly framework, along with a library of local inelastic, damage, and failure models. Further, application of simulated thermo-mechanical loading, generation of output results, and selection of architectures to represent the composite material have been automated to increase the user friendliness, as well as to make it more robust in terms of input preparation and code execution. Finally, classical lamination theory has been implemented within the software, wherein GMC is used to model the composite material response of each ply. Thus, the full range of GMC composite material capabilities is available for analysis of arbitrary laminate configurations as well. The pre-processing tasks include generation of a multitude of different repeating unit cells (RUCs) for CMCs and PMCs, visualization of RUCs from MACGMC input and output files and generation of the RUC section of a MACGMC input file. The post-processing tasks include visualization of the predicted composite response, such as local stress and strain contours, damage initiation and progression, stress-strain behavior, and fatigue response. In addition to the above, several miscellaneous scripts have been developed that can be used to perform repeated Monte-Carlo simulations to enable probabilistic simulations with minimal manual intervention. This document is formatted to provide MATLAB source files and descriptions of how to utilize them. It is assumed that the user has a basic understanding of how MATLAB scripts work and some MATLAB programming experience.

Materials Engineering↗