Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Adaptation of G-TAG Software for Validating Touch-and-Go Comet Surface Sampling Design Methodology

The G-TAG software tool was developed under the R&TD on Integrated Autonomous Guidance, Navigation, and Control for Comet Sample Return, and represents a novel, multi-body dynamics simulation software tool for studying TAG sampling. The G-TAG multi-body simulation tool provides a simulation environment in which a Touch-and-Go (TAG) sampling event can be extensively tested. TAG sampling requires the spacecraft to descend to the surface, contact the surface with a sampling collection device, and then to ascend to a safe altitude. The TAG event lasts only a few seconds but is mission-critical with potentially high risk. Consequently, there is a need for the TAG event to be well characterized and studied by simulation and analysis in order for the proposal teams to converge on a reliable spacecraft design. This adaptation of the G-TAG tool was developed to support the Comet Odyssey proposal effort, and is specifically focused to address comet sample return missions. In this application, the spacecraft descends to and samples from the surface of a comet. Performance of the spacecraft during TAG is assessed based on survivability and sample collection performance. For the adaptation of the G-TAG simulation tool to comet scenarios, models are developed that accurately describe the properties of the spacecraft, approach trajectories, and descent velocities, as well as the models of the external forces and torques acting on the spacecraft. The adapted models of the spacecraft, descent profiles, and external sampling forces/torques were more sophisticated and customized for comets than those available in the basic G-TAG simulation tool. Scenarios implemented include the study of variations in requirements, spacecraft design (size, locations, etc. of the spacecraft components), and the environment (surface properties, slope, disturbances, etc.). The simulations, along with their visual representations using G-View, contributed to the Comet Odyssey New Frontiers proposal effort by indicating problems and/or benefits of different approaches and designs.

Mandic, Milan↗
NASA NTRS2020

AERoBOND Project Summary

Under NASA’s Convergent Aeronautics Solutions (CAS) project, the Adhesive-Free Bonding of Complex Composites (AERoBOND) project investigated off-stoichiometric epoxy polymers for fast, reliable assembly of epoxy matrix composite structures. The project goal was to demonstrate feasibility of the AERoBOND joining method by demonstrating mechanical properties greater than 80% of conventional co-cured materials while reducing structure weight by 1%. The project consisted of three convergent research areas: material and process development, systems analysis, and material and process modeling. Material and process development was the largest component of AERoBOND with approximately 6 FTE and 1WYE of support to formulate and characterize new resins, prepare carbon fiber prepregs, fabricate laminates, measure mechanical properties, analyze failure results, and select material and process improvements. The systems analysis activity estimated the potential reduction in part count and aircraft weight by comparing models of composite wing boxes with no fasteners (co-cured structure), fasteners in major joints (co-cured stringers), and fasteners in all joints. The materials and process modeling activity included a molecular model of the AERoBOND materials system to predict mechanical properties of resins with offset stoichiometry and a process model to predict the effect of resin formulation and processing conditions on the extent of mixing and degree of cure in a finished joint. As the number of airline passenger trips doubles in the next 20 years (IATA/Tourism Economics Air Passenger Forecasts, April 2019), the increased demand for new commercial aircraft is now the single greatest technical challenge to the airframe manufacturing industry. To meet efficiency requirements, new aircraft must be fabricated primarily from high performance structural composites, but manufacturing processes are inherently slow with the largest bottleneck attributed to assembly and installation of fasteners (NASA/TM–2019-220428). Manufactures of commercial transport aircraft are compelled to install more than 100,000 redundant fasteners into bonded joints to prevent failures due to unpredictable weak bonds. In structural adhesive bonds, the interface between adherend and adhesive is nearly two-dimensional making it susceptible to minute quantities of contamination, which can cause weak bonds. Currently, bond strength assessment is only possible through destructive testing (i.e., breaking the joint). For these reasons, regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded, primary-structures to alleviate concerns with bond performance. The AERoBOND process enables reflow of matrix resin during assembly to eliminate the material discontinuity at the interface, thereby eliminating the dependence of mechanical performance on interfacial adhesion. The AERoBOND joint is equivalent to the interlaminar region obtained during a co-cure process, so joint performance depends on the cohesive properties of the matrix resin. Conventional co-cured structures, although too costly and complex for large-scale manufacturing, are trusted by manufacturers and regulators, and are certified for flight with few or no redundant fasteners.Systems analysis performed on a composite wing model at the scale of a single-aisle commercial transport aircraft indicated that >20,000 redundant fasteners per wing could be eliminated by implementing the AERoBOND joining method. A total weight reduction of 15% was predicted in a wing box by eliminating fasteners and thinning components that must no longer support localized fastener loads and accommodate fastener dimensions. Interlaminar shear fracture toughness measured by the end-notched flexure test was greater than 1 kJ/m2 (nearly 140% of the co-cured benchmark property), which is greatly in excess of the project goals for mechanical properties. Testing was planned to measure interlaminar tensile fracture toughness as well as interlaminar tensile and shear strengths using the same AERoBOND configuration, but was delayed due to closure of LaRC facilities during the COVID-19 pandemic. The AERoBOND process model is partially validated and available for experimental use. It allows the user to input AERoBOND process parameters such as material composition, laminate configuration, and cure cycle to predict the final cure state of the AERoBOND joint. A preliminary, multi-scale material model was developed to predict AERoBOND joint mechanical properties (stiffness and strength) based on the cure state of the joint provided by the process model. The timing for transition of this technology within NASA is excellent as NASA initiates new enduring projects to address composites manufacturing rate challenges. AERoBOND technology is well suited to AAVP/AATT objectives for rapid manufacturing of a composite wing. A minimal effort (1 FTE/$15k procurement/0 WYE) is proposed in FY21 to continue a minor mechanical testing effort and maintain a SAA with ASX composites to develop commercial quality prepreg material. An RFI with the composites industry is suggested to quantify the technology gap between the current TRL and the TRL needed for transition to industry. A moderate effort [3-4 FTE/$150k/1 WYE (~$115k)] is proposed in FY22 for the “high rate composites manufacturing” project currently in planning. The partnership with ASX Composites will be expanded to produce material for sub-element/element-scale “panel-off” activities. Industry partnerships with airframe manufacturers is an expected component to explore damage tolerance and environmental stability. Further development of multi-scale modeling tools (process model, meso-scale model, and molecular model) is planned to enhance and deliver tools for rapid manufacturing infusion.

Frank Louis Palmieri↗
NASA NTRS2018

Overview of the Spacecraft Design for the Psyche Mission Concept

In January 2017, Psyche and a second mission concept were selected by NASA for flight as part of the 14th Discovery mission competition. Assigned for an initial launch date in 2023, the Psyche team was given direction shortly after selection to research the possibility for earlier opportunities. Ultimately, the team was able to identify a launch opportunity in 2022 with a reduced flight time to its destination. This was accomplished in large part to crosscutting trades centered on the electrical power subsystem. These trades were facilitated through the Psyche mission's planned use of Solar Electric Propulsion (SEP), which enables substantial flexibility with respect to trajectory design. In combination with low-thrust trajectory analysis tools, the team was able to robustly converge to solutions with a higher fidelity and accuracy of results. These trades also took advantage of the 1300 series product line produced by Space Systems Loral (SSL), which enabled power growth while maintaining strong system-level heritage through its modular design that has been utilized on a large number of geostationary (GEO) communications satellites. This paper presents an overview of the Psyche mission concept, and the unique architecture that enables the use of commercially developed electric propulsion and space power systems from Space Systems Loral to provide flexibility in mission design. This paper then discusses the trades that allowed the Psyche team to meet a 2022 launch date.

Prikl, Zachary↗
NASA NTRS2013

Statistical Analysis of CFD Solutions From the Fifth AIAA Drag Prediction Workshop

A graphical framework is used for statistical analysis of the results from an extensive N-version test of a collection of Reynolds-averaged Navier-Stokes computational fluid dynamics codes. The solutions were obtained by code developers and users from North America, Europe, Asia, and South America using a common grid sequence and multiple turbulence models for the June 2012 fifth Drag Prediction Workshop sponsored by the AIAA Applied Aerodynamics Technical Committee. The aerodynamic configuration for this workshop was the Common Research Model subsonic transport wing-body previously used for the 4th Drag Prediction Workshop. This work continues the statistical analysis begun in the earlier workshops and compares the results from the grid convergence study of the most recent workshop with previous workshops.

Statistical analysis↗
NASA NTRS2010

Accuracy Analysis for Finite-Volume Discretization Schemes on Irregular Grids

A new computational analysis tool, downscaling test, is introduced and applied for studying the convergence rates of truncation and discretization errors of nite-volume discretization schemes on general irregular (e.g., unstructured) grids. The study shows that the design-order convergence of discretization errors can be achieved even when truncation errors exhibit a lower-order convergence or, in some cases, do not converge at all. The downscaling test is a general, efficient, accurate, and practical tool, enabling straightforward extension of verification and validation to general unstructured grid formulations. It also allows separate analysis of the interior, boundaries, and singularities that could be useful even in structured-grid settings. There are several new findings arising from the use of the downscaling test analysis. It is shown that the discretization accuracy of a common node-centered nite-volume scheme, known to be second-order accurate for inviscid equations on triangular grids, degenerates to first order for mixed grids. Alternative node-centered schemes are presented and demonstrated to provide second and third order accuracies on general mixed grids. The local accuracy deterioration at intersections of tangency and in flow/outflow boundaries is demonstrated using the DS tests tailored to examining the local behavior of the boundary conditions. The discretization-error order reduction within inviscid stagnation regions is demonstrated. The accuracy deterioration is local, affecting mainly the velocity components, but applies to any order scheme.

Diskin, Boris↗
NASA NTRS2002

Higher Order, Hybrid BEM/FEM Methods Applied to Antenna Modeling

In this presentation, the authors address topics relevant to higher order modeling using hybrid BEM/FEM formulations. The first of these is the limitation on convergence rates imposed by geometric modeling errors in the analysis of scattering by a dielectric sphere. The second topic is the application of an Incomplete LU Threshold (ILUT) preconditioner to solve the linear system resulting from the BEM/FEM formulation. The final tOpic is the application of the higher order BEM/FEM formulation to antenna modeling problems. The authors have previously presented work on the benefits of higher order modeling. To achieve these benefits, special attention is required in the integration of singular and near-singular terms arising in the surface integral equation. Several methods for handling these terms have been presented. It is also well known that achieving ~he high rates of convergence afforded by higher order bases may als'o require the employment of higher order geometry models. A number of publications have described the use of quadratic elements to model curved surfaces. The authors have shown in an EFIE formulation, applied to scattering by a PEC .sphere, that quadratic order elements may be insufficient to prevent the domination of modeling errors. In fact, on a PEC sphere with radius r = 0.58 Lambda(sub 0), a quartic order geometry representation was required to obtain a convergence benefi.t from quadratic bases when compared to the convergence rate achieved with linear bases. Initial trials indicate that, for a dielectric sphere of the same radius, - requirements on the geometry model are not as severe as for the PEC sphere. The authors will present convergence results for higher order bases as a function of the geometry model order in the hybrid BEM/FEM formulation applied to dielectric spheres. It is well known that the system matrix resulting from the hybrid BEM/FEM formulation is ill -conditioned. For many real applications, a good preconditioner is required to obtain usable convergence from an iterative solver. The authors have examined the use of an Incomplete LU Threshold (ILUT) preconditioner . to solver linear systems stemming from higher order BEM/FEM formulations in 2D scattering problems. Although the resulting preconditioner provided aD excellent approximation to the system inverse, its size in terms of non-zero entries represented only a modest improvement when compared with the fill-in associated with a sparse direct solver. Furthermore, the fill-in of the preconditioner could not be substantially reduced without the occurrence of instabilities. In addition to the results for these 2D problems, the authors will present iterative solution data from the application of the ILUT preconditioner to 3D problems.

Fink, P. W.↗
NASA NTRS1971

An analysis and comparison of several trajectory optimization methods

The sensitivities of the convergence characteristics of the methods to initially assumed parameters and trial solution, convergence times, computer logic, and storage requirements are discussed. Numerical comparison of the convergence characteristics is made by considering a minimum time, low thrust, Earth-Mars transfer trajectory. A modified quasi-linearization method reduces convergence time by approximately 70% when compared with the generalized Newton-Raphson method and allows the terminal boundary to be specified by a general function of the problem variables. A uniquely specified and easily determined, time dependent weighting matrix for the gradient techniques accelerates the shaping of the optimal control program and improves the convergence characteristics during the terminal iterations. Convergence envelopes, indicating how sensitive the convergence characteristics are to initially assumed parameters, are plotted for the perturbation and quasi-linearization methods. Several iteration schemes are proposed which increase the size of the convergence envelopes and decrease the sensitivity of the method to initially assumed parameters.

Lewallen, J. M.↗
NASA NTRS1988

Analysis of the spectral vanishing viscosity method for periodic conservation laws

The convergence of the spectral vanishing method for both the spectral and pseudospectral discretizations of the inviscid Burgers' equation is analyzed. It is proven that this kind of vanishing viscosity is responsible for a spectral decay of those Fourier coefficients located toward the end of the computed spectrum; consequently, the discretization error is shown to be spectrally small independent of whether the underlying solution is smooth or not. This in turn implies that the numerical solution remains uniformly bounded and convergence follows by compensated compactness arguments.

Maday, Yvon↗
NASA NTRS2023

Designing an Observing System to Study the Surface Biology and Geology (SBG) of the Earth in the 2020s

Observations of planet Earth from space are a critical resource for science and society. Satellite measurements represent very large investments and United States (US) agencies organize their effort to maximize the return on that investment. The US National Research Council conducts a survey of Earth science and applications to prioritize observations for the coming decade. The most recent survey prioritized a visible to shortwave infrared imaging spectrometer and a multispectral thermal infrared imager to meet a range of needs for studying Surface Biology and Geology (SBG). SBG will be the premier integrated observatory for observing the emerging impacts of climate change by characterizing the diversity of plant life and resolving chemical and physiological signatures. It will address wildfire risk, behavior, and recovery as well as responses to hazards such as oil spills, toxic minerals in minelands, harmful algal blooms, landslides, and other geological hazards. The SBG team analyzed needed instrument characteristics (spatial, temporal, and spectral resolutions, measurement uncertainty) and assessed the cost, mass, power, volume, and risk of different architectures. We present an overview of the Research and Applications trade-study analysis of algorithms, calibration and validation needs, and societal applications with specifics of substudies detailed in other articles in this special collection. We provide a value framework to converge from hundreds down to three candidate architectures recommended for development. The analysis identified valuable opportunities for international collaboration to increase the revisit frequency, adding value for all partners, leading to a clear measurement strategy for an observing system architecture.

Observing System↗
NASA NTRS1988

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗
NASA NTRS2005

Convergence Acceleration and Documentation of CFD Codes for Turbomachinery Applications

The development and analysis of turbomachinery components for industrial and aerospace applications has been greatly enhanced in recent years through the advent of computational fluid dynamics (CFD) codes and techniques. Although the use of this technology has greatly reduced the time required to perform analysis and design, there still remains much room for improvement in the process. In particular, there is a steep learning curve associated with most turbomachinery CFD codes, and the computation times need to be reduced in order to facilitate their integration into standard work processes. Two turbomachinery codes have recently been developed by Dr. Daniel Dorney (MSFC) and Dr. Douglas Sondak (Boston University). These codes are entitled Aardvark (for 2-D and quasi 3-D simulations) and Phantom (for 3-D simulations). The codes utilize the General Equation Set (GES), structured grid methodology, and overset O- and H-grids. The codes have been used with success by Drs. Dorney and Sondak, as well as others within the turbomachinery community, to analyze engine components and other geometries. One of the primary objectives of this study was to establish a set of parametric input values which will enhance convergence rates for steady state simulations, as well as reduce the runtime required for unsteady cases. The goal is to reduce the turnaround time for CFD simulations, thus permitting more design parametrics to be run within a given time period. In addition, other code enhancements to reduce runtimes were investigated and implemented. The other primary goal of the study was to develop enhanced users manuals for Aardvark and Phantom. These manuals are intended to answer most questions for new users, as well as provide valuable detailed information for the experienced user. The existence of detailed user s manuals will enable new users to become proficient with the codes, as well as reducing the dependency of new users on the code authors. In order to achieve the objectives listed, the following tasks were accomplished: 1) Parametric Study Of Preconditioning Parameters And Other Code Inputs; 2) Code Modifications To Reduce Runtimes; 3) Investigation Of Compiler Options To Reduce Code Runtime; and 4) Development/Enhancement of Users Manuals for Aardvark and Phantom

Marquart, Jed E.↗
NASA NTRS1987

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗
NASA NTRS1994

Enhanced simulation software for rocket turbopump, turbulent, annular liquid seals

One of the main objectives of this work is to develop a new dynamic analysis for liquid annular seals with arbitrary profile and to analyze a general distorted interstage seal of the space shuttle main engine high pressure oxygen turbopump (SSME-ATD-HPOTP). The dynamic analysis developed is based on a method originally proposed by Nelson and Nguyen. A simpler scheme based on cubic splines is found to be computationally more efficient and has better convergence properties at higher eccentricities. The first order solution of the original analysis is modified by including a more exact solution that takes into account the variation of perturbed variables along the circumference. A new set of equations for dynamic analysis are derived based on this more general model. A unified solution procedure that is valid for both Moody's and Hirs' friction models is presented. Dynamic analysis is developed for three different models: constant properties, variable properties, and thermal effects with variable properties. Arbitrarily varying seal profiles in both axial and circumferential directions are considered. An example case of an elliptical seal with varying degrees of axial curvature is analyzed in detail. A case study based on predicted clearances of an interstage seal of the SSME-ATD-HPOTP is presented. Dynamic coefficients based on external specified load are introduced to analyze seals that support a preload. The other objective of this work is to study the effect of large rotor displacements of SSME-ATD-HPOTP on the dynamics of the annular seal and the resulting transient motion. One task is to identify the magnitude of motion of the rotor about the centered position and establish limits of effectiveness of using current linear models. This task is accomplished by solving the bulk flow model seal governing equations directly for transient seal forces for any given type of motion, including motion with large eccentricities. Based on the above study, an equivalence is established between linearized coefficients based transient motion and the same motion as predicted by the original governing equations. An innovative method is developed to model nonlinearities in an annular seal based on dynamic coefficients computed at various static eccentricities. This method is thoroughly tested for various types of transient motion using bulk flow model results as a benchmark.

Padavala, Satya↗
NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗
NASA NTRS2017

Statistical Analysis of CFD Solutions from the 6th AIAA CFD Drag Prediction Workshop

A graphical framework is used for statistical analysis of the results from an extensive N-version test of a collection of Reynolds-averaged Navier-Stokes computational fluid dynamics codes. The solutions were obtained by code developers and users from North America, Europe, Asia, and South America using both common and custom grid sequences as well as multiple turbulence models for the June 2016 6th AIAA CFD Drag Prediction Workshop sponsored by the AIAA Applied Aerodynamics Technical Committee. The aerodynamic configuration for this workshop was the Common Research Model subsonic transport wing-body previously used for both the 4th and 5th Drag Prediction Workshops. This work continues the statistical analysis begun in the earlier workshops and compares the results from the grid convergence study of the most recent workshop with previous workshops.

Statistical analysis↗
NASA NTRS1996

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗
NASA NTRS1989

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix, which are updated only in the direction of the last move. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗
NASA NTRS2010

Statistical Analysis of CFD Solutions from the Fourth AIAA Drag Prediction Workshop

A graphical framework is used for statistical analysis of the results from an extensive N-version test of a collection of Reynolds-averaged Navier-Stokes computational fluid dynamics codes. The solutions were obtained by code developers and users from the U.S., Europe, Asia, and Russia using a variety of grid systems and turbulence models for the June 2009 4th Drag Prediction Workshop sponsored by the AIAA Applied Aerodynamics Technical Committee. The aerodynamic configuration for this workshop was a new subsonic transport model, the Common Research Model, designed using a modern approach for the wing and included a horizontal tail. The fourth workshop focused on the prediction of both absolute and incremental drag levels for wing-body and wing-body-horizontal tail configurations. This work continues the statistical analysis begun in the earlier workshops and compares the results from the grid convergence study of the most recent workshop with earlier workshops using the statistical framework.

Statistical analysis↗