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 199 records · Page 11

A quasi-Newton versus a homotopy method for nonlinear structural analysis

The globally convergent quasi-Newton minimization algorithm and the homotopy algorithms are discussed in detail and their effectiveness in solving certain classes of highly nonlinear problems of structural analysis is demonstrated. The application of the double dogleg strategy controls the directions and step-lengths of the quasi-Newtonian algorithm and overcomes the problem of nonpositive definite Hessians being produced during the iteration process. The algorithms are applied to a centrally loaded clamped beam, the snap-through of a shallow arch, and a shallow reticulated dome.

Kamat, M. P.↗

System Sensitivity Analysis Applied to the Conceptual Design of a Dual-Fuel Rocket SSTO

This paper reports the results of initial efforts to apply the System Sensitivity Analysis (SSA) optimization method to the conceptual design of a single-stage-to-orbit (SSTO) launch vehicle. SSA is an efficient, calculus-based MDO technique for generating sensitivity derivatives in a highly multidisciplinary design environment. The method has been successfully applied to conceptual aircraft design and has been proven to have advantages over traditional direct optimization methods. The method is applied to the optimization of an advanced, piloted SSTO design similar to vehicles currently being analyzed by NASA as possible replacements for the Space Shuttle. Powered by a derivative of the Russian RD-701 rocket engine, the vehicle employs a combination of hydrocarbon, hydrogen, and oxygen propellants. Three primary disciplines are included in the design - propulsion, performance, and weights & sizing. A complete, converged vehicle analysis depends on the use of three standalone conceptual analysis computer codes. Efforts to minimize vehicle dry (empty) weight are reported in this paper. The problem consists of six system-level design variables and one system-level constraint. Using SSA in a 'manual' fashion to generate gradient information, six system-level iterations were performed from each of two different starting points. The results showed a good pattern of convergence for both starting points. A discussion of the advantages and disadvantages of the method, possible areas of improvement, and future work is included.

Olds, John R.↗

Comparative genomic analysis of thermophilic fungi reveals convergent evolutionary adaptations and gene losses

Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.

59 BASIC BIOLOGICAL SCIENCES↗

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

An inexact semismooth Newton method with application to adaptive randomized sketching for dynamic optimization

In many applications, one can only access the inexact gradients and inexact hessian times vector products. Thus it is essential to consider algorithms that can handle such inexact quantities with a guaranteed convergence to solution. An inexact adaptive and provably convergent semismooth Newton method is considered to solve constrained optimization problems. In particular, dynamic optimization problems, which are known to be highly expensive, are the focus. A memory efficient semismooth Newton algorithm is introduced for these problems. The source of efficiency and inexactness is the randomized matrix sketching. Further, applications to optimization problems constrained by partial differential equations are also considered.

97 MATHEMATICS AND COMPUTING↗

Real-Time Local Volt/VAR Control Under External Disturbances with High PV Penetration

Volt/var control (VVC) of smart PV inverter is becoming one of the most popular solutions to address the voltage challenges associated with high PV penetration. This work focuses on the local droop VVC recommended by the grid integration standards IEEE1547, rule21 and addresses their major challenges i.e. appropriate parameters selection under changing conditions, and the control being vulnerable to instability (or voltage oscillations) and significant steady state error (SSE). This is achieved by proposing a two-layer local real-time adaptive VVC that has two major features i.e. a) it is able to ensure both low SSE and control stability simultaneously without compromising either; and b) it dynamically adapts its parameters to ensure good performance in a wide range of external disturbances such as sudden cloud cover, cloud intermittency, and substation voltage changes. A theoretical analysis and convergence proof of the proposed control is also discussed. The proposed control is implementation friendly as it fits well within the integration standard framework and depends only on the local bus information. Furthermore, the performance is compared with the existing droop VVC methods in several scenarios on a large unbalanced 3-phase feeder with detailed secondary side modeling.

14 SOLAR ENERGY↗

Parametric decay of an extraordinary electromagnetic wave in relativistic plasma

Parametric instability of an extraordinary electromagnetic wave in plasma preheated to a relativistic temperature is considered. A set of self-similar nonlinear differential equations taking into account the electron “thermal” mass is derived and investigated. Small perturbations of the parameters of the heated plasma are analyzed in the linear approximation by using the dispersion relation determining the phase velocities of the fast and slow extraordinary waves. In contrast to cold plasma, the evanescence zone in the frequency range above the electron upper hybrid frequency vanishes and the asymptotes of both branches converge. Theoretical analysis of the set of nonlinear equations shows that the growth rate of decay instability increases with increasing initial temperature of plasma electrons. This result is qualitatively confirmed by numerical simulations of plasma heating by a laser pulse injected from vacuum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗

Percept

The Percept software package is a collection of libraries and executables that provide tools for verifying computer simulations of engineering components and systems. Percept is useful for simulations using the finite element or finite volume methods on unstructured meshes. Percept includes API's for adaptive mesh refinement, geometry representation, the method of manufactured solutions, analysis of convergence including the convergence of vibrational eigenmodes, and metrics for analyzing the difference between fields represented on two different overlapping unstructured grids.

Carnes, Brian↗

Validity test for linear error analysis

To determine whether estimation process simulated by linear error analysis will converge, criterion has been developed based on extension of classical observability. Particular application of technique is with groups of batched navigation data where statistics of estimation errors are derived with classical minimum-variance methods.

Diamant, L. S.↗

Adaptive finite element analysis based on p-convergence

The results of numerical experiments are presented in which a posteriori estimators of error in strain energy were examined on the basis of a typical problem in linear elastic fracture mechanics. Two estimators were found to give close upper and lower bounds for the strain energy error. The potential significance of this is that the same estimators may provide a suitable basis for adaptive redistribution of the degrees of freedom in finite element models.

Szabo, B. A.↗

Prediction of subsonic aircraft flows with jet exhaust interactions

A numerical procedure to calculate the flow fields resulting from the viscous inviscid interactions that occur when a strong jet exhaust and aircraft flow field coupling exists was developed. The approach divides the interaction region into zones which are either predominantly viscous or inviscid. The flow in the inviscid zone, which surrounds most of the aircraft, is calculated using an existing potential flow code. The viscous flow zone, which encompasses the jet plume, is modeled using a parabolized Navier-Stokes code. The procedure features the coupling of the zonal solutions such that sufficient information is transferred between the zones to preserve the effects of the interactions. The zonal boundaries overlap and the boundary conditions are the information link between zones. An iteration scheme iterates the coupled analysis until convergence has been obtained.

Roberts, D. W.↗

Flow analysis of nozzle installations with strong airplane flow interactions

A numerical procedure has been developed to calculate the flow fields resulting from the viscous-inviscid interactions that occur when a strong jet exhaust and aircraft flow field coupling exists. The approach used in the current procedure is to divide the interaction region into zones which are either predominantly viscous or inviscid. The flow in the inviscid zone, which surrounds most of the aircraft, is calculated using an existing linearized potential flow code. The viscous flow zone, which encompasses the jet plume, is modeled using a parabolized Navier-Stokes code. The key feature of the present procedure is the coupling of the zonal solutions such that sufficient information is transferred between the zones to preserve the effects of the interactions. The zonal boundaries overlap with the boundary conditions being the information link between zones. An iteraction scheme iterates the coupled analysis until convergence has been obtained. The procedure has been successfully used for several test cases for which the computed results are presented.

Roberts, D. W.↗

A zonal method for modeling 3-D aircraft flow fields with jet plume effects

A coupled zonal method has been developed for modeling flows about aircraft with highly integrated advanced nozzle installations in which the exhaust plume interacts with the aircraft flow field to enhance performance. A cost effective flow analysis is achieved by dividing the flow domain into zones and coupling the zonal solutions. PAN AIR is used to model the flow in the inviscid zone which surrounds the aircraft. The plume zone flow field is solved with a PNS code. The flow in the optional internal nozzle zone is modeled using a full Navier-Stokes code. An automated iterative procedure couples the inviscid and plume zones by passing updated boundary conditions between the zones. The coupled analysis provides converged solutions for the overall flow field. The method is demonstrated for a highly integrated nonaxisymmetric vectored nozzle installed on a swept-wing aircraft.

Roberts, D. W.↗

Applications of a transonic wing design method

A method for designing wings and airfoils at transonic speeds using a predictor/corrector approach was developed. The procedure iterates between an aerodynamic code, which predicts the flow about a given geometry, and the design module, which compares the calculated and target pressure distributions and modifies the geometry using an algorithm that relates differences in pressure to a change in surface curvature. The modular nature of the design method makes it relatively simple to couple it to any analysis method. The iterative approach allows the design process and aerodynamic analysis to converge in parallel, significantly reducing the time required to reach a final design. Viscous and static aeroelastic effects can also be accounted for during the design or as a post-design correction. Results from several pilot design codes indicated that the method accurately reproduced pressure distributions as well as the coordinates of a given airfoil or wing by modifying an initial contour. The codes were applied to supercritical as well as conventional airfoils, forward- and aft-swept transport wings, and moderate-to-highly swept fighter wings. The design method was found to be robust and efficient, even for cases having fairly strong shocks.

Campbell, Richard L.↗

Percept

The Percept software package is a collection of libraries and executables that provide tools for verifying computer simulations of engineering components and systems. Percept is useful for simulations using the finite element or finite volume methods on unstructured meshes. Percept includes API's for adaptive mesh refinement, geometry representation, the method of manufactured solutions, analysis of convergence including the convergence of vibrational eigenmodes, and metrics for analyzing the difference between fields represented on two different overlapping unstructured grids.

Source record↗