SEARCH · Search NASA

Results for “Complete linearization method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Analytic Nuclear Gradients for Complete Active Space Linearized Pair-Density Functional Theory

Accurately modeling photochemical reactions is difficult due to the presence of conical intersections and locally avoided crossings, as well as the inherently multiconfigurational character of excited states. As such, one needs a multistate method that incorporates state interaction in order to accurately model the potential energy surface at all nuclear coordinates. The recently developed linearized pair-density functional theory (L-PDFT) is a multistate extension of multiconfiguration PDFT, and it has been shown to be a cost-effective post-MCSCF method (as compared to more traditional and expensive multireference many-body perturbation methods or multireference configuration interaction methods) that can accurately model potential energy surfaces in regions of strong nuclear–electronic coupling in addition to accurately predicting Franck–Condon vertical excitations. Here, in this paper, we report the derivation of analytic gradients for L-PDFT and their implementation in the PySCF-forge software, and we illustrate the utility of these gradients for predicting ground- and excited-state equilibrium geometries and adiabatic excitation energies for formaldehyde, s-trans-butadiene, phenol, and cytosine.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part I—Analysis Methods

Use of corrosion-resistant cladding can greatly extend the design life of structural components in many advanced reactor systems. However, there are currently no ASME design rules for cladded components to guard against creep-fatigue failure and ratcheting strain accumulation in elevated temperature nuclear service. This paper, presented in two parts, addresses this gap by proposing a design strategy for cladded components that do not require long-term testing of clad materials. The proposed approach relies on approximate design analysis methods for two types of clad materials—a soft clad that creeps much faster than and has lower yield stress than the class A base material and a hard clad that creeps much slower than and has higher yield stress than the class A base material. Part I discusses the approximate analysis strategies for the clad materials—treat a soft clad as perfectly compliant and a hard clad as linear elastic—and Part II develops a complete set of design rules for each of the two types of cladded components. In conclusion, finite element analyses of representative high temperature reactor components show that the proposed design analysis methods can bound the design quantities in soft cladded components and approximate the design quantities in hard cladded components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1976

Some inconsistencies of the finite element method as applied to inelastic response

The inadequacy of a two noded beam-column element with a linear axial and a cubic transverse displacement field for inelastic analysis is demonstrated. For complete equilibrium satisfaction in the linear elastic range a three noded beam-column element is shown to be consistent. Next, the sensitivity of the inelastic response to numerical solutions of the inelastic response of a cantilever beam resulting from approximate integration of strain energy are brought out and finally, consequences of this on the nonlinear transient response of structures are considered.

Kamat, M. P.↗

NASA NTRS1974

Calculation of unsteady transonic aerodynamics for oscillating wings with thickness (computer program)

A computer program has been developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program of Rodemich and Andrew, and accounts for the nonuniform flow caused by finite thickness by application of the local linearization concept. A brief description of each subroutine is given, and the method of input is shown in detail. A sample problem as well as a complete listing of the computer program are presented.

Ruo, S. Y.↗

NASA NTRS1994

A grid generation and flow solution method for the Euler equations on unstructured grids

A grid generation and flow solution algorithm for the Euler equations on unstructured grids is presented. The grid generation scheme utilizes Delaunay triangulation and self-generates the field points for the mesh based on cell aspect ratios and allows for clustering near solid surfaces. The flow solution method is an implicit algorithm in which the linear set of equations arising at each time step is solved using a Gauss Seidel procedure which is completely vectorizable. In addition, a study is conducted to examine the number of subiterations required for good convergence of the overall algorithm. Grid generation results are shown in two dimensions for a National Advisory Committee for Aeronautics (NACA) 0012 airfoil as well as a two-element configuration. Flow solution results are shown for two-dimensional flow over the NACA 0012 airfoil and for a two-element configuration in which the solution has been obtained through an adaptation procedure and compared to an exact solution. Preliminary three-dimensional results are also shown in which subsonic flow over a business jet is computed.

Anderson, W. Kyle↗

DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

NASA NTRS1976

Simplified solution for point contact deformation between two elastic solids

A linear-regression by the method of least squares is made on the geometric variables that occur in the equation for point contact deformation. The ellipticity and the complete eliptic integrals of the first and second kind are expressed as a function of the x, y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). These simplified equations enable one to calculate easily the point-contact deformation to within 3 percent without resorting to charts or numerical methods.

Brewe, D. E.↗

NASA NTRS1980

Free-stream disturbances, continuous eigenfunctions, boundary-layer instability and transition

A rational foundation is provided for the application of the linear stability theory of parallel shear flows to transition prediction. An explicit method is given for carrying out the necessary calculations. The expansions are shown to be complete. Sample calculations show that a typical boundary layer is very sensitive to vorticity disturbance in the inner boundary layer near the critical layer. Vorticity disturbances three or four boundary layer thicknesses above the boundary are nearly uncoupled from the boundary layer, in that the amplitudes of the discrete Tollmein-Schlichting waves are an extremely small fraction of the amplitude of the disturbance.

Salwen, H.↗

NASA NTRS1992

Grid generation and flow solution method for Euler equations on unstructured grids

A grid generation and flow solution algorithm for the Euler equations on unstructured grids is presented. The grid generation scheme, which uses Delaunay triangulation, generates the field points for the mesh based on cell aspect ratios and allows clustering of grid points near solid surfaces. The flow solution method is an implicit algorithm in which the linear set of equations arising at each time step is solved using a Gauss-Seidel procedure that is completely vectorizable. Also, a study is conducted to examine the number of subiterations required for good convergence of the overall algorithm. Grid generation results are shown in two dimensions for an NACA 0012 airfoil as well as a two element configuration. Flow solution results are shown for a two dimensional flow over the NACA 0012 airfoil and for a two element configuration in which the solution was obtained through an adaptation procedure and compared with an exact solution. Preliminary three dimensional results also are shown in which the subsonic flow over a business jet is computed.

Anderson, W. Kyle↗

NASA NTRS2018

Realistic Covariance Generation for the GPM Spacecraft

A covariance realism process for NASA's Global Precipitation Measurement (GPM) spacecraft is detailed. The GPM spacecraft is in a low earth orbit, and performs collision avoidance maneuvers few times a year. Currently GPM is below the International Space Station (ISS). So, in addition to cataloged debris objects, GPM must contend with smallsat/cubesat objects that are deployed from the ISS. Both operational scenarios require complete knowledge of the expected GPM prediction errors as a function of time. In this study, we present a method for generating realistic predicted covariance that uses linear propagation of the covariance with the addition of process noise. Further analyses are presented for the process noise ''tuning'' that generates an inflation factor based on the observed error statistics of the predictive satellite trajectories when compared to the definitive ones. Different tuning strategies are considered and compared via a Goodness-of-Fit testing for the Gaussian properties of the scaled covariance. SpaceNav's realistic covariance generation approach takes into account the contribution of predicted maneuver errors in the increased propagation uncertainty. Corresponding maneuver uncertainty is injected into the state uncertainty, and is used within the collision avoidance process to determine the collision risk for close approach events that follow a maneuver. This is a critical step in the maneuver planning process that provides the satellite operator with an accurate quantification of the collision probability for planned maneuvers. Using this information, an informed decision can be made to proceed with a maneuver if the collision risk is acceptable. This approach is validated by Monte-Carlo simulations and results are presented.

spacecraft uncertainty propagation↗

NASA NTRS2019

Electromagnetic Scattering by Discrete Random Media. IV: Coherent Backscattering

The problem of backscattering of light by a discrete random medium illuminated by an obliquely incident plane electromagnetic wave is considered.The analysis is performed in a linear-polarization basis and includes a complete derivation of the cross reflection matrix for a layer with densely and sparsely distributed particles, the design of an approximate method for computing the ladder and cross reflection matrices in the case of a semi-infinite medium with a sparse distribution of particles, the derivation of the relations between the elements of the ladder and cross reflection matrices in the exact backscattering direction for dense and sparse media, and the development of practical algorithms for solving the underlying integral equations by the method of Picard iterations and the discrete ordinate method. Simulation results for particles with large size parameters are also presented.

Adrian Doicu↗

NASA NTRS1967

Computer program provides linear sampled- data analysis for high order systems

Computer program performs transformations in the order S-to W-to Z to allow arithmetic to be completed in the W-plane. The method is based on a direct transformation from the S-plane to the W-plane. The W-plane poles and zeros are transformed into Z-plane poles and zeros using the bilinear transformation algorithm.

Bunn, D. B.↗

NASA NTRS1975

Initialization procedures for primitive equation models

A linear analysis and comparison of the damping properties of six dynamic initialization schemes is presented, indicating that the Okamura-Rivas scheme has the most efficient damping properties over the whole frequency range, and suggesting that it should be faster than the other methods and given more stable results. The results obtained with a nonlinear shallow water equations model agree well with the linear analysis. The Okamura-Rivas scheme attains complete balance in the equivalent of 5 to 6 hours of leapfrog forecasting, and requires in this model an order of magnitude less computation than the balance equation solution.

Grant, W. K. F.↗

NASA NTRS1977

Simplified solution for elliptical-contact deformation between two elastic solids

A linear regression by the method of least squares is made on the geometric variables that occur in the equation for elliptical-contact deformation. The ellipticity and the complete elliptic integrals of the first and second kind are expressed as a function of the x,y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). The procedure for solving these variables without the use of charts or a high-speed computer would be quite tedious. These simplified equations enable one to calculate easily the elliptical-contact deformation to within 3 percent accuracy without resorting to charts or numerical methods.

Brewe, D. E.↗

NASA NTRS1992

Effect of surface tension on the onset of convection in a double-diffusive layer

The effect of surface tension on the stability of a double-diffusive layer is considered using linear stability analysis. The surface tension is assumed to vary linearly with temperature and solute concentration. The eigenvalue problem is solved by the Galerkin method. Results show that the predicted stability boundary based on Marangoni effects alone is completely altered in the presence of buoyancy effects induced by low gravity levels (about 10 exp -5 g). At reduced gravity levels, salt-finger instability may onset in the overstable mode due to the stabilizing effect of surface tension. Fluid properties in terms of the Prandtl and the Lewis numbers have a profound effect on the stability conditions; opposite stability characteristics are found in salt solutions and in molten metals.

Chen, C. F.↗

NASA NTRS2008

Design and Field Test of a Mass Efficient Crane for Lunar Payload Handling and Inspection: The Lunar Surface Manipulation System

Devices for lifting, translating and precisely placing payloads are critical for efficient Earthbased construction operations. Both recent and past studies have demonstrated that devices with similar functionality will be needed to support lunar outpost operations. Lunar payloads include: a) prepackaged hardware and supplies which must be unloaded from landers and then accurately located at their operational site, b) sensor packages used for periodic inspection of landers, habitat surfaces, etc., and c) local materials such as regolith which require grading, excavation and placement. Although several designs have been developed for Earth based applications, these devices lack unique design characteristics necessary for transport to and use on the harsh lunar surface. These design characteristics include: a) composite components, b) compact packaging for launch, c) simple in-field reconfiguration and repair, and d) support for tele-operated or automated operations. Also, in contrast to Earth-based construction, where special purpose devices dominate a construction site, a lunar outpost will require versatile devices which provide operational benefit from initial construction through sustained operations. This paper will detail the design of a unique, high performance, versatile lifting device designed for operations on the lunar surface. The device is called the Lunar Surface Manipulation System to highlight the versatile nature of the device which supports conventional cable suspended crane operations as well as operations usually associated with a manipulator such as precise positioning where the payload is rigidly grappled by a tool attached to the tip of the device. A first generation test-bed to verify design methods and operational procedures is under development at the NASA Langley Research Center and recently completed field tests at Moses Lake Washington. The design relied on non-linear finite element analysis which is shown to correlate favorably with laboratory experiments. A key design objective, reviewed in this paper, is the device s simplicity, resulting from a focus on the minimum set of functions necessary to perform payload offload. Further development of the device has the potential for significant mass savings, with a high performance device incorporating composite elements estimated to have a mass less than 3% of the mass of the maximum lunar payload lifted at the tip. The paper will conclude with future plans for expanding the operational versatility of the device.

Doggett, William R.↗

NASA NTRS1989

An interpretation and solution of ill-conditioned linear equations

Data insufficiency, poorly conditioned matrices and singularities in equations occur regularly in complex optimization, correlation, and interdisciplinary model studies. This work concerns itself with two methods of obtaining certain physically realistic solutions to ill-conditioned or singular algebraic systems of linear equations arising from such studies. Two efficient computational solution procedures that generally lead to locally unique solutions are presented when there is insufficient data to completely define the model, or a least-squares error formulation of this system results in an ill-conditioned system of equations. If it is assumed that a reasonable estimate of the uncertain data is available in both cases cited above, then we shall show how to obtain realistic solutions efficiently, in spite of the insufficiency of independent data. The proposed methods of solution are more efficient than singular-value decomposition for dealing with such systems, since they do not require solutions for all the non-zero eigenvalues of the coefficient matrix.

Ojalvo, I. U.↗

DOE OSTI2023

Reconstructing cosmological initial conditions from late-time structure with convolutional neural networks

ABSTRACT We present a method to reconstruct the initial linear-regime matter density field from the late-time non-linearly evolved density field in which we channel the output of standard first-order reconstruction to a convolutional neural network (CNN). Our method shows dramatic improvement over the reconstruction of either component alone. We show why CNNs are not well-suited for reconstructing the initial density directly from the late-time density: CNNs are local models, but the relationship between initial and late-time density is not local. Our method leverages standard reconstruction as a preprocessing step, which inverts bulk gravitational flows sourced over very large scales, transforming the residual reconstruction problem from long-range to local and making it ideally suited for a CNN. We develop additional techniques to account for redshift distortions, which warp the density fields measured by galaxy surveys. Our method improves the range of scales of high-fidelity reconstruction by a factor of 2 in wavenumber above standard reconstruction, corresponding to a factor of 8 increase in the number of well-reconstructed modes. In addition, our method almost completely eliminates the anisotropy caused by redshift distortions. As galaxy surveys continue to map the Universe in increasingly greater detail, our results demonstrate the opportunity offered by CNNs to untangle the non-linear clustering at intermediate scales more accurately than ever before.

Astronomy & Astrophysics↗