Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient method”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Microwave Imaging on Metal Objects

This final report for the project discusses the attempts to model, using different methods, microwave image reconstruction. Maximum Entropy Method was not successful. Attempts to use Singular Value Decomposition (SVD) got some good results after initial failure. SVD is based upon a theory of linear algebra, to the effect that any M X N Matrix A whose number of rows M is greater than or equal to its number of columns, N can be written as the product of an M X N column-orthogonal matrix U, an N X N diagonal Matrix, W, with m positive or zero elements (the singular values) and the transposition of an N X N orthogonal matrix V. In microwave imaging, the scattered fields can be expressed by the induced current distribution. The SVD method required more contiguous computer memory than was available. Work was also done on the Conjugate Gradient Method (CGM), which didn't work well when tried earlier. It was found that separation of the imaginary part and the real part during calculation may work. This work was considered incomplete as of the end of the grant period.

Tolliver, C. L.↗
NASA NTRS2026

AC and DC Fault Management for Megawatt Electrified Aircraft Electrical Powertrains Task 3: Lifetime and Reliability of Electrical Insulators

This research project was a collaborative investigation between researchers at the RTX Technology Research Center (RTRC) and the University of Texas at Austin and made a significant contribution to enabling electric aircraft. The transport of electric power between the points of generation and use requires power cables. These cables must be smaller, lighter and provide a more predictable life than power cables used in stationary applications. Consequently, this investigation provided heretofore unavailable information supporting the safety and reliability of smaller lighter power cables for electrified aircraft. In addition, the research identified key additional engineering data needed to support quantitative reliability assessments. Important advances included: • Demonstrated that at least one manufacturer can make a novel, smaller, lighter power cable that is free from serious defects. • Developed and published an appropriate analytical construct to describe the life of this novel cable. This is a necessary step for use in aviation where the understanding of remaining life is critical. • Demonstrated thermal-mechanical aging that suggested 1000+ flights before the thermal-mechanical processes produced defects large enough that the defect growth was accelerated electrically. • Showed that electrical aging took place at two rates. The first possibly lasting weeks to months and the second possibly days to weeks. If robust, this provides a good diagnostic for cable replacement. • Demonstrated that the traditional electrical testing of cable materials using manufactured voids can be misleading due to the size of the voids. Emerging laser drilling technology permitted demonstration that the physics of failure in realistically small voids is different from that in the unrealistically large voids used in earlier research, which is very important for high-quality, high-performance, small aircraft cables. Although this project represents a significant contribution to the specifics of cable aging in the aircraft environment, important additional research remains to be completed, including: • Non-uniform thermal cycling by applying the heat from the center conductor to maximize thermal stress next to the core area where the electric gradient is the strongest. This builds on the uniform thermal cycling that has been completed. • The augmentation of the thermal-mechanical failure rate by electrical processes. Better understanding of these time constants strongly affects the ability to predict life. • Termination design: Terminations provide not only electrical reflection potential, but a location for a series arc fault and an area where ozone can diffuse into the center conductor and negatively affect cable insulation. • The abrasion and ozone resistance of the cable jacket. • Pressure cycling as an accelerant of thermal, mechanical, and/or electrical aging. • Possible methods for online PD detection and offline PD localization

model↗
NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗
NASA NTRS1987

Laser schlieren crystal monitor

A system and method for monitoring the state of a crystal which is suspended in a solution is described which includes providing a light source for emitting a beam of light along an optical axis. A collimating lens is arranged along the optical axis for collimating the emitted beam to provide a first collimated light beam consisting of parallel light rays. By passing the first collimated light beam through a transparent container, a number of the parallel light rays are deflected off the surfaces of said crystal being monitored according to the refractive index gradient to provide a deflected beam of deflected light rays. A focusing lens is arranged along optical axis for focusing the deflected rays towards a desired focal point. A knife edge is arranged in a predetermined orientation at the focal point; and a screen is provided. A portion of the deflected beam is blocked with the knife edge to project only a portion of the deflected beam. A band is created at one edge of the image of the crystal which indicates the state of change of the surface of the crystal being monitored.

Owen, Robert B.↗
NASA NTRS1997

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗
NASA NTRS1990

Research on diabatic initialization

The objective of this research project is to contribute to improvement in the synoptic analyses in the tropics for numerical weather prediction and climate research. In addition to a prediction model, four dimensional data assimilation systems have two principal components. One is objective analysis and the other is initialization. Various methods of objective analysis are designed primarily to analyze the mass and rotational wind fields. Methods of initialization are developed to obtain the irrotational wind and its associated vertical velocity field which are balanced with the mass field and free from meteorological noise. There are essentially three approaches to the problem of initialization: quasi-geostrophic theory, bounded derivative method and nonlinear normal mode method. In the midlatitudes, these approaches generally produce satisfactory results even without diabatic effects for large-scale motions. In the tropics, the situation is quite different from that in the midlatitudes. Because of a small magnitude of the Coriolis parameter and a weak horizontal temperature gradient in the tropics, any method of initialization must incorporate diabatic effects. In fact, it can be said that understanding the problem of diabatic initialization is the key to improving the analysis and weather forecasting in the tropics.

Kasahara, Akira↗
NASA NTRS1996

Cascade Error Projection: A Learning Algorithm for Hardware Implementation

In this paper, we workout a detailed mathematical analysis for a new learning algorithm termed Cascade Error Projection (CEP) and a general learning frame work. This frame work can be used to obtain the cascade correlation learning algorithm by choosing a particular set of parameters. Furthermore, CEP learning algorithm is operated only on one layer, whereas the other set of weights can be calculated deterministically. In association with the dynamical stepsize change concept to convert the weight update from infinite space into a finite space, the relation between the current stepsize and the previous energy level is also given and the estimation procedure for optimal stepsize is used for validation of our proposed technique. The weight values of zero are used for starting the learning for every layer, and a single hidden unit is applied instead of using a pool of candidate hidden units similar to cascade correlation scheme. Therefore, simplicity in hardware implementation is also obtained. Furthermore, this analysis allows us to select from other methods (such as the conjugate gradient descent or the Newton's second order) one of which will be a good candidate for the learning technique. The choice of learning technique depends on the constraints of the problem (e.g., speed, performance, and hardware implementation); one technique may be more suitable than others. Moreover, for a discrete weight space, the theoretical analysis presents the capability of learning with limited weight quantization. Finally, 5- to 8-bit parity and chaotic time series prediction problems are investigated; the simulation results demonstrate that 4-bit or more weight quantization is sufficient for learning neural network using CEP. In addition, it is demonstrated that this technique is able to compensate for less bit weight resolution by incorporating additional hidden units. However, generation result may suffer somewhat with lower bit weight quantization.

Duong, Tuan A.↗
NASA NTRS2019

Optimal Control within the Context of Multidisciplinary Design, Analysis, and Optimization

Multidisciplinary design, analysis and optimization involves modeling the interactions of complex systems across a variety of disciplines. The optimization of such systems can be a computationally expensive exercise with multiple levels of nested nonlinear solvers running under an optimizer.The application of optimal control in project development often involves performing trajectory optimization for fixed vehicle designs or parametric sweeps across some key vehicle properties.This information is then relayed to the subsystem design teams who update their designs and relay some bulk characteristics back to the trajectory optimization procedure.This iteration is then repeated until the design closes.However, with increasing interest in more tightly coupled systems, such as electric and hybrid-electric aircraft propulsion and boundary layer ingestion, this process is prone to ignore subtle coupling between vehicle subsystem designs and vehicle operation on a given mission.Integrating trajectory optimization into a tightly coupled multidisciplinary design procedure can be computationally prohibitive, depending on the complexity of the subsystem analyses and the optimal control technique applied.To address these issues a new optimal control software tool, Dymos, has been developed.Dymos is built upon NASA's OpenMDAO software and can leverage its capabilities to efficiently compute gradients for the optimization and optimize complex models in parallel on distributed memory systems.This report provides some explanation into the numerical methods employed in Dymos and provides several use cases that demonstrate its performance on traditional optimal control problems and improvements ino techniques have been used extensively in recent decades to solve a variety of optimal control problems, typically in the form of aerospace vehicle trajectory optimization.

pseudospectral↗
NASA NTRS2013

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗
NASA NTRS1999

Project Description and Publications List for UAH CMMR

This research combines a state of the art X-ray Transmission Microscope, XTM, with a specially designed x-ray transparent horizontal Bridgman furnace to image (with resolutions up to 3 micrometers) the solidification of metal alloys in real-time. The objective is to obtain real-time dynamic data to provide direct measure of the solute profile in the liquid, phase coalescence and growth in the liquid, and the detailed interface morphology (e,g., dendrites and cells) during solidification. We are also enhancing the XTM data with precise solid-liquid interfacial temperature and the thermal gradient measurement techniques, and working on the application of this technology to the study of the fundamentals of solidification in microgravity. Over the last several years we have successfully imaged in real-time: interfacial-morphologies, phase growth, coalescence, incorporation of phases into the growing interface, and the solute boundary layer in the liquid at the solid-liquid interface. We have also measured true local growth rates and can evaluate segregation structures in the solid. Interfacial undercoolings are being measured either with a special Seebeck furnace or with micro-thermocouple arrays we are developing. These later techniques are presently being incorporated with the XTM furnace. This last year emphasized the investigation of the solute layer in the melt during solidification. Methods were developed to quantify the solute concentrations using x-ray absorption and to compare to predictions from simulations. In addition, work is being completed on a brass-board portable XTM that incorporates a vertical Bridgman furnace.

Kaukler, William F.↗
NASA NTRS2001

International Space Station Radiation Shielding Model Development

The projected radiation levels within the International Space Station (ISS) have been criticized by the Aerospace Safety Advisory Panel in their report to the NASA Administrator. Methods for optimal reconfiguration and augmentation of the ISS shielding are now being developed. The initial steps are to develop reconfigurable and realistic radiation shield models of the ISS modules, develop computational procedures for the highly anisotropic radiation environment, and implement parametric and organizational optimization procedures. The targets of the redesign process are the crew quarters where the astronauts sleep and determining the effects of ISS shadow shielding of an astronaut in a spacesuit. The ISS model as developed will be reconfigurable to follow the ISS. Swapping internal equipment rack assemblies via location mapping tables will be one option for shield optimization. Lightweight shield augmentation materials will be optimally fit to crew quarter areas using parametric optimization procedures to minimize the augmentation shield mass. The optimization process is being integrated into the Intelligence Synthesis Environment s (ISE s) immersive simulation facility at the Langley Research Center and will rely on High Performance Computing and Communication (HPCC) for rapid evaluation of shield parameter gradients.

Qualls, G. D.↗
NASA NTRS1971

Angular momentum desaturation for Skylab using gravity gradient torques

An angular momentum desaturation method for momentum exchange devices of orbiting spacecraft is described. The specific application of the method is to the Skylab which contains three double-gimbaled control moment gyros for precise attitude control and maneuvering. It is assumed that the attitude reference is inertially fixed and that two of the vehicle principal moments of inertia are much larger than the third. Gravity gradient torques and resultant angular momentum accumulation are developed for small deviations from the reference. The assumed moment-of-inertia distribution allows desaturation about all axes with only two attitude angles each for the two axes with large moments of inertia. The necessary desaturation maneuvers can be decoupled for a special set of orbital coordinates. All maneuvers are made during the night portion of the orbit, and the percentage utilized for desaturation is selectable. Expressions for the attitude angle commands are developed assuming infinite vehicle rates. The effect of finite rates introduces an efficiency into the desaturation. Expressions for this efficiency are developed and means for compensation are treated. Arbitrary misalignments between geometric vehicle axes and principal moment-of-inertia axes are permissible. An angle bias about the sun line minimizes the angular momentum accumulation about the sun line projection into the orbital plane. Adaptive desaturation maneuver limiting consistent with the available maneuver momentum is included.

Kennel, H. F.↗
NASA NTRS1995

The reliability of X-ray constraints of intrinsic cluster shapes

Using the simulation of Katz & White (1993) we have tested the viability of X-ray analysis for constraining the intrinsic shapes of clusters of galaxies considering the effects of both substructure and steep temperature gradients. We restrict our analysis to the aggregate shapes of clusters on scales of r approximately 1-2 Mpc in order to reduce our sensitivity to subclustering in the core. For low redshifts (z less than or approximately = 0.25) the X-ray method accurately measures the true ellipticity of the three-dimensional cluster dark matter provided the inclination of the cluster is known to within approximately 30 deg; assuming the gas is isothermal adds only small errors to the derived shapes. At higher redshifts the X-ray method yields unreliable results since the gas does not trace the cluster gravitational potential. We proffer some necessary conditions for the reliability of X-ray methods characterized by both the amount of substructure in the X-ray surface brightness images and the shapes of the isophotes. We conclude that measurements of the aggregate shapes of clusters on scales r approximately 1-2 Mpc are insensitive to core substructure representing scales of a few hundred kpc. Therefore our results suggest that the X-ray measurements of aggregate cluster shapes by Fabricant, Rybicki, & Gorenstein (1984) and Buote & Canizares (1992) are valid provided that they do not suffer from serious projection effects. A substantial number of Abell clusters observed with the ROSAT PSPC will be amenable to X-ray shape analysis.

Buote, David A.↗
NASA NTRS2002

Large Coded Aperture Mask for Spaceflight Hard X-ray Images

The 2.6 square meter coded aperture mask is a vital part of the Burst Alert Telescope on the Swift mission. A random, but known pattern of more than 50,000 lead tiles, each 5 mm square, was bonded to a large honeycomb panel which projects a shadow on the detector array during a gamma ray burst. A two-year development process was necessary to explore ideas, apply techniques, and finalize procedures to meet the strict requirements for the coded aperture mask. Challenges included finding a honeycomb substrate with minimal gamma ray attenuation, selecting an adhesive with adequate bond strength to hold the tiles in place but soft enough to allow the tiles to expand and contract without distorting the panel under large temperature gradients, and eliminating excess adhesive from all untiled areas. The largest challenge was to find an efficient way to bond the > 50,000 lead tiles to the panel with positional tolerances measured in microns. In order to generate the desired bondline, adhesive was applied and allowed to cure to each tile. The pre-cured tiles were located in a tool to maintain positional accuracy, wet adhesive was applied to the panel, and it was lowered to the tile surface with synchronized actuators. Using this procedure, the entire tile pattern was transferred to the large honeycomb panel in a single bond. The pressure for the bond was achieved by enclosing the entire system in a vacuum bag. Thermal vacuum and acoustic tests validated this approach. This paper discusses the methods, materials, and techniques used to fabricate this very large and unique coded aperture mask for the Swift mission.

Vigneau, Danielle N.↗
NASA NTRS1998

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.↗
NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗
NASA NTRSDate not supplied

Impact of Thermoplastic Composites: Testing and Modeling

Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.

Composite Materials↗
NASA NTRS1999

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗