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 109 records · Page 6

Optimal Control using Composite Bernstein Approximants

In this work, we present composite Bernstein polynomials as a direct collocation method for approximating optimal control problems. An analysis of the convergence properties of composite Bernstein polynomials is provided, and beneficial properties of composite Bernstein polynomials for the solution of optimal control problems are discussed. The efficacy of the proposed approximation method is demonstrated through a bang-bang example. Lastly, we apply this method to a motion planning problem, offering a practical solution that emphasizes the ability of this method to solve complex optimal control problems.

Gage MacLin↗

The proximity of hotspots to convergent and divergent plate boundaries

An analysis of four different hotspot distributions, ranging from Morgan's (1972) original list of 19 to Vogt's (1981) list of 117 reveals that the hotspots are preferentially located near divergent plate boundaries. The probability of this proximity occurring by chance alone is quite remote, less than 0.01 for all four hotspot distributions. The same analysis also reveals that the hotspots are preferentially excluded from regions near convergent plate boundaries. The probability of this exclusion occurring by chance alone is 0.1 or less for three out of the four distributions examined. We interpret this behavior as being a consequence of the effects of large scale convective circulation on ascending mantle plumes. Mantle thermal plumes, the most probable source of hotspots, arise from instabilities in a basal thermal boundary layer. Plumes are suppressed from regions beneath convergent boundaries by descending flow and are entrained into the upwelling flow beneath spreading centers. Plate-scale convective circulation driven by subduction may also advect mantle thermal plumes toward spreading centers.

Weinstein, Stuart A.↗

Recent advances in the structual analysis of scientific balloons

The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

Schur, W. W.↗

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

Formal verification of a fault tolerant clock synchronization algorithm

A formal specification and mechanically assisted verification of the interactive convergence clock synchronization algorithm of Lamport and Melliar-Smith is described. Several technical flaws in the analysis given by Lamport and Melliar-Smith were discovered, even though their presentation is unusally precise and detailed. It seems that these flaws were not detected by informal peer scrutiny. The flaws are discussed and a revised presentation of the analysis is given that not only corrects the flaws but is also more precise and easier to follow. Some of the corrections to the flaws require slight modifications to the original assumptions underlying the algorithm and to the constraints on its parameters, and thus change the external specifications of the algorithm. The formal analysis of the interactive convergence clock synchronization algorithm was performed using the Enhanced Hierarchical Development Methodology (EHDM) formal specification and verification environment. This application of EHDM provides a demonstration of some of the capabilities of the system.

Rushby, John↗

An extended field method solution for a pulse loaded thin plate

This paper deals with a new extension of a weighted residual method of analysis called the extended field method. The extended field method is applied for the first time to the problem of the transient vibration of a uniformly thin elastic plate. Numerical results have been obtained which validate the analysis procedure and show better solution convergence than is obtainable by standard methods of analysis for the same number of degrees of freedom. Further studies are necessary to obtain still better convergence.

Donaldson, B. K.↗

A zonal method for modeling powered-lift aircraft flow fields

A zonal method for modeling powered-lift aircraft flow fields is based on the coupling of a three-dimensional Navier-Stokes code to a potential flow code. By minimizing the extent of the viscous Navier-Stokes zones the zonal method can be a cost effective flow analysis tool. The successful coupling of the zonal solutions provides the viscous/inviscid interations that are necessary to achieve convergent and unique overall solutions. The feasibility of coupling the two vastly different codes is demonstrated. The interzone boundaries were overlapped to facilitate the passing of boundary condition information between the codes. Routines were developed to extract the normal velocity boundary conditions for the potential flow zone from the viscous zone solution. Similarly, the velocity vector direction along with the total conditions were obtained from the potential flow solution to provide boundary conditions for the Navier-Stokes solution. Studies were conducted to determine the influence of the overlap of the interzone boundaries and the convergence of the zonal solutions on the convergence of the overall solution. The zonal method was applied to a jet impingement problem to model the suckdown effect that results from the entrainment of the inviscid zone flow by the viscous zone jet. The resultant potential flow solution created a lower pressure on the base of the vehicle which produces the suckdown load. The feasibility of the zonal method was demonstrated. By enhancing the Navier-Stokes code for powered-lift flow fields and optimizing the convergence of the coupled analysis a practical flow analysis tool will result.

Roberts, D. W.↗

Improving the convergence rate to steady state of parabolic ADI methods

The present, residuals' L(2)-norms analysis of the rate of convergence to steady state for parabolic ADI solvers allows the prediction of the number of iterations required for convergence, as a function of the Courant number alpha. A modification of current ADI codes is presented which significantly improves the convergence rate and is insensitive to the Courant number over a large range of alpha. This corrected algorithm is tested for the cases of Dirichlet problems for uniform grids of many mesh sizes, mixed Dirichlet-Neumann problems, and problems defined on stretched grids and/or problems with variable coefficients.

Abarbanel, Saul S.↗

Spectroscopic analysis of transition state energy levels - Bending-rotational spectrum and lifetime analysis of H3 quasibound states

Converged quantum mechanical calculations of scattering matrices and transition probabilities are reported for the reaction of H with H2 with total angular momentum 0, 1, and 4 as functions of total energy in the range 0.85-1.15 eV on an accurate potential energy surface. The resonance structure is illustrated with Argand diagrams. State-to-state reactive collision delay times and lifetimes are presented. For J = 0, 1, and 4, the lowest-energy H3 resonance is at total energies of 0.983, 0.985, and 1.01 eV, respectively, with lifetimes of about 16-17 fs. For J = 1 and 4 there is a higher-energy resonance at 1.10-1.11 eV. For J = 1 the lifetime is about 4 fs and for J = 4 it is about 1 fs.

Zhao, Meishan↗

Severe storm initiation and development from satellite infrared imagery and Rawinsonde data

The geographical distribution of potential temperatures, mixing ratio, and streamlines of flow patterns at 850, 700, and 500 mb heights are used to understand the prestorm convection and the horizontal convergence of moisture. From the analysis of 21 tornadoes the following conclusions are reached: (1) Strong horizontal convergence of moisture appeared at the 850, 700, and 500 mb levels in the area 12 hours before the storm formation; (2) An abundantly moist atmosphere below 3 km (700 mb) becomes convectively unstable during the time period between 12 and 24 hours before the initiation of the severe storms; (3) Strong winds veering with height with direction parallel to the movement of a dryline, surface fronts, etc; (4) During a 36-hour period, a tropopause height in the areas of interest is lowest at the time of tornadic cloud formation; (5) A train of gravity waves is detected before and during the cloud formation period. Rapid-scan infrared imagery provides near real-time information on the life cycle of the storm which can be summarized as follows: (1) Enhanced convection produced an overshooting cloud top penetrating above the tropopause, making the mass density of the overshooting cloud much greater than the mass density of the surrounding air; (2) The overshooting cloud top collapsed at the end of the mature stage of the cloud development; (3) The tornado touchdown followed the collapse of the overshooting cloud top.

Hung, R. J.↗

Wavenumber-frequency deconvolution of aeroacoustic microphone phased array data of arbitrary coherence

Deconvolution of aeroacoustic data acquired with microphone phased arrays is a computationally-challenging task for distributed sources with arbitrary coherence. A new technique for performing such deconvolution is proposed. This technique relies on analysis of the array data in the wavenumber-frequency domain, allowing for fast convolution and reduced storage requirements when compared to traditional coherent deconvolution. A positive semidefinite constraint for the iterative deconvolution procedure is implemented and shows improved behavior in terms of quantifiable convergence metrics when compared to a standalone covariance inequality constraint. A series of simulations validates the method's ability to resolve coherence and phase angle relationships between partially-coherent sources, as well as determines convergence criteria for deconvolution analysis. Simulations for point sources near the microphone phased array show potential for handling such data in the wavenumber-frequency domain. In particular, a physics-based integration boundary calculation is described, and can successfully isolate sources and track the appropriate integration bounds with and without the presence of flow. Magnitude and phase relationships between multiple sources are successfully extracted. Limitations of the deconvolution technique are determined from the simulations, particularly in the context of a simulated acoustic field in a closed test section wind tunnel with strong boundary layer contamination. A final application to a trailing edge noise experiment conducted in an open-jet wind tunnel matches best estimates of acoustic levels from traditional calculation methods and qualitatively assesses the coherence characteristics of the trailing edge noise source.

phased array↗

Tectonics and volcanism of Eastern Aphrodite Terra, Venus - No subduction, no spreading

Eastern Aphrodite Terra, a deformed region with high topographic relief on Venus, has been interpreted as analogous to a terrestrial extensional or convergent plate boundary. However, analysis of geological and structural relations indicates that the tectonics of eastern Aphrodite Terra is dominated by blistering of the crust by magma diapirs. The findings imply that, within this region, vertical tectonism dominates over horizontal tectonism and, consequently, that this region is neither a divergent nor a convergent plate boundary.

Hansen, Vicki L.↗

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗

Thermal analysis and design of a cooling system for a Mach 14 nozzle

The analysis and design of a Mach 14 converging diverging nozzle wall liner is provided. The analysis indicates that: no fin on the coolant side of the nozzle wall is optimum, the thermal stresses are dominant, and the critical area is very near the throat. The molybdenum alloy TZM, with a wall thickness of 2.0 mm in the throat area, appears to be the only material capable of meeting design requirements. Additionally, cooling water at 2000 psia with a flow velocity of 25 m/s in the coolant passages is required.

Mullisen, Ronald↗

Quantitative validation of Ames IR intensity and new line lists for (32/33/34)S(16)O2, (32)S(18)O2 and (16)O(32)S(18)O

The quality of Ames-296K SO2 Infrared (IR) line list intensities is first validated by quantitative exploration of several dipole moment surfaces (DMSs) and partition sum convergence. The DMSs are computed with several of Dunning’s correlation-consistent basis sets and their vibrational dependence are compared to the empirical model derived from Stark effect experiments reported by D. Patel, D. Margolese, and T.R. Dykea [J.Chem.Phys. 70, 2740 (1979)]. The effective dipole deviations from the DMS adopted in the Ames IR lists is 0.2-0.4% for vibrational states up to 3v3. The vibrational dependence of the dipole moment is also in good agreement, except for nv1. Partition sum convergence at 296K is confirmed by new calculations with rotational quantum number J up to 150 and upper state E’ up to 8000 cm-1. The isotopologue consistency of the Ames IR line lists is superior relative to the regular Effective Hamiltonian (EH) models and Effective Dipole Moment (EDM) models. The v1+v2 and v2+v3 intensity consistency check reveals the recently reported experimental intensities need significant improvement or re-analysis. After the accuracy, convergence, and isotopologue consistency have been confirmed, the theoretical Ames-296K intensities are combined with the experimental line positions or EH models that experimental spectroscopists published after 2009. Three high-resolution IR line sets are reported for the (32/33/34)S(16)O2, (32)S(18)O2 and (16)O(32)S(18)O isotopologues: (1) the “New Lines Sets” include experimentally measured line positions; (2) the “Expanded Line Sets” include possible transitions among new rovibrational levels assigned in experiments and ground state (GS) levels predicted by reliable EH models; (3) the “Ames + MARVEL Sets” include possible transitions among all those levels reported in a recent MARVEL analysis. [Tóbiás et al, JQSRT 208, 152 (2018)]. Compared to the limited data in High-resolution TRANsmission molecular absorption database (HITRAN), these line sets have significantly improved the data coverage up to 4000 cm-1. Some missing bands can be traced to the unpublished experimental data. The isotopologue consistency of these line sets will help identify the uncertainties and defects in the experimental EH and EDM models. These line sets are good candidates for the next HITRAN update, if line shape parameters are available. The line sets can be downloaded from supplementary files or from the Ames Molecular Spectroscopic Database at http://huang.seti.org.

Xinchuan Huang↗

An all-at-once reduced Hessian SQP scheme for aerodynamic design optimization

This paper introduces a computational scheme for solving a class of aerodynamic design problems that can be posed as nonlinear equality constrained optimizations. The scheme treats the flow and design variables as independent variables, and solves the constrained optimization problem via reduced Hessian successive quadratic programming. It updates the design and flow variables simultaneously at each iteration and allows flow variables to be infeasible before convergence. The solution of an adjoint flow equation is never needed. In addition, a range space basis is chosen so that in a certain sense the 'cross term' ignored in reduced Hessian SQP methods is minimized. Numerical results for a nozzle design using the quasi-one-dimensional Euler equations show that this scheme is computationally efficient and robust. The computational cost of a typical nozzle design is only a fraction more than that of the corresponding analysis flow calculation. Superlinear convergence is also observed, which agrees with the theoretical properties of this scheme. All optimal solutions are obtained by starting far away from the final solution.

Feng, Dan↗

High Fidelity System Simulation of Multiple Components in Support of the UEET Program

The High Fidelity System Simulation effort has addressed various important objectives to enable additional capability within the NPSS framework. The scope emphasized High Pressure Turbine and High Pressure Compressor components. Initial effort was directed at developing and validating intermediate fidelity NPSS model using PD geometry and extended to high-fidelity NPSS model by overlaying detailed geometry to validate CFD against rig data. Both "feedforward" and feedback" approaches of analysis zooming was employed to enable system simulation capability in NPSS. These approaches have certain benefits and applicability in terms of specific applications "feedback" zooming allows the flow-up of information from high-fidelity analysis to be used to update the NPSS model results by forcing the NPSS solver to converge to high-fidelity analysis predictions. This apporach is effective in improving the accuracy of the NPSS model; however, it can only be used in circumstances where there is a clear physics-based strategy to flow up the high-fidelity analysis results to update the NPSS system model. "Feed-forward" zooming approach is more broadly useful in terms of enabling detailed analysis at early stages of design for a specified set of critical operating points and using these analysis results to drive design decisions early in the development process.

Plybon, Ronald C.↗

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗