Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling 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 19 records

Derivation of improved load transformation matrices for launchers-spacecraft coupled analysis, and direct computation of margins of safety

Load and stress recovery from transient dynamic studies are improved upon using an extended acceleration vector in the modal acceleration technique applied to structural analysis. Extension of the normal LTM (load transformation matrices) stress recovery to automatically compute margins of safety is presented with an application to the Hubble space telescope.

Klein, M.↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗

Shape optimization of single- and two-element airfoils on multiblock grids

A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Lacasse, James M.↗

An assessment of airfoil design by numerical optimization

A practical procedure for optimum design of aerodynamic shapes is demonstrated. The proposed procedure uses an optimization program based on the method of feasible directions coupled with an analysis program that uses a relaxation solution of the inviscid, transonic, small-disturbance equations. Results are presented for low-drag, nonlifting transonic airfoils. Extension of the method to lifting airfoils, other speed regimes, and to three dimensions if feasible.

Hicks, R. M.↗

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗

Implementation of a multiblock sensitivity analysis method in numerical aerodynamic shape optimization

A multiblock sensitivity analysis method is applied in a numerical aerodynamic shape optimization technique. The Sensitivity Analysis Domain Decomposition (SADD) scheme which is implemented in this study was developed to reduce the computer memory requirements resulting from the aerodynamic sensitivity analysis equations. Discrete sensitivity analysis offers the ability to compute quasi-analytical derivatives in a more efficient manner than traditional finite-difference methods, which tend to be computationally expensive and prone to inaccuracies. The direct optimization procedure couples CFD analysis based on the two-dimensional thin-layer Navier-Stokes equations with a gradient-based numerical optimization technique. The linking mechanism is the sensitivity equation derived from the CFD discretized flow equations, recast in adjoint form, and solved using direct matrix inversion techniques. This investigation is performed to demonstrate an aerodynamic shape optimization technique on a multiblock domain and its applicability to complex geometries. The objectives are accomplished by shape optimizing two aerodynamic configurations. First, the shape optimization of a transonic airfoil is performed to investigate the behavior of the method in highly nonlinear flows and the effect of different grid blocking strategies on the procedure. Secondly, shape optimization of a two-element configuration in subsonic flow is completed. Cases are presented for this configuration to demonstrate the effect of simultaneously reshaping interfering elements. The aerodynamic shape optimization is shown to produce supercritical type airfoils in the transonic flow from an initially symmetric airfoil. Multiblocking effects the path of optimization while providing similar results at the conclusion. Simultaneous reshaping of elements is shown to be more effective than individual element reshaping due to the inclusion of mutual interference effects.

Lacasse, James M.↗

Mode coupling in partial reflections from the ionosphere at vertical incidence.

The effect of mode coupling has been calculated for radio waves partially reflected from the lower ionosphere at vertical incidence. Calculations were made for parameters appropriate to the partial reflection experiments (Belrose, 1970). The four reflection coefficients (two direct and two coupled) are calculated for a plane discontinuity using a matrix perturbation analysis. The ratios of coupled to direct reflection coefficients are calculated for a general one-dimensional reflecting layer using the coupled equations of Clemmow and Heading (1954). At the highest electron density considered, 10,000 per cu cm, about 10% of the reflected extraordinary wave amplitude might be due to coupling from the incident ordinary wave. The effect is smaller at lower electron density.

Connolly, D. J.↗

Application of numerical optimization to the design of low speed airfoils

A practical procedure for the optimum design of low-speed airfoils is demonstrated. The procedure uses an optimization program based on the method of feasible directions coupled with an aerodynamic analysis program that uses a relaxation solution of the inviscid, full potential equation. Results are presented for airfoils designed to have small adverse pressure gradients, high maximum lift, and low pitching moment.

Hicks, R. M.↗

Computational optimization and wind tunnel test of transonic wing designs

A practical procedure for the optimum design of transonic wings is demonstrated. The procedure uses an optimization program based on the method of feasible directions coupled with an aerodynamic analysis program which solves the three-dimensional potential equation for subsonic through transonic flow. Two new wings for the A-7 aircraft were designed by using the optimization procedure to achieve specified surface pressure distributions. The new wings, along with the existing A-7 wing, were tested in the Ames 11-foot transonic wind tunnel. The experimental data show that all of the performance goals were met.

Haney, H. P.↗

Application of numerical optimization to the design of wings with specified pressure distributions

A practical procedure for the optimum design of transonic wings is demonstrated. The procedure uses an optimization program based on the method of feasible directions coupled with an aerodynamic analysis program which solves the three-dimensional potential equation for subsonic through transonic flow. Two new wings for the A-7 aircraft were designed by using the optimization procedure to achieve specified surface pressure distributions. The new wings, along with the existing A-7 wing, were tested in the Ames 11 ft transonic wind tunnel. The experimental data show that all of the performance goals were met. However, comparisons of the wind tunnel results with the theoretical predictions indicate some differences at conditions for which strong shock waves occur.

Haney, H. P.↗

Design optimization of axisymmetric bodies in nonuniform transonic flow

An inviscid transonic code capable of designing an axisymmetric body in a uniform or nonuniform flow was developed. The design was achieved by direct optimiation by coupling an analysis code with an optimizer. Design examples were provided for axisymmetric bodies with fineness ratios of 8.33 and 5 at different Mach numbers. It was shown that by reducing the nose radius and increasing the afterbody thickness of initial shapes obtained from symmetric NACA four-digit airfoil contours, wave drag could be reduced by 29 percent for a body of fineness ratio 8.33 in a nonuniform transonic flow of M = 0.98 to 0.995. The reduction was 41 percent for a body of fineness ratio 5 in a uniform transonic flow of M = 0.925 and 65 percent for the same body but in a nonuniform transonic flow of M = 0.90 to 0.95.

Lan, C. Edward↗

Efficient Two-Way Coupled Analysis of Steady-State Particle-Laden Hypersonic Flows

A direct solution approach for surface erosion in particle-laden hypersonic flows is extended for use in low-cost two-way coupled solutions of dilute gas-particle flows. The Trajectory Control Volume method, which uses a sparse set of probe particles to predict surface erosion distributions on general vehicles, is reformulated for the solution of source terms by mean trajectory subdivision and computing a flux differencing. The approach is verified successfully against a boundary layer solution and shown to agree well with experimental measurements. A representative Mars entry case, with conditions and geometry based on the ExoMars Schiaparelli capsule, is solved with the approach to study the impact of two-way coupling on surface heating and erosion. Results indicate that for realistic loading conditions, heating is largely unmodified compared to one-way coupled results at peak heating trajectory conditions, and no measureable difference is observed in the surface erosion rate. At exaggerated loading conditions high enough to observe coupling effects, the worst case collisional heating can increase heating by up to 60%.

Particle Laden Flow↗

Concurrent optimization of airframe and engine design parameters

An integrated system for the multidisciplinary analysis and optimization of airframe and propulsion design parameters is being developed. This system is known as IPAS, the Integrated Propulsion/Airframe Analysis System. The traditional method of analysis is one in which the propulsion system analysis is loosely coupled to the overall mission performance analysis. This results in a time consuming iterative process. First, the engine is designed and analyzed. Then, the results from this analysis are used in a mission analysis to determine the overall aircraft performance. The results from the mission analysis are used as a guide as the engine is redesigned and the entire process repeated. In IPAS, the propulsion system, airframe, and mission are closely coupled. The propulsion system analysis code is directly integrated into the mission analysis code. This allows the propulsion design parameters to be optimized along with the airframe and mission design parameters, significantly reducing the time required to obtain an optimized solution.

Lavelle, Thomas M.↗

An in-depth probabilistic study of external tank attach ring

This report deals with conducting a probabilistic study of the external tank attach ring (ETA) used as an interface between the external tank attach struts and the solid rocket booster. The ideas were to use probabilistic distributions for material, geometric, and load properties; to calculate probabilistic margins of safety; and then to compare results against the deterministic factors of safety that were used in the actual design process. The report describes how this was done and discusses some of the road blocks and data problems that were encountered during the study and provides some conclusions. A further refinement of this study is being considered for future work which would make more direct use of finite element analysis data coupled with Monte Carlo simulation. The basic conclusion herein indicates that the probabilistic margins of safety for the cases analyzed (by use of existing data) appear to support deterministic results and actually indicate higher reliabilities.

Pizzano, Frank↗

Concurrent optimization of airframe and engine design parameters

An integrated system for the multidisciplinary analysis and optimization of airframe and propulsion design parameters is being developed. This system is known as IPAS, the Integrated Propulsion/Airframe Analysis System. The traditional method of analysis is one in which the propulsion system analysis is loosely coupled to the overall mission performance analysis. This results in a time consuming iterative process. First, the engine is designed and analyzed. Then, the results from this analysis are used in a mission analysis to determine the overall aircraft performance. The results from the mission analysis are used as a guide as the engine redesigned and the entire process repeated. In IPAS, the propulsion system, airframe, and mission are closely coupled. The propulsion system analysis code is directly integrated into the mission analysis code. This allows the propulsion design parameters to be optimized along with the airframe and mission design parameters, significantly reducing the time required to obtain an optimized solution.

Lavelle, Thomas M.↗

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

Efficient sensitivity analysis and optimization of a helicopter rotor

Aeroelastic optimization of a system essentially consists of the determination of the optimum values of design variables which minimize the objective function and satisfy certain aeroelastic and geometric constraints. The process of aeroelastic optimization analysis is illustrated. To carry out aeroelastic optimization effectively, one needs a reliable analysis procedure to determine steady response and stability of a rotor system in forward flight. The rotor dynamic analysis used in the present study developed inhouse at the University of Maryland is based on finite elements in space and time. The analysis consists of two major phases: vehicle trim and rotor steady response (coupled trim analysis), and aeroelastic stability of the blade. For a reduction of helicopter vibration, the optimization process requires the sensitivity derivatives of the objective function and aeroelastic stability constraints. For this, the derivatives of steady response, hub loads and blade stability roots are calculated using a direct analytical approach. An automated optimization procedure is developed by coupling the rotor dynamic analysis, design sensitivity analysis and constrained optimization code CONMIN.

Lim, Joon W.↗