Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer simulation model development”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,729 records · Page 96

Diurnal Forcing of Planetary Atmospheres

This project investigated and developed models which will assist in developing an understanding of the diurnal forcings in the atmospheres of Mars, Venus, Earth and Jupiter. The differences in the temperatures during the day, and night, generate very deep diurnal boundary layers and strong slope-related winds. A general circulation model (GCM) was adapted for computation of Martian atmospheric dynamics, to simulate the global transport of these atmospheric tracers. The model was applied to the study of the Martian seasonal water cycle. The principles have been applied to venus. The region of study is the upper cloud level, and above rather than the area close to the surface. The model was also applied in efforts to study the evolution of volcanic clouds in the Earth's stratosphere. In particular this was used to model the dispersion of sulfuric acid particles in the atmosphere. A brief study of the application of the tidal theory to the planet Jupiter has been undertaken. The goal of the study is to deduce the source of the tidal dissipation which is thought to have led to the orbital resonances amongst the Galilean satellites and considerable heating of the interiors of Io, Europa, and Ganymede.

Houben, Howard C.↗

Testing and Analysis of Composite Skin/Stringer Debonding Under Multi-Axial Loading

A consistent step-wise approach is presented to investigate the damage mechanism in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions. The approach uses experiments to detect the failure mechanism, computational stress analysis to determine the location of first matrix cracking and computational fracture mechanics to investigate the potential for delamination growth. In a first step, tests were performed on specimens, which consisted of a tapered composite flange, representing a stringer or frame, bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending to evaluate the debonding mechanisms between the skin and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both in-plane tension and out-of-plane bending loads simultaneously. Specimen edges were examined on the microscope to document the damage occurrence and to identify typical damage patterns. For all three load cases, observed failure initiated in the flange, near the flange tip, causing the flange to almost fully debond from skin. In a second step, a two dimensional plane-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, computed principal stresses exceeded the transverse strength of the material in those areas of the flange where the matrix cracks had developed during the tests. In a third step, delaminations of various lengths were simulated in two locations where delaminations were observed during the tests. The analyses showed that at the loads corresponding to matrix ply crack initiation computed strain energy release rates exceeded the values obtained from a mixed mode failure criterion in one location, Hence. Unstable delamination propagation is likely to occur as observed in the experiments.

Krueger, Ronald↗

Testing and Analysis of Composite Skin/Stringer Debonding under Multi-Axial Loading

A consistent step-wise approach is presented to investigate the damage mechanism in composite bonded skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions. The approach uses experiments to detect the failure mechanism, computational stress analysis to determine the location of first matrix cracking and computational fracture mechanics to investigate the potential for delamination growth. In a first step, tests were performed on specimens, which consisted of a tapered composite flange, representing a stringer or frame, bonded onto a composite skin. Tests were performed under monotonic loading conditions in tension, three-point bending, and combined tension/bending to evaluate the debonding mechanisms between the skin and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load frame was used that was capable of applying both in-plane tension and out-of-plane bending loads simultaneously. Specimen edges were examined on the microscope to document the damage occurrence and to identify typical damage patterns. For all three load cases, observed failure initiated in the flange, near the flange tip, causing the flange to almost fully debond from the skin. In a second step, a two-dimensional plane-strain finite element model was developed to analyze the different test cases using a geometrically nonlinear solution. For all three loading conditions, computed principal stresses exceeded the transverse strength of the material in those areas of the flange where the matrix cracks had developed during the tests. In a third step, delaminations of various lengths were simulated in two locations where delaminations were observed during the tests. The analyses showed that at the loads corresponding to matrix ply crack initiation computed strain energy release rates exceeded the values obtained from a mixed mode failure criterion in one location. Hence, unstable delamination propagation is likely to occur as observed in the experiments.

Krueger, Ronald↗

Tectonic-Climate Interactions in Action Orogenic Belts: Quantification of Dynamic Topography with SRTM data

This project was undertaken to examine the approach to steady state in collisional mountain belts. Although the primary thrust of this grant was to look at larger collisional mountain belts, such as the Himalaya, the Tien Shan, and Southern Alps, we began by looking at smaller structures represented by growing and propagating folds. Like ranges that are evolving toward a topographic steady state, these folds undergo a series of morphologic changes as they are progressively uplifted and eroded. We wanted to document the nature of these changes and to try to discern some of the underlying controls on them. We initially focused on the Wheeler Ridge anticline in southern California. Subsequently, we progressed to looking at the topographic development and the effects of differential uplift and glaciation on the Kyrgyz Range in the northern Tien Shan. This range is unusual inasmuch as it is transformed along its length from a simple uplift with a largely preserved Mesozoic erosion surface arching across it to a highly dissected and heavily glaciated uplift in the region where uplift has been sustained at higher rates over longer intervals. In efforts to understand the distribution of erosion rates at 10(exp 3) - 10(exp 5) year time scales, cosmogenic radionuclide (CRN) concentrations have been gaining increasingly widespread usage (Brown et al., 1995; Riebe et al., 2004; Riebe et al., 2001; Vance et al., 2003). Most studies to date, however, have been conducted in slowly eroding ranges. In rapidly eroding mountains where landslides deliver most of the sediments to the rivers, we hypothesized that CRN concentrations could be highly perturbed by the stochastic processes of landsliding. Therefore, we undertook the development of a numerical model that simulated the effects of both landsliding and grain-by-grain attrition within fluvial catchments. This modeling effort has shown the effects of catchment size and erosion rate on CRN concentrations and allows a prediction of where to sample to obtain the optimal erosion rate estimates using CRN techniques. Finally, we developed computational techniques to operate on DEMs to extract useful information that would enable quantification of climate-erosion interactions. In particular, we worked on rapid techniques to define catchments of any given range of sizes, to extract channel gradients, to combine precipitation information to calculate discharge, and to utilize various stream-power models to determine the erosional energy within any given catchment within a transect. We briefly describe results from Wheeler Ridge, the Kyrgyz Range, the Nepal Himalaya, and our numerical modeling.

Burbank, Douglas W.↗

HFGMC Enhancement of MAC/GMC

Additional information about a mathematical model denoted the high-fidelity generalized method of cells (HFGMC) and implementation of the HFGMC within version 4.0 of the MAC/GMC software has become available. MAC/GMC (Micromechanics Analysis Code With Generalized Method of Cells) was a topic of several prior NASA Tech Briefs articles, version 4.0 having been described in "Comprehensive Micromechanics-Analysis Code - Version 4.0" (LEW-17495-1), NASA Tech Briefs, Vol. 29, No. 9 (September 2005), page 54. MAC/GMC predicts elastic and inelastic thermomechanical responses of composite materials. MAC/GMC utilizes the generalized method of cells (GMC) - a model of micromechanics that predicts macroscopic responses of a composite material as functions of the properties, sizes, shapes, and responses of its constituents (e.g., matrix and fibers). The accuracy of the GMC is limited by neglect of coupling between normal and shear stresses. The HFGMC was developed by combining elements of the GMC and a related model, denoted the higher-order theory for functionally graded materials (HOTFGM), that can account for this coupling. Hence, the HFGMC enables simulation of stress and strain with greater accuracy. Some alterations of the MAC/GMC data structure were necessitated by the greater computational complexity of the HFGMC.

Arnold, Steven M.↗

Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models

A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.

Saether, Erik↗

Experimental Database with Baseline CFD Solutions: 2-D and Axisymmetric Hypersonic Shock-Wave/Turbulent-Boundary-Layer Interactions

A database compilation of hypersonic shock-wave/turbulent boundary layer experiments is provided. The experiments selected for the database are either 2D or axisymmetric, and include both compression corner and impinging type SWTBL interactions. The strength of the interactions range from attached to incipient separation to fully separated flows. The experiments were chosen based on criterion to ensure quality of the datasets, to be relevant to NASA's missions and to be useful for validation and uncertainty assessment of CFD Navier-Stokes predictive methods, both now and in the future. An emphasis on datasets selected was on surface pressures and surface heating throughout the interaction, but include some wall shear stress distributions and flowfield profiles. Included, for selected cases, are example CFD grids and setup information, along with surface pressure and wall heating results from simulations using current NASA real-gas Navier-Stokes codes by which future CFD investigators can compare and evaluate physics modeling improvements and validation and uncertainty assessments of future CFD code developments. The experimental database is presented tabulated in the Appendices describing each experiment. The database is also provided in computer-readable ASCII files located on a companion DVD.

Database↗

On the Comparison of the Long Penetration Mode (LPM) Supersonic Counterflowing Jet to the Supersonic Screech Jet

The authors provide a brief overview of the classic tonal screech noise problem created by underexpanded supersonic jets, briefly describing the fluid dynamic-acoustics feedback mechanism that has been long established as the basis for this well-known aeroacoustics problem. This is followed by a description of the Long Penetration Mode (LPM) supersonic underexpanded counterflowing jet phenomenon which has been demonstrated in several wind tunnel tests and modeled in several computational fluid dynamics (CFD) simulations. The authors provide evidence from test and CFD analysis of LPM that indicates that acoustics feedback and fluid interaction seen in LPM are analogous to the aeroacoustics interactions seen in screech jets. Finally, the authors propose applying certain methodologies to LPM which have been developed and successfully demonstrated in the study of screech jets and mechanically induced excitation in fluid oscillators for decades. The authors conclude that the large body of work done on jet screech, other aeroacoustic phenomena, and fluid oscillators can have direct application to the study and applications of LPM counterflowing supersonic cold flow jets.

Farr, Rebecca A.↗

Computer Simulation and Modeling of CO2 Removal Systems for Exploration 2013-2014

The Atmosphere Revitalization Recovery and Environmental Monitoring (ARREM) project was initiated in September of 2011 as part of the Advanced Exploration Systems (AES) program. Under the ARREM project and the follow-on Life Support Systems (LSS) project, testing of sub-scale and full-scale systems has been combined with multiphysics computer simulations for evaluation and optimization of subsystem approaches. In particular, this paper will describes the testing and 1-D modeling of the combined water desiccant and carbon dioxide sorbent subsystems of the carbon dioxide removal assembly (CDRA). The goal is a full system predictive model of CDRA to guide system optimization and development.

Coker, R.↗

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jet flow control used to suppress boundary layer separation and increase the maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate streamwise direction. It also generates turbulent eddies at the oscillation frequency, which are typically large relative to the scales of boundary layer turbulence, and which augment mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from a nozzle downstream of the oscillator, and an array of sweep jets which suppresses boundary layer separation are performed. Simulation results are compared to data from a dedicated validation experiment of a single oscillator and its sweep jet, and from a wind tunnel test of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets. A critical step in the work is the development of realistic time-dependent sweep jet inflow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the computational fluid dynamics (CFD) solver Overow, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗

A comparison between implicit and hybrid methods for the calculation of steady and unsteady inlet flows

Numerical simulation of steady and unsteady transonic diffuser flows using two different computer codes are discussed and compared with experimental data. The codes solve the Reynolds-averaged, compressible, Navier-Stokes equations using various turbulence models. One of the codes has been applied extensively to diffuser flows and uses the hybrid method of MacCormack. This code is relatively inefficient numerically. The second code, which was developed more recently, is fully implicit and is relatively efficient numerically. Simulations of steady flows using the implicit code are shown to be in good agreement with simulations using the hybrid code. Both simulations are in good agreement with experimental results. Simulations of unsteady flows using the two codes are in good qualitative agreement with each other, although the quantitative agreement is not as good as in the steady flow cases. The implicit code is shown to be eight times faster than the hybrid code for unsteady flow calculations and up to 32 times faster for steady flow calculations. Results of calculations using alternative turbulence models are also discussed.

Coakley, T. J.↗

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗

Computer-Aided Modeling and Analysis of Power Processing Systems (CAMAPPS), phase 1

The large-signal behaviors of a regulator depend largely on the type of power circuit topology and control. Thus, for maximum flexibility, it is best to develop models for each functional block a independent modules. A regulator can then be configured by collecting appropriate pre-defined modules for each functional block. In order to complete the component model generation for a comprehensive spacecraft power system, the following modules were developed: solar array switching unit and control; shunt regulators; and battery discharger. The capability of each module is demonstrated using a simplified Direct Energy Transfer (DET) system. Large-signal behaviors of solar array power systems were analyzed. Stability of the solar array system operating points with a nonlinear load is analyzed. The state-plane analysis illustrates trajectories of the system operating point under various conditions. Stability and transient responses of the system operating near the solar array's maximum power point are also analyzed. The solar array system mode of operation is described using the DET spacecraft power system. The DET system is simulated for various operating conditions. Transfer of the software program CAMAPPS (Computer Aided Modeling and Analysis of Power Processing Systems) to NASA/GSFC (Goddard Space Flight Center) was accomplished.

Kim, S.↗

Subgrid Scale Modeling in Solar Convection Simulations using the ASH Code

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Young, Y.-N.↗

Upset susceptibility study employing circuit analysis and digital simulation

An approach to predict the susceptibility of digital systems to signal disturbances is described. Electrical disturbances on a digital system's input and output lines can be induced by activities and conditions including static electricity, lightning discharge, electromagnetic interference (EMI), and electromagnetic pulsation (EMP). The electrical signal disturbances employed for the susceptibility study were limited to nondestructive levels, i.e., the system does not sustain partial or total physical damage and reset and/or reload brings the system to an operational status. The front-end transition from the electrical disturbances to the equivalent digital signals was accomplished by computer-aided circuit analysis. The super-sceptre (system for circuit evaluation of transient radiation effects) programs was used. Gate models were developed according to manufacturers' performance specifications and parameters resulting from construction processes characteristic of the technology. Digital simulation at the gate and functional level was employed to determine the impact of the abnormal signals on system performance and to study the propagation characteristics of these signals through the system architecture. Example results are included for an Intel 8080 processor configuration.

Carreno, V. A.↗

Upset susceptibility study employing circuit analysis and digital simulation

This paper describes an approach to predicting the susceptibility of digital systems to signal disturbances. Electrical disturbances on a digital system's input and output lines can be induced by activities and conditions including static electricity, lightning discharge, Electromagnetic Interference (EMI) and Electromagnetic Pulsation (EMP). The electrical signal disturbances employed for the susceptibility study were limited to nondestructive levels, i.e., the system does not sustain partial or total physical damage and reset and/or reload will bring the system to an operational status. The front-end transition from the electrical disturbances to the equivalent digital signals was accomplished by computer-aided circuit analysis. The Super-Sceptre (system for circuit evaluation of transient radiation effects) Program was used. Gate models were developed according to manufacturers' performance specifications and parameters resulting from construction processes characteristic of the technology. Digital simulation at the gate and functional level was employed to determine the impact of the abnormal signals on system performance and to study the propagation characteristics of these signals through the system architecture. Example results are included for an Intel 8080 processor configuration.

Carreno, V. A.↗

A Three-Dimensional, Time-Dependent Numerical Modeling of Super-Sonic, Super-Alfvenic MHD Flow

A three-dimensional, time-dependent, MHD disturbance propagating through a steady inhomogeneous MHD flow is analyzed by a numerical model based on the Lax-Wendroff finite difference method. The initial steady MHD flow is assumed to be supersonic and super-Alfvenic. The background steady MHD flow is obtained by dropping the time-dependent terms in the unsteady MHD equations. A transient MHD disturbance corresponding to a MHD shock is generated at the lower radial boundary surface. The resulting development and propagation of 3-D MHD disturbances interacting with the background steady MHD flow are numerically simulated in this first application to the field of interplanetary physics. Detailed descriptions of the mathematical formulation and numerical methodology are presented. Example computations are performed to demonstrate the validity of the formulation and the capabilities of the model.

Han, S. M.↗

Advanced Turbulence Modeling Concepts

The ZCET program developed at NASA Glenn Research Center is to study hydrogen/air injection concepts for aircraft gas turbine engines that meet conventional gas turbine performance levels and provide low levels of harmful NOx emissions. A CFD study for ZCET program has been successfully carried out. It uses the most recently enhanced National combustion code (NCC) to perform CFD simulations for two configurations of hydrogen fuel injectors (GRC- and Sandia-injector). The results can be used to assist experimental studies to provide quick mixing, low emission and high performance fuel injector designs. The work started with the configuration of the single-hole injector. The computational models were taken from the experimental designs. For example, the GRC single-hole injector consists of one air tube (0.78 inches long and 0.265 inches in diameter) and two hydrogen tubes (0.3 inches long and 0.0226 inches in diameter opposed at 180 degree). The hydrogen tubes are located 0.3 inches upstream from the exit of the air element (the inlet location for the combustor). To do the simulation, the single-hole injector is connected to a combustor model (8.16 inches long and 0.5 inches in diameter). The inlet conditions for air and hydrogen elements are defined according to actual experimental designs. Two crossing jets of hydrogen/air are simulated in detail in the injector. The cold flow, reacting flow, flame temperature, combustor pressure and possible flashback phenomena are studied. Two grid resolutions of the numerical model have been adopted. The first computational grid contains 0.52 million elements, the second one contains over 1.3 million elements. The CFD results have shown only about 5% difference between the two grid resolutions. Therefore, the CFD result obtained from the model of 1.3-million grid resolution can be considered as a grid independent numerical solution. Turbulence models built in NCC are consolidated and well tested. They can handle both coarse and fine grids near the wall. They can model the effect of anisotropy of turbulent stresses and the effect of swirling. The chemical reactions of Magnusson model and ILDM method were both used in this study.

Shih, Tsan-Hsing↗