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 955 records · Page 53

Photo sensor array technology development

The development of an improved capability photo sensor array imager for use in a Viking '75 type facsimile camera is presented. This imager consists of silicon photodiodes and lead sulfide detectors to cover a spectral range from 0.4 to 2.7 microns. An optical design specifying filter configurations and convergence angles is described. Three electronics design approaches: AC-chopped light, DC-dual detector, and DC-single detector, are investigated. Experimental and calculated results are compared whenever possible using breadboard testing and tolerance analysis techniques. Results show that any design used must be forgiving of the relative instability of lead sulfide detectors. A final design using lead sulfide detectors and associated electronics is implemented by fabrication of a hybrid prototype device. Test results of this device show a good agreement with calculated values.

Rossman, M. W.↗

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The time-dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard-wall duct without flow and with spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the steady solutions of the acoustic field. The transient method did not converge to the steady-state solution for cutoff acoustic duct modes. This has implications as to its use in a variable-area duct, where modes may become cutoff in the smal-area portion of the duct. For single cutoff mode propagation the steady-state impedance boundary condition produced acoustic reflections during the initial transient that caused finite instabilities in the numerical calculations. The stability problem is resolved by reformulating the exit boundary condition. Example calculations show good agreement with exact analytical and numerical results for forcing frequencies above, below, and nearly at the cutoff frequency.

Baumeister, K. J.↗

A new methodology for free wake analysis using curved vortex elements

A method using curved vortex elements was developed for helicopter rotor free wake calculations. The Basic Curve Vortex Element (BCVE) is derived from the approximate Biot-Savart integration for a parabolic arc filament. When used in conjunction with a scheme to fit the elements along a vortex filament contour, this method has a significant advantage in overall accuracy and efficiency when compared to the traditional straight-line element approach. A theoretical and numerical analysis shows that free wake flows involving close interactions between filaments should utilize curved vortex elements in order to guarantee a consistent level of accuracy. The curved element method was implemented into a forward flight free wake analysis, featuring an adaptive far wake model that utilizes free wake information to extend the vortex filaments beyond the free wake regions. The curved vortex element free wake, coupled with this far wake model, exhibited rapid convergence, even in regions where the free wake and far wake turns are interlaced. Sample calculations are presented for tip vortex motion at various advance ratios for single and multiple blade rotors. Cross-flow plots reveal that the overall downstream wake flow resembles a trailing vortex pair. A preliminary assessment shows that the rotor downwash field is insensitive to element size, even for relatively large curved elements.

Bliss, Donald B.↗

Hierarchic models for laminated plates

Structural plates and shells are three-dimensional bodies, one dimension of which happens to be much smaller than the other two. Thus, the quality of a plate or shell model must be judged on the basis of how well its exact solution approximates the corresponding three-dimensional problem. Of course, the exact solution depends not only on the choice of the model but also on the topology, material properties, loading and constraints. The desired degree of approximation depends on the analyst's goals in performing the analysis. For these reasons models have to be chosen adaptively. Hierarchic sequences of models make adaptive selection of the model which is best suited for the purposes of a particular analysis possible. The principles governing the formulation of hierarchic models for laminated plates are presented. The essential features of the hierarchic models described models are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. The formulation is based on one parameter (beta) which characterizes the hierarchic sequence of models, and a set of constants whose influence was assessed by a numerical sensitivity study. The recommended selection of these constants results in the number of fields increasing by three for each increment in the power of beta. Numerical examples analyzed with the proposed sequence of models are included and good correlation with the reference solutions was found. Results were obtained for laminated strips (plates in cylindrical bending) and for square and rectangular plates with uniform loading and with homogeneous boundary conditions. Cross-ply and angle-ply laminates were evaluated and the results compared with those of MSC/PROBE. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Actis, Ricardo Luis↗

Real-Time Parameter Estimation in the Frequency Domain

A method for real-time estimation of parameters in a linear dynamic state-space model was developed and studied. The application is aircraft dynamic model parameter estimation from measured data in flight. Equation error in the frequency domain was used with a recursive Fourier transform for the real-time data analysis. Linear and nonlinear simulation examples and flight test data from the F-18 High Alpha Research Vehicle were used to demonstrate that the technique produces accurate model parameter estimates with appropriate error bounds. Parameter estimates converged in less than one cycle of the dominant dynamic mode, using no a priori information, with control surface inputs measured in flight during ordinary piloted maneuvers. The real-time parameter estimation method has low computational requirements and could be implemented

Morelli, Eugene A.↗

The Energy Budget of the Polar Atmosphere in MERRA

Components of the atmospheric energy budget from the Modern Era Retrospective-analysis for Research and Applications (MERRA) are evaluated in polar regions for the period 1979-2005 and compared with previous estimates, in situ observations, and contemporary reanalyses. Closure of the energy budget is reflected by the analysis increments term, which results from virtual enthalpy and latent heating contributions and averages -11 W/sq m over the north polar cap and -22 W/sq m over the south polar cap. Total energy tendency and energy convergence terms from MERRA agree closely with previous study for northern high latitudes but convergence exceeds previous estimates for the south polar cap by 46 percent. Discrepancies with the Southern Hemisphere transport are largest in autumn and may be related to differences in topography with earlier reanalyses. For the Arctic, differences between MERRA and other sources in TOA and surface radiative fluxes maximize in May. These differences are concurrent with the largest discrepancies between MERRA parameterized and observed surface albedo. For May, in situ observations of the upwelling shortwave flux in the Arctic are 80 W/sq m larger than MERRA, while the MERRA downwelling longwave flux is underestimated by 12 W/sq m throughout the year. Over grounded ice sheets, the annual mean net surface energy flux in MERRA is erroneously non-zero. Contemporary reanalyses from the Climate Forecast Center (CFSR) and the Interim Re-Analyses of the European Centre for Medium Range Weather Forecasts (ERA-I) are found to have better surface parameterizations, however these collections are also found to have significant discrepancies with observed surface and TOA energy fluxes. Discrepancies among available reanalyses underscore the challenge of reproducing credible estimates of the atmospheric energy budget in polar regions.

Cullather, Richard I.↗

TEM Analyses of Itokawa Regolith Grains and Lunar Soil Grains to Directly Determine Space Weathering Rates on Airless Bodies

Samples returned from the moon and Asteroid Itokawa by NASA's Apollo Missions and JAXA's Hayabusa Mission, respectively, provide a unique record of their interaction with the space environment. Space weathering effects result from micrometeorite impact activity and interactions with the solar wind. While the effects of solar wind interactions, ion implantation and solar flare particle track accumulation, have been studied extensively, the rate at which these effects accumulate in samples on airless bodies has not been conclusively determined. Results of numerical modeling and experimental simulations do not converge with observations from natural samples. We measured track densities and rim thicknesses of three olivine grains from Itokawa and multiple olivine and anorthite grains from lunar soils of varying exposure ages. Samples were prepared for analysis using a Leica EM UC6 ultramicrotome and an FEI Quanta 3D dual beam focused ion beam scanning electron microscope (FIB-SEM). Transmission electron microscope (TEM) analyses were performed on the JEOL 2500SE 200kV field emission STEM. The solar wind damaged rims on lunar anorthite grains are amorphous, lack inclusions, and are compositionally similar to the host grain. The rim width increases as a smooth function of exposure age until it levels off at approximately 180 nm after approximately 20 My (Fig. 1). While solar wind ion damage can only accumulate while the grain is in a direct line of sight to the Sun, solar flare particles can penetrate to mm-depths. To assess whether the track density accurately predicts surface exposure, we measured the rim width and track density in olivine and anorthite from the surface of rock 64455, which was never buried and has a surface exposure age of 2 My based on isotopic measurements. The rim width from 64455 (60-70nm) plots within error of the well-defined trend for solar wind amorphized rims in Fig. 1. Measured solar flare track densities are accurately reflecting the surface exposure of the grains. Track densities correlate with the amorphous rim thicknesses. While the space-weathered rims of anorthite grains are amorphous, the space-weathered rims on both Itokawa and lunar olivine grains show solar wind damaged rims that are not amorphous. Instead, the rims are nanocrystalline with high dislocation densities and sparse inclusions of nanophase Fe metal. The rim thicknesses on the olivine grains also correlate with track density. The Itokawa olivine grains have track densities that indicate surface exposures of approximately 10(exp 5) years. Longer exposures (up to approximately 10(exp 7) years) do not amorphize the rims, as evidenced by lunar soil olivines with high track densities (approximately 10(exp 11) cm(exp -2)). From the combined data, shown in Fig. 1, it is clear that olivine is damaged (but not amorphized) more rapidly by the solar wind compared to anorthite. The olivine damaged rim forms quickly (in approximately 10(exp 6) y) and saturates at approximately 120nm with longer exposure time. The anorthite damaged rims form more slowly, amorphize, and grow thicker than the olivine rims. This is in agreement with numerical modeling data which predicts that solar wind damaged rims on anorthite will be thicker than olivine. However, the models predict that both olivine and anorthite rims will amorphize and reach equilibrium widths in less than 10(exp 3) y, in contrast to what is observed for natural samples. Laboratory irradiation experiments, which show rapid formation of fully amorphous and blistered surfaces from simulated solar wind exposures are also in contrast to observations of natural samples. These results suggest that there is a flux dependence on the type and extent of irradiation damage that develops in olivine. This flux dependence suggests that great caution be used in extrapolating between high-flux laboratory experiments and the natural case, as demonstrated by. We constrain the space weathering rate through analysis of returned samples. Provided that the track densities and the solar wind damaged rim widths exhibited by the Itokawa grains are typical of the fine-grained regions of Itokawa, then the space weathering rate is on the order of 10(exp 5) y. Space weathering effects in lunar soils saturate within a few My of exposure while those in Itokawa regolith grains formed in approximately 10(exp 5) y. Olivine and anorthite respond differently to solar wind irradiation. The space weathering effects in olivine are particularly difficult to reconcile with laboratory irradiation studies and numerical models. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind amorphization of different minerals on airless bodies.

Berger, Eve L.↗

Multilevel Monte Carlo Estimation of Unbiased Expectation via Sample Reuse and the Low Variance Estimation of Asymptotic Rates

A new variant of the multilevel Monte Carlo estimator [5, 3, 9, 12] is presented for the estimation of expectation statistics that utilizes sample reuse in specified levels, explicitly removes approximation error bias associated with numerically computed output quantities of interest that have an asymptotic limit behavior, and permits a low variance estimate of the asymptotic rate of convergence to that limit. In addition, it is shown that this new multilevel Monte Carlo variant can yield a computational cost savings. A review of Monte Carlo and multilevel Monte Carlo estimators is presented that includes analysis of expected value, expected mean squared error, and the calculation of optimized multilevel sample size parameters. The multilevel Monte Carlo estimator produces estimates of expectation for numerically approximated output quantities of interest that are biased by approximation error. When the quantity of interest can be modeled as the asymptotic limit of numerically approximated output quantities of interest, it is theoretically possible to remove this approximation error bias in the multilevel Monte Carlo estimator. In actual implementations, however, this procedure is unreliable due to statistical variability and inaccuracy in estimating the needed asymptotic limit. Analysis and numerical experiment show that the proposed variant of the multilevel Monte Carlo method greatly reduces (in some cases eliminates) the statistical variability in this limit estimation.

Barth, Timothy↗

Enhancing ACPF Analysis: Integrating Newton-Raphson Method with Gradient Descent and Computational Graphs

This paper presents a new method for enhancing Alternating Current Power Flow (ACPF) analysis. The method integrates the Newton-Raphson (NR) method with Enhanced-Gradient Descent (GD) and computational graphs. The integration of renewable energy sources in power systems introduces variability and unpredictability, and this method addresses these challenges. It leverages the robustness of NR for accurate approximations and the flexibility of GD for handling variable conditions, all without requiring Jacobian matrix inversion. Furthermore, computational graphs provide a structured and visual framework that simplifies and systematizes the application of these methods. The goal of this fusion is to overcome the limitations of traditional ACPF methods and improve the resilience, adaptability, and efficiency of modern power grid analyses. We validate the effectiveness of our advanced algorithm through comprehensive testing on established IEEE benchmark systems. Furthermore, our findings demonstrate that our approach not only speeds up the convergence process but also ensures consistent performance across diverse system states, representing a significant advancement in power flow computation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A comparison between residential relocation timing of Sydney and Chicago residents: A Bayesian survival analysis

We report that understanding households' behaviour in residential relocation timing is of great importance in the field of transport engineering and economics. This research aims to develop a residential relocation model by considering the potential dynamic impacts of other households' decisions and variables, including economic and demographic attributes, housing features, intra-household decision-making structures, travel mode choice, and other life-course attributes. A multivariate parametric survival model with both fixed and time-varying covariates is developed. To the best of the authors' knowledge, this study is the first paper in the literature of residential relocation timing to propose the use of a Bayesian model in contrast to the widely used classic frequentist approach and have conducted a discussion on its advantages. An emerging residential relocation dataset collected for two cities in Australia and the USA (Sydney and Chicago cities) has been used, which covers residence, vehicle ownership, occupation, education, economic and demographic attributes of respondents. A comprehensive comparison between the results of two cities and a comparison between two Bayesian and frequentist approaches are made. This study confirms the impact of life-course variables, intra-household decision-making behaviours, and sociodemographic attributes on home mobility. According to the model outputs, the accelerating or decelerating impact of explanatory variables on the relocation timing has been almost the same in the two cities. The Bayesian model was confirmed to have some advantages over the frequentist model, including being straightforward to interpret, availability of making inferences on the results, and ease of handling complex models, and optimisation convergence complexities.

99 GENERAL AND MISCELLANEOUS↗

Unsteady analysis of rotor blade tip flow

The development of the VSAERO-TS and VSAERO-H computer programs for calculating the unsteady aerodynamic characteristics of arbitrarily shaped wings oscillating in pitch is presented. The effect of several wake parameters on chordwise pressure distribution in VSAERO-TS is given and the convergence characteristics of both programs are discussed. In the program, the influence coefficient for each panel is formulated for a planar surface and so a skewed panel is represented by a projected flat quadrilateral lying in the mean plane. Since panels in the extreme roll-up region of the tip vortex are highly skewed the program was modified to treat each highly skewed panel as a pair of triangles. The programs are validated by comparing the chordwise pressure distribution of several blade tip planforms with experimental data. The comparison, for the most part, is good. The triangular panel representation improved the chordwise pressure distribution near the tip region for higher mean angle of attack.

Maskey, B.↗

An adaptive remeshing method for finite element thermal analysis

A finite element remeshing approach that makes use of quadrilateral and triangular elements is described. The approach uses the solution on a previous mesh to create a new mesh. Meshes are completely unstructured with highly refined elements in regions of steep gradients and larger elements where gradients are smaller. Studies of convergence rates for heat conduction problems with exact solutions show that for problems with highly localized solution variations, the remeshing approach gives smaller solution errors with fewer unknowns than refinement of uniform, structured meshes.

Thornton, Earl A.↗

A Six-Node Curved Triangular Element and a Four-Node Quadrilateral Element for Analysis of Laminated Composite Aerospace Structures

Mathematical development and some computed results are presented for Mindlin plate and shell elements, suitable for analysis of laminated composite and sandwich structures. These elements use the conventional 3 (plate) or 5 (shell) nodal degrees of freedom, have no communicable mechanisms, have no spurious shear energy (no shear locking), have no spurious membrane energy (no membrane locking) and do not require arbitrary reduction of out-of-plane shear moduli or under-integration. Artificial out-of-plane rotational stiffnesses are added at the element level to avoid convergence problems or singularity due to flat spots in shells. This report discusses a 6-node curved triangular element and a 4-node quadrilateral element. Findings show that in regular rectangular meshes, the Martin-Breiner 6-node triangular curved shell (MB6) is approximately equivalent to the conventional 8-node quadrilateral with integration. The 4-node quadrilateral (MB4) has very good accuracy for a 4-node element, and may be preferred in vibration analysis because of narrower bandwidth. The mathematical developments used in these elements, those discussed in the seven appendices, have been applied to elements with 3, 4, 6, and 10 nodes and can be applied to other nodal configurations.

Martin, C. Wayne↗

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

Steady-state off-design modeling of the supercritical carbon dioxide recompression cycle for concentrating solar power applications with two-tank sensible-heat storage

Concentrating solar power researchers are evaluating the potential of the supercritical carbon dioxide recompression cycle to improve the thermal efficiency and decrease the capital costs of next-generation systems. This analysis investigates the steady-state off-design performance of a recompression cycle integrated with two-tank sensible-heat thermal energy storage as the ambient temperature and heat-transfer fluid (HTF) inlet conditions to the cycle change. This paper presents off-design component models and then cycle convergence and control models to maximize net cycle power output while constraining the high-side pressure and air-cooler fan power to their respective design values and fixing the cycle HTF outlet temperature to its design value. Results show that inventory control and air-cooler fan power are important control parameters that can be optimized to maximize off-design power output as the ambient temperature and HTF mass flow rate diverge from design. The high-side pressure and fan power constraints cause the cycle net power to degrade when the ambient temperature is warmer than the design value, and this study calculates a maximum mass flow rate for hot days above which the cycle cannot achieve the design HTF outlet temperature. Finally, the analysis shows that optimizing compressor shaft speeds can improve cold-day performance by around 1.5 percentage points and part-load performance by up to 2.5 percentage points versus a baseline case with no active compressor control.

14 SOLAR ENERGY↗

Cubic AgMnSbTe 3 Semiconductor with a High Thermoelectric Performance

The reaction of MnTe with AgSbTe 2 in an equimolar ratio (ATMS) provides a new semiconductor, AgMnSbTe 3 . AgMnSbTe 3 crystallizes in an average rock-salt NaCl structure with Ag, Mn, and Sb cations statistically occupying the Na sites. AgMnSbTe 3 is a p-type semiconductor with a narrow band gap of ~0.33 eV. A pair distribution function analysis indicates that local distortions are associated with the location of the Ag atoms in the lattice. Density functional theory calculations suggest a specific electronic band structure with multi-peak valence band maxima prone to energy convergence. In addition, Ag 2 Te nanograins precipitate at grain boundaries of AgMnSbTe 3 . Here, the energy offset of the valance band edge between AgMnSbTe 3 and Ag 2 Te is ~0.05 eV, which implies that Ag 2 Te precipitates exhibit a negligible effect on the hole transmission. As a result, ATMS exhibits a high power factor of ~12.2 μWcm -1 K -2 at 823 K, ultralow lattice thermal conductivity of ~0.34 Wm -1 K -1 (823 K), high peak ZT of ~1.46 at 823 K, and high average ZT of ~0.87 in the temperature range of 400–823 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗

Voltage Calculations in Secondary Distribution Networks via Physics-Inspired Neural Network Using Smart Meter Data

The increasing penetration of distributed energy resources (DERs) leads to voltage issues across distribution networks, necessitating voltage calculations by utilities. Electric model-free voltage calculation offers an enticing solution. However, most researches mainly focus on primary distribution networks ignoring secondary distribution networks and commonly overlook extreme voltage case calculations, which require the model’s extrapolation abilities. Here, in addressing the gaps, this paper presents a customized physics-inspired neural network (PINN) model, the structure of which is inspired by the derived coupled power flow model of primary-secondary distribution networks. To ensure precision and rapid convergence, a crafted training framework for the PINN model is proposed. The PINN’s “structure-mimetic” design enables superior extrapolation for unseen scenarios and enhances physical information awareness. We demonstrate this through two applications: hosting capacity analysis and customer-transformer connectivity. The effectiveness and advantages of the proposed PINN model are validated on two public testing systems and one utility distribution feeder model.

Distribution network↗