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 1,639 records · Page 91

Improved methods of vibration analysis of pretwisted, airfoil blades

Vibration analysis of pretwisted blades of asymmetric airfoil cross section is performed by using two mixed variational approaches. Numerical results obtained from these two methods are compared to those obtained from an improved finite difference method and also to those given by the ordinary finite difference method. The relative merits, convergence properties and accuracies of all four methods are studied and discussed. The effects of asymmetry and pretwist on natural frequencies and mode shapes are investigated. The improved finite difference method is shown to be far superior to the conventional finite difference method in several respects. Close lower bound solutions are provided by the improved finite difference method for untwisted blades with a relatively coarse mesh while the mixed methods have not indicated any specific bound.

Subrahmanyam, K. B.↗

Jacobi-Bessel analysis of reflector antennas with elliptical apertures

Although many reflector antennas possess circular projected apertures, there are recent satellite and ground antenna applications for which it is desirable to employ reflectors with elliptical apertures. Here a modification of the Jacobi-Bessel expansion is presented for the diffraction analysis of reflectors with elliptical apertures. A comparative study is also performed between this modified Jacobi-Bessel algorithm and the one which uses the Jacobi-Bessel expansion over a circumscribing circular region. Numerical results are presented for offset reflectors with elliptical and circular apertures and the improved convergence properties of the modified algorithm are highlighted.

Rahmat-Samii, Yahya↗

Quotient-difference type generalizations of the power method and their analysis

The recursion relations that were proposed by W. F. Ford and A. Sidi (Appl. Numer. Math, 4 (1988), pp. 477-489) for implementing vector extrapolation methods are used for devising generalizations of the power method for linear operators. These generalizations are shown to produce approximations to largest eigenvalues of a linear operator under certain conditions. They are similar in form to the quotient-difference algorithm and share similar convergence properties with the latter. These convergence properties also resemble those obtained for the basic LR and QR algorithms. Finally, it is shown that the convergence rate produced by one fo these generalizations is twice as fast for normal operators as it is for nonnormal operators.

Sidi, Avram↗

Towards industrial-strength Navier-Stokes codes

In this paper we discuss our experiences with Navier-Stokes (NS) codes using central differencing (CD) and scalar artificial dissipation (SAD). The NS-CDSAD codes have been developed by several researchers. Our results confirm that for typical commercial transport wing and wing/body configurations flying at transonic conditions with all turbulent boundary layers, NS-CDSAD codes, when used with the Johnson-King turbulence model, are capable of computing pressure distributions in excellent agreement with experimental data. However, results are not as good when laminar boundary layers are present. Exhaustive 2-D grid refinement studies supported by detailed analysis suggest that the numerical errors associated with SAD severely contaminate the solution in the laminar portion of the boundary layer. It is left as a challenge to the CFD community to find and fix the problems with Navier-Stokes codes and to produce a NS code which converges reliably and properly captures the laminar portion of the boundary layer on a reasonable grid.

Jou, Wen-Huei↗

Preconditioning for the Navier-Stokes equations with finite-rate chemistry

The extension of Van Leer's preconditioning procedure to generalized finite-rate chemistry is discussed. Application to viscous flow is begun with the proper preconditioning matrix for the one-dimensional Navier-Stokes equations. Eigenvalue stiffness is resolved and convergence-rate acceleration is demonstrated over the entire Mach-number range from nearly stagnant flow to hypersonic. Specific benefits are realized at the low and transonic flow speeds typical of complete propulsion-system simulations. The extended preconditioning matrix necessarily accounts for both thermal and chemical nonequilibrium. Numerical analysis reveals the possible theoretical improvements from using a preconditioner for all Mach number regimes. Numerical results confirm the expectations from the numerical analysis. Representative test cases include flows with previously troublesome embedded high-condition-number areas. Van Leer, Lee, and Roe recently developed an optimal, analytic preconditioning technique to reduce eigenvalue stiffness over the full Mach-number range. By multiplying the flux-balance residual with the preconditioning matrix, the acoustic wave speeds are scaled so that all waves propagate at the same rate, an essential property to eliminate inherent eigenvalue stiffness. This session discusses a synthesis of the thermochemical nonequilibrium flux-splitting developed by Grossman and Cinnella and the characteristic wave preconditioning of Van Leer into a powerful tool for implicitly solving two and three-dimensional flows with generalized finite-rate chemistry. For finite-rate chemistry, the state vector of unknowns is variable in length. Therefore, the preconditioning matrix extended to generalized finite-rate chemistry must accommodate a flexible system of moving waves. Fortunately, no new kind of wave appears in the system. The only existing waves are entropy and vorticity waves, which move with the fluid, and acoustic waves, which propagate in Mach number dependent directions. The nonequilibrium vibrational energies and species densities in the unknown state vector act strictly as convective waves. The essential concept for extending the preconditioning to generalized chemistry models is determining the differential variables which symmetrize the flux Jacobians. The extension is then straight-forward. This algorithm research effort will be released in a future version of the production level computational code coined the General Aerodynamic Simulation Program (GASP), developed by Walters, Slack, and McGrory.

Godfrey, Andrew G.↗

Investigating On-Orbit Attitude Determination Anomalies for the Solar Dynamics Observatory Mission

The Solar Dynamics Observatory (SDO) was launched on February 11, 2010 from Kennedy Space Center on an Atlas V launch vehicle into a geosynchronous transfer orbit. SDO carries a suite of three scientific instruments, whose observations are intended to promote a more complete understanding of the Sun and its effects on the Earth's environment. After a successful launch, separation, and initial Sun acquisition, the launch and flight operations teams dove into a commissioning campaign that included, among other things, checkout and calibration of the fine attitude sensors and checkout of the Kalman filter (KF) and the spacecraft s inertial pointing and science control modes. In addition, initial calibration of the science instruments was also accomplished. During that process of KF and controller checkout, several interesting observations were noticed and investigated. The SDO fine attitude sensors consist of one Adcole Digital Sun Sensor (DSS), two Galileo Avionica (GA) quaternion-output Star Trackers (STs), and three Kearfott Two-Axis Rate Assemblies (hereafter called inertial reference units, or IRUs). Initial checkout of the fine attitude sensors indicated that all sensors appeared to be functioning properly. Initial calibration maneuvers were planned and executed to update scale factors, drift rate biases, and alignments of the IRUs. After updating the IRU parameters, the KF was initialized and quickly reached convergence. Over the next few hours, it became apparent that there was an oscillation in the sensor residuals and the KF estimation of the IRU bias. A concentrated investigation ensued to determine the cause of the oscillations, their effect on mission requirements, and how to mitigate them. The ensuing analysis determined that the oscillations seen were, in fact, due to an oscillation in the IRU biases. The low frequencies of the oscillations passed through the KF, were well within the controller bandwidth, and therefore the spacecraft was actually following the oscillating biases, resulting in movement of the spacecraft on the order of plus or minus 20 arcsec. Though this level of error met the ACS attitude knowledge requirement of [35, 70, 70] arcsec, 3 sigma, the desire of the ACS and instrument teams was to remove as much of the oscillation as possible. The Kearfott IRUs have an internal temperature controller, designed to maintain the IRU temperature at a constant temperature of approximately 70 C, thus minimizing the change in the bias drift and scale factors of the mechanical gyros. During ground testing of the observatory, it was discovered that the 83-Hz control cycle of the IRU heaters put a tremendous amount of stress on the spacecraft battery. Analysis by the power systems team indicated that the constant charge/discharge on the battery due to the IRU thermal control cycle could potentially limit the life of the battery. After much analysis, the decision was made not to run the internal IRU heaters. Analysis of on orbit data revealed that the oscillations in the IRU bias had a connection to the temperature of the IRU; changes in IRU temperature resulted in changes in the amplitude and period of the IRU biases. Several mitigating solutions were investigated, the result of which was to tune the KF with larger IRU noise assumptions which allows the KF to follow and correct for the time-varying IRU biases.

Vess, Melissa F.↗

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗

Aeroelastic Analysis of a Distributed Electric Propulsion Wing

An aeroelastic analysis of a prototype distributed electric propulsion wing is presented. Results using MSC Nastran (Registered Trademark) doublet lattice aerodynamics are compared to those based on FUN3D Reynolds Averaged Navier- Stokes aerodynamics. Four levels of grid refinement were examined for the FUN3D solutions and solutions were seen to be well converged. It was found that no oscillatory instability existed, only that of divergence, which occurred in the first bending mode at a dynamic pressure of over three times the flutter clearance condition.

Massey, Steven J.↗

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations over the Indian Ocean

Observational analyses and a hierarchy of ocean general circulation model (OGCM) experiments were performed to understand the influence of oceanic processes on the warm sea surface temperature anomalies (SSTAs) prior to the convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden-Julian Oscillation (MJO), in the equatorial Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool region during the November-April season of the 2001-2012 period. 17/39 ISO events initiated in the Seychelles-Chagos Thermocline Ridge (SCTR) before propagating eastward; the remaining events initiated in the southern Arabian Sea (6) or Warm Pool (16) regions. The SCTR event set was notable in that it contained more global-scale MJOs (71-76%), as defined by the RMM and OMI indices, than the WP events (25-44%). Additionally, ~24% (44%) of the SCTR (Warm Pool) events were preceded by strong oceanic process-induced SSTAs of similar magnitude to those of shortwave radiative and turbulent heat fluxes. The Arabian Sea events, however, were not associated with statistically significant SSTA signals prior to convection. Based on a mixed layer heat budget analysis, entrainment and upwelling reduction were the dominant oceanic processes contributing to the warming, in contrast with boreal summer, when horizontal advection dominated. We examined several case studies, including primary MJO events, where oceanic Rossby waves were associated with the entrainment and upwelling reduction. Two simple atmospheric boundary layer convergence models revealed that the SSTAs contributed at least half of the total convergence and suggested that the ocean dynamical effect was responsible for the majority of SSTA-forced convergence for those case studies. These results underscore the need for climate prediction models to accurately represent the ocean structure and processes to include the effects of oceanic predictors.

West, B. Jason↗

A Combined Variable-Temperature Neutron Diffraction and Thermogravimetric Analysis Study on a Promising Oxygen Electrode, SrCo 0.9 Nb 0.1 O 3–δ , for Reversible Solid Oxide Fuel Cells

The present study reported herein investigates the temperature–structure–stoichiometry relationship of a promising oxygen electrode SrCo 0.9 Nb 0.1 O 3-δ over a temperature (T) range from room temperature (RT) to 900 °C. The techniques employed are variable-temperature neutron diffraction (VTND) and thermogravimetric analysis (TGA). At $T$ < 75 °C, VTND reveals a tetragonal (P4/mmm) structure with a G-type magnetic ordering. Above 75 °C, the nucleus structure remains the same, while the magnetic ordering disappears. A phase transition from tetragonal ($P4/mmm$) to cubic ($Pm\overline3m$) is observed at 412 °C, where the two Co sites and three O sites in the $P4/mmm$ phase converge to one equivalent site, respectively. The phase transition temperature coincides with the peak temperature of oxygen uptake obtained by TGA. It is also observed that the Nb dopant has no preferred Co site to occupy. The oxygen vacancies are mostly located at the O3 site surrounding the Co2 site in the $P4/mmm$ structure. The intermediate-spin state of Co 3+ at the Co2 site is responsible for the observed distortions of CoO 6 octahedra, i.e., elongation of Co2O 6 octahedra and shortening of Co1O 6 octahedra along the c-axis, which is a phenomenon known as Jahn–Teller distortion. At high temperatures, large thermal displacement factor for O 2– is observed with high concentration of oxygen vacancies, providing a structural environment favorable to high O 2– conductivity in Nb-doped SrCoO 3 -based oxygen electrode materials.

25 ENERGY STORAGE↗

Nanothermodynamics of large iron clusters by means of a flat histogram Monte Carlo method

The thermodynamics of iron clusters of various sizes, from 76 to 2452 atoms, typical of the catalyst particles used for carbon nanotubes growth, has been explored by a flat histogram Monte Carlo (MC) algorithm (called the σ-mapping), developed by Soudan et al. [J. Chem. Phys. 135, 144109 (2011), Paper I]. This method provides the classical density of states, g{sub p}(E{sub p}) in the configurational space, in terms of the potential energy of the system, with good and well controlled convergence properties, particularly in the melting phase transition zone which is of interest in this work. To describe the system, an iron potential has been implemented, called “corrected EAM” (cEAM), which approximates the MEAM potential of Lee et al. [Phys. Rev. B 64, 184102 (2001)] with an accuracy better than 3 meV/at, and a five times larger computational speed. The main simplification concerns the angular dependence of the potential, with a small impact on accuracy, while the screening coefficients S{sub ij} are exactly computed with a fast algorithm. With this potential, ergodic explorations of the clusters can be performed efficiently in a reasonable computing time, at least in the upper half of the solid zone and above. Problems of ergodicity exist in the lower half of the solid zone but routes to overcome them are discussed. The solid-liquid (melting) phase transition temperature T{sub m} is plotted in terms of the cluster atom number N{sub at}. The standard N{sub at}{sup −1/3} linear dependence (Pawlow law) is observed for N{sub at} >300, allowing an extrapolation up to the bulk metal at 1940 ±50 K. For N{sub at} <150, a strong divergence is observed compared to the Pawlow law. The melting transition, which begins at the surface, is stated by a Lindemann-Berry index and an atomic density analysis. Several new features are obtained for the thermodynamics of cEAM clusters, compared to the Rydberg pair potential clusters studied in Paper I.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmark solution of the dynamic response of a spherical shell at finite strain

Our paper describes the development of high fidelity solutions for the study of homogeneous (elastic and inelastic) spherical shells subject to dynamic loading and undergoing finite deformations. The goal of the activity is to provide high accuracy results that can be used as benchmark solutions for the verification of computational physics codes. Furthermore, the equilibrium equations for the geometrically non-linear problem are solved through mode expansion of the displacement field and the boundary conditions are enforced in a strong form. Time integration is performed through high-order implicit Runge–Kutta schemes. Finally, we evaluate accuracy and convergence of the proposed method by means of numerical examples with finite deformations and material non-linearities and inelasticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusion Synthetic Acceleration for High Order S{sub N} Transport on Meshes with Curved Surfaces

Thermal radiation transport is integral to multiphysics applications such as inertial confinement fusion and astrophysics. Materials in these problems can be exceptionally optically thick and diffusive. A transport solver must be capable of accurate and efficient solutions to such problems. Previously, Woods et al. established a proof-of-concept for high order (HO) discontinuous Galerkin finite element (DGFEM) transport on meshes with curved surfaces. They numerically demonstrated the spatial convergence is order p + 1. Remarkably, these convergence rates are relatively invariant to the amount of curvature of the spatial mesh. We present here an extension of this work that demonstrates a rapidly-convergent iterative solver for the FEM transport equations on meshes with curved sides. We use the modified interior penalty (MIP) diffusion synthetic acceleration (DSA) established by Wang and Ragusa. The partially consistent MIP equations stem from the symmetric interior penalty (IP) derivation of the diffusion equation. The IP method is not stable for optically thick cells so it was combined with the diffusion conforming form (DCF), a discretized diffusion equation derived from the discretized S N equations that is not stable for intermediate optical thicknesses. The MIP equations adapt a 'switch' to the IP method in optically thinner regions and the DCF method in optically thicker regions. The literature does not reveal any Fourier analyses performed for iterative transport on meshes with curved surfaces. Wang and Ragusa performed a Fourier analysis for orthogonal meshes and our results appear to be as expected. Previous literature has reported on some of the sensitivities that affect the spectral radius of DSA schemes: polynomial order, and scattering ratios. Presently, we focus on the sensitivity of the spectral radius to polynomial order, scattering ratio, the constant C, and mesh curvature for varying cell thicknesses. We also numerically study the behavior of the spectral radius to verify our DSA implementation. A sequence of problems is designed to exercise the acceleration. As the problems become more optically thick and diffusive, the accelerated solution continues to converge toward the analytic solution. We employ the Modular Finite Element Methods (MFEM) package, developed at Lawrence Livermore National Laboratory, to generate the linear algebraic system of equations for each angle in the discrete ordinates quadrature set, and use UMFPack to solve this system of equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis and application of minimum variance discrete linear system identification

An on-line minimum variance (MV) parameter identifier is developed which embodies both accuracy and computational efficiency. The formulation results in a linear estimation problem with both additive and multiplicative noise (AMN). The resulting filter which utilizes both the covariance of the parameter vector itself and the covariance of the error in identification is proven to be mean-square convergent and mean-square consistent. The MV parameter identification scheme is then used to construct a stable state and parameter estimation algorithm.

Kotob, S.↗

High-Fidelity Analysis of Six-Passenger Quadrotor Air Taxi Concept

A code-to-code comparison has been performed for high-fidelity simulations ofNASA’s six-passenger quadrotor air taxi con- cept vehicle. Themultidisciplinary simulations combine comprehensive rotorcraft dynamics with high-fidelity fluid dynamics obtained from an unsteady Navier-Stokes computational fluid dynamics code. The results are then compared with those in the literature. A quantitative comparison of converged trim angles has been performed and normal force, chord force, and pitching-moment coefficients arepresented for qualitativecomparison.

eVTOL↗

Failure statistics for commercial lithium ion batteries: A study of 24 pouch cells

There are relatively few publications that assess capacity decline in enough commercial cells to quantify cell-to-cell variation, but those that do show a surprisingly wide variability. Capacity curves cross each other often, a challenge for efforts to measure the state of health and predict the remaining useful life (RUL) of individual cells. We analyze capacity fade statistics for 24 commercial pouch cells, providing an estimate for the time to 5% failure. Our data indicate that RUL predictions based on remaining capacity or internal resistance are accurate only once the cells have already sorted themselves into “better” and “worse” ones. Analysis of our failure data, using maximum likelihood techniques, provide uniformly good fits for a variety of definitions of failure with normal and with 2- and 3-parameter Weibull probability density functions, but we argue against using a 3-parameter Weibull function for our data. pdf fitting parameters appear to converge after about 15 failures, although business objectives should ultimately determine whether data from a given number of batteries provides sufficient confidence to end lifecycle testing. Increased efforts to make batteries with more consistent lifetimes should lead to improvements in battery cost and safety.

25 ENERGY STORAGE↗

A TWO-MOMENT RADIATION HYDRODYNAMICS MODULE IN ATHENA USING A TIME-EXPLICIT GODUNOV METHOD

We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas and radiation variables including the non-stiff source terms, and a local implicit method to integrate the stiff source terms. We adopt the M{sub 1} closure relation and include all leading source terms to O({beta}{tau}). We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. Our code is dimensionally unsplit in one, two, and three space dimensions and is parallelized using MPI. The streaming and diffusion limits are well described by the M{sub 1} closure model, and our implementation shows excellent behavior for a problem with a concentrated radiation source containing both regimes simultaneously. Our operator-split method is ideally suited for problems with a slowly varying radiation field and dynamical gas flows, in which the effect of the RSLA is minimal. We present an analysis of the dispersion relation of RHD linear waves highlighting the conditions of applicability for the RSLA. To demonstrate the accuracy of our method, we utilize a suite of radiation and RHD tests covering a broad range of regimes, including RHD waves, shocks, and equilibria, which show second-order convergence in most cases. As an application, we investigate radiation-driven ejection of a dusty, optically thick shell in the ISM. Finally, we compare the timing of our method with other well-known iterative schemes for the RHD equations. Our code implementation, Hyperion, is suitable for a wide variety of astrophysical applications and will be made freely available on the Web.

79 ASTRONOMY AND ASTROPHYSICS↗

A low-mass hub–filament with double centre revealed in NGC 2071-North

We present the first analysis in NGC 2071-North as a resolved hub–filament featuring a double centre. This ∼1.5 × 1.5 pc 2 scale filament hub contains ∼500 M ⊙ . Seen from Planck, magnetic field lines may have facilitated the gathering of material at this isolated location. The energy balance analysis, supported by infalling gas signatures, reveals that these filaments are currently forming stars. Herschel 100 $\mu$m emission concentrates in the hub, at IRAS 05451+0037 and LkHα 316, and presents diffuse lobes and loops around them. We suggest that such a double centre could be formed, because the converging locations of filament pairs are offset, by 2.3 arcmin (0.27 pc). This distance also matches the diameter of a hub ring, seen in column density and molecular tracers, such as HCO + (1–0) and HCN (1–0), that may indicate a transition and the connection between the hub and the radiating filaments. We argue that all of the three components of the emission star LkHα 316 are in physical association. We find that a ∼0.06-pc-sized gas loop, attached to IRAS 05451+0037, can be seen at wavelengths all the way from the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS)-i to Herschel 100 $\mu$m. These observations suggest that both protostars at the double hub centre are interacting with the cloud material. In our 13 CO data, we do not seem to find the outflow of this region that was identified in the 1980s with much lower resolution.

79 ASTRONOMY AND ASTROPHYSICS↗