Search NASA⌕ Search

SEARCH · Search NASA

Results for “Problem structuring methods”

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,405 records · Page 78

A perturbation technique for parameter identification in distributed structures

Structures are often characterized by parameters, such as mass and stiffness, that are spatially distributed. Parameter identification of distributed structures is subject to many of the difficulties involved in the modeling problem, and the choice of the model can greatly affect the results of the parameter identification process. Analogously to control spillover in the control of distributed-parameter systems, identification spillover is shown to exist as well and its effect is to degrade the parameter estimates. Moreover, as in modeling by the Rayleigh-Ritz method, it is shown that, for a Rayleigh-Ritz type identification algorithm, an inclusion principle exists in the identification of distributed-parameter systems as well, so that the identified natural frequencies approach the actual natural frequencies monotonically from above.

Meirovitch, L.↗

Use of the quasilinearization algorithm for the simulation of LSS slewing

The use of the Maximum Principle for the large angle slewing of large space structures (LSS) usually results in the so-called two-point boundary-value problem, in which many requirements (e.g., minimum time, small amplitude, and limited control power, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical algorithm. There are many candidate algorithms available for this problem (e.g., quasilinearization, gradient, etc.). Here researchers discuss only the quasilinearization method which has been used for several cases of large angle slewing of LSS. The basic idea of this algorithm is to make a series of successive approximations of the solution from a particular solvable case (linear or nonlinear) to a more general practical case. For the rigid spacecraft slewing problem with no constraints on the controls, the solution procedure can be found in the literature. This procedure needs to be modified if a minimum time for the slewing problem is desired with control limits given. Recently, an indirect method for finding the minimum time was developed to meet all these requirements. For the general mixed (including both rigid and flexible parts) problem, an additional constraint of small vibrational amplitude on the flexible parts is imposed.

Li, Feiyue↗

On the problem of modeling for parameter identification in distributed structures

Structures are often characterized by parameters, such as mass and stiffness, that are spatially distributed. Parameter identification of distributed structures is subject to many of the difficulties involved in the modeling problem, and the choice of the model can greatly affect the results of the parameter identification process. Analogously to control spillover in the control of distributed-parameter systems, identification spillover is shown to exist as well and its effect is to degrade the parameter estimates. Moreover, as in modeling by the Rayleigh-Ritz method, it is shown that, for a Rayleigh-Ritz type identification algorithm, an inclusion principle exists in the identification of distributed-parameter systems as well, so that the identified natural frequencies approach the actual natural frequencies monotonically from above.

Norris, Mark A.↗

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗

Searching for Planets with the Space Interferometry Mission

The Space Interferometry Mission (SIM) will be the first space-based long baseline Michelson interferometer designed for precision astrometry. It will address a wide range of problems in stellar astrophysics and Galactic structure, delivering precision astrometry of stars down to 20 magnitude throughout the entire Galaxy. SIM uses a 10-m Michelson interferometer in Earth-trailing solar orbit to provide 4 microarcsecond (gas) precision astrometry. With a 5-year mission lifetime, SIM will be a powerful tool for discovering planets around nearby stars, through detection of the stellar reflex motion. The astrometric method complements the radial velocity technique which as already yielded many new planets, with an important benefit of directly measuring planetary masses. SIM will have a single-measurement precision of 1 microarcsecond in a frame defined by nearby reference stars, enabling searches for planets with masses as small as a few earth masses around the nearest stars. More massive planets will be detectable to much larger distances. In addition to precision astrometry SIM will also serve an important role as a technology precursor for future astrophysics missions using interferometers. Two technologies demonstrated will be high dynamic-range aperture synthesis imaging at 10-milliarcsec resolution in the optical, and fringe nulling to 10 (exp -4).

Unwin, Stephen↗

Analytical Modal Analysis for Thin-Film Flat Lenses

Due to strong potential applications and more demanding requirements imposed upon thin-film structures for space deployable, there has been increasing research and development activities during recent years in the field of vibration analysis of these types of structures. Moreover, interests in employing these structural components have received renewed emphasis in recent years within NASA and the Air Force. This is due to their inherent lightweight, low packaging and launch volume, and relative simplicity of deployment. Among the potential mission concepts for which these structural elements are included, one can mention solar sails, space solar power generation systems, solar thermal propulsion vehicles, large space telescopes, and inflatable communication antennas. This paper presents analytical procedures to determine vibration and physical characteristics of thin film lenses with circular and elliptical shapes membranes considered in design of a solar concentrator. In general, three methods are used to obtain approximate solutions of Helmholtz boundary value problems. One method requires that solution satisfy the differential equation exactly and the boundary condition approximately. Another method demands a solution that satisfies the boundary conditions exactly and the governing equations approximately. The third method sees a solution that satisfies both the governing equation and boundary conditions approximately. Extensive reviews of vibrations of membrane and plates are provided by Leissa and Mazumdar.

Hamid R. Hamidzadeh↗

Antenna development for the Space Shuttle Orbiter vehicle

A number of difficulties have to be overcome in connection with the design and the mounting of the antennas of the Space Shuttle Orbiter. These difficulties, in combination with a service life goal of 100,000 h, made it necessary to develop antennas with special design features. The particular problems associated with developing antennas for the Shuttle Orbiter are discussed. It is pointed out that this vehicle must function as both a spacecraft and as an aircraft. This feature imposes some rather severe design constraints on the Shuttle antennas. The various structural, operational, and environmental design requirements are considered and design techniques which are used to meet these requirements are examined. A description is presented of the methods employed for testing the various antennas, both with and without the thermal protection system.

Cubley, H. D.↗

A family of solution algorithms for nonlinear structural analysis based on relaxation equations

A family of hierarchical algorithms for nonlinear structural equations are presented. The algorithms are based on the Davidenko-Branin type homotopy and shown to yield consistent hierarchical perturbation equations. The algorithms appear to be particularly suitable to problems involving bifurcation and limit point calculations. An important by-product of the algorithms is that it provides a systematic and economical means for computing the stepsize at each iteration stage when a Newton-like method is employed to solve the systems of equations. Some sample problems are provided to illustrate the characteristics of the algorithms.

Park, K. C.↗

Early surface habitation elements for planetary exploration missions

Rigid pressure vessel concepts appropriate for initial, crew-supporting exploration activities on planetary surfaces are discussed and compared. The problem of early exploration crew support is first bounded and outlined. The studied trade-space spans five crew sizes from 4 to 12, and three diameters of module structure from 4.4 to 10 m. Other considerations of gravity, orientation, structure and topology are also used as parameters. The paper illustrates application of a technique which uses interpolation across a field of point-designs to yield trade results. Methods to reduce the option set from 1480 to 30 are presented. Reproducible evaluation criteria are explained and applied, to enable selecting the most promising candidates. A final concept is synthesized for further study use, and presented via configuration drawings and weight statements. It is a 7.6 m diameter module, with two floors running lengthwise and a cross-sectioning pressure bulkhead. Versions for four, six and perhaps eight people can be landed fully integrated on the moon.

Sherwood, Brent↗

Modeling dynamically coupled fluid-duct systems with finite line elements

Structural analysis of piping systems, especially dynamic analysis, typically considers the duct structure and the contained fluid column separately. Coupling of these two systems, however, forms a new dynamic system with characteristics not necessarily described by the superposition of the two component system's characteristics. Methods for modeling the two coupled components simultaneously using finite line elements are presented. Techniques for general duct intersections, area or direction changes, long radius bends, hydraulic losses, and hydraulic impedances are discussed. An example problem and results involving time transients are presented. Additionally, a program to enhance post-processing of line element models is discussed.

Saxon, J. B.↗

Residual life and strength estimates of aircraft structural components with MSD/MED

Economic and safe operation of the flight vehicles flying beyond their initial design life calls for an in-depth structural integrity evaluation of all components with potential for catastrophic damages. Fuselage panels with cracked skin and/or stiffening elements is one such example. A three level analytical approach is developed to analyze the pressurized fuselage stiffened shell panels with damaged skin or stiffening elements. A global finite element analysis is first carried out to obtain the load flow pattern through the damaged panel. As an intermediate step, the damaged zone is treated as a spatially three-dimensional structure modeled by plate and shell finite elements, with all the neighboring elements that can alter the stress state at the crack tip. This is followed by the Schwartz-Neumann alternating method for local analysis to obtain the relevant crack tip parameters that govern the onset of fracture and the crack growth. The methodology developed is generic in nature and aims at handling a large fraction of problem areas identified by the Industry Committee on Wide-Spread Fatigue Damage.

Singh, Ripudaman↗

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗

Lunar Habitat Optimization Using Genetic Algorithms

Long-duration surface missions to the Moon and Mars will require bases to accommodate habitats for the astronauts. Transporting the materials and equipment required to build the necessary habitats is costly and difficult. The materials chosen for the habitat walls play a direct role in protection against each of the mentioned hazards. Choosing the best materials, their configuration, and the amount required is extremely difficult due to the immense size of the design region. Clearly, an optimization method is warranted for habitat wall design. Standard optimization techniques are not suitable for problems with such large search spaces; therefore, a habitat wall design tool utilizing genetic algorithms (GAs) has been developed. GAs use a "survival of the fittest" philosophy where the most fit individuals are more likely to survive and reproduce. This habitat design optimization tool is a multiobjective formulation of up-mass, heat loss, structural analysis, meteoroid impact protection, and radiation protection. This Technical Publication presents the research and development of this tool as well as a technique for finding the optimal GA search parameters.

SanScoucie, M. P.↗

Experimental Flight Validation of the Prandtl 1933 Bell Spanload

This report describes the validation of the 1933 Prandtl bell spanload. This spanload is the minimum induced drag of a wing for a given structural weight with properties that eliminate adverse yaw. Aircraft using the Prandtl bell spanload were flown and investigated. The results of this research show that many previously held assumptions should be rethought, and the creation of aircraft using the Prandtl bell spanload will require considerable new techniques. Part of this work centered on the use of inverse methods. The usual first-step approach to a computational fluids problem is to create a geometry of the aircraft. Once this geometry exists, the computational fluids solution has been solved; however, the problem of creating the geometry still exists. A very small segment of the computational fluids world has concentrated on inverse solutions. This design approach begins with an end result - the computational fluids solution; from this end result, the geometry is sought. A more generic inverse tool was desired that would allow for the design of wings - specifically Prandtl bell spanload wings. Such a tool has been developed and allows for variable taper, aspect ratio, sweep, airfoils, and design-lift coefficients; output from the tool results in twist distribution of wings.

Albion H. Bowers↗

Thin-film temperature sensors for gas turbine engines Problems and prospects

The erosion and corrosion of thermocouples used to measure the temperature in turbine engines are studied. Structural and metallurgical interactions and instabilities at thermocouple interfaces are analyzed. Consideration is given to the adhesion, dielectric quality, surface topography, and hardness of the thermal oxides; it is observed that the structural and thermoelectric stability of thin-film thermocouple elements depends on adhesion, surface topography, and dielectric strength. The electrical conductivity and impurity content of the oxide scale are evaluated. Methods for improving the adhesion of thermocouples on the alumina surfaces are described. Compositional inhomogeneities in the sensors and contamination of the thermocouple elements are examined. The fabrication of the thermocouples is discussed. It is noted that Al2O3 and Si3N4 are useful for developing stable thermocouple elements on the surface of the blades and vanes.

Budhani, R. C.↗

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗

Derivation of Effective Properties Based on Porous Scale Simulations Using Filtering Techniques

This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.

Alexsander Zibitsker↗

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗