Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximation & 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Investigation of MMOD Impact on STS-115 Shuttle Payload Bay Door Radiator

The Orbiter radiator system consists of eight individual 4.6 m x 3.2 m panels located with four on each payload bay door. Forward panels #1 and #2 are 2.3 cm thick while the aft panels #3 and #4 have a smaller overall thickness of 1.3 cm. The honeycomb radiator panels consist of 0.028 cm thick Aluminum 2024-T81 facesheets and Al5056-H39 cores. The face-sheets are topped with 0.005 in. (0.127 mm) silver-Teflon tape. The radiators are located on the inside of the shuttle payload bay doors, which are closed during ascent and reentry, limiting damage to the on-orbit portion of the mission. Post-flight inspections at the Kennedy Space Center (KSC) following the STS-115 mission revealed a large micrometeoroid/orbital debris (MMOD) impact near the hinge line on the #4 starboard payload bay door radiator panel. The features of this impact make it the largest ever recorded on an orbiter payload bay door radiator. The general location of the damage site and the adjacent radiator panels can be seen in Figure 2. Initial measurements of the defect indicated that the hole in the facesheet was 0.108 in. (2.74 mm) in diameter. Figure 3 shows an image of the front side damage. Subsequent observations revealed exit damage on the rear facesheet. Impact damage features on the rear facesheet included a 0.03 in. diameter hole (0.76 mm), a approx.0.05 in. tall bulge (approx.1.3 mm), and a larger approx.0.2 in. tall bulge (approx.5.1 mm) that exhibited a crack over 0.27 in. (6.8 mm) long. A large approx.1 in. (25 mm) diameter region of the honeycomb core was also damaged. Refer to Figure 4 for an image of the backside damage to the panel. No damage was found on thermal blankets or payload bay door structure under the radiator panel. Figure 5 shows the front facesheet with the thermal tape removed. Ultrasound examination indicated a maximum facesheet debond extent of approximately 1 in. (25 mm) from the entry hole. X-ray examinations revealed damage to an estimated 31 honeycomb cells with an extent of 0.85 in. x 1.1 in. (21.6 x 27.9 mm). Pieces of the radiator at and surrounding the impact site were recovered during the repair procedures at KSC. They included the thermal tape, front facesheet, honeycomb core, and rear facesheet. These articles were examined at JSC using a scanning electron microscope (SEM) with an energy dispersive x-ray spectrometer (EDS). Figure 6 shows SEM images of the entry hole in the facesheet. The asymmetric height of the lip may be attributed to projectile shape and impact angle. Numerous instances of a glass-fiber organic matrix composite were observed in the facesheet tape sample. The fibers were approximately 10 micrometers in diameter and variable lengths. EDS analysis indicated a composition of Mg, Ca, Al, Si, and O. Figures 7 and 8 present images of the fiber bundles, which were believed to be circuit board material based on similarity in fiber diameter, orientation, consistency, and composition. A test program was initiated in an attempt to simulate the observed damage to the radiator facesheet and honeycomb. Twelve test shots were performed using projectiles cut from a 1.6 mm thick fiberglass circuit board substrate panel. Results from test HITF07017, shown in figures 9 and 10, correlates with the observed impact features reasonably well. The test was performed at 4.14 km/sec with an impact angle of 45 degrees using a cylindrical projectile with a diameter and length of 1.25 mm. The fiberglass circuit board material had a density of 1.65 g/cu cm, giving a projectile mass of 2.53 mg. An analysis was performed using the Bumper code to estimate the probability of impact to the shuttle from a 1.25 mm diameter particle. Table 1 shows a 1.6% chance (impact odds = 1 in 62) of a 1.25 mm or larger MMOD impact on the radiators of the vehicle during a typical ISS mission. There is a 0.4% chance (impact odds = 1 in 260) that a 1.25 mm or larger MMOD particle would impact the RCC wing leading edge and nose cap during a typical miion. Figure 11 illustrates the vulnerable areas of the wing leading edge reinforced carbon-carbon (RCC), an area of the vehicle that is very sensitive to impact damage. The highlighted red, orange, yellow, and light green areas would be expected to experience critical damage if impacted by an OD particle such as the one that hit the RH4 radiator panel on STS-115.

Hyde, J.↗
NASA NTRS1996

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗
NASA NTRS2001

The Large-Scale Plasmaspheric Density Trough Associated With the 24 May 2000 Geomagnetic Storm: IMAGE EUV Observations and Global Core Plasma Modeling

The IMAGE EUV imager observed a plasmaspheric density, trough in association with a geomagnetically active period on 24 May 2000. In EUV, this density, trough appears as an Archimedes spiral extending from Earth's shadow to approximately 1900 MLT. We present an analysis of this density trough using simulated EUV images, Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density, depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M. L.↗
DOE OSTI2022

Multiscale analysis in solids with unseparated scales: fine-scale recovery, error estimation, and coarse-scale adaptivity

There are several engineering applications in which the assumptions of homogenization and scale separation may be violated, in particular, for metallic structures constructed through additive manufacturing. Instead of resorting to direct numerical simulation of the macroscale system with an embedded fine scale, an alternative approach is to use an approximate macroscale constitutive model, but then estimate the model-form error using a posteriori error estimation techniques and subsequently adapt the macroscale model to reduce the error for a given boundary value problem and quantity of interest. Here, we investigate this approach to multiscale analysis in solids with unseparated scales using the example of an additively manufactured metallic structure consisting of a polycrystalline microstructure that is neither periodic nor statistically homogeneous. As a first step to the general nonlinear case, we focus here on linear elasticity in which each grain within the polycrystal is linear elastic but anisotropic.

42 ENGINEERING↗
NASA NTRS1988

Equilibrium gas flow computations. II - An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, equilibrium gas laws. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer flux-vector splittings, and Roe's approximate Riemann solver are presented for three-dimensional, time-varying grids. The approximations inherent in previous generalizations are discussed.

Vinokur, Marcel↗
NASA NTRS2000

IMAGE EUV Observations and Modeling of the Plasmaspheric Density Trough Associated with the 24 May 2000 Geomagnetic Storm

The IMAGE EUV imager observed a plasmaspheric density trough in association with a geomagnetically active period on 24 May 2000. At EUV wavelengths, this density trough appeared as an Archimedes spiral extending from Earth's shadow to approximately 1800 MLT. We present an analysis of this density trough using simulated EUV images. Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) [Gallagher et al., 20001 and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M.L.↗
NASA NTRS1985

Finite element engine blade structural optimization

The Structural Tailoring of Engine Blades (STAEBL) computer program was developed to perform engine fan blade numerical optimizations. These blade optimizations seek a minimum weight or cost design that satisfies realistic blade design constraints, by tuning one to twenty design variables. The STAEBL system has been generalized to include both fan and compressor blade numerical optimizations. The system analyses have been significantly improved through the inclusion of an efficient plate finite element analysis for blade stress and frequency determinations. Additionally, a finite element based approximate severe foreign object damage (FOD) analysis has been included. The new FOD analysis gives very accurate estimates of the full nonlinear bird ingestion solution. Optimizations of fan and compressor blades have been performed using the system, showing significant cost and weight reductions, while comparing very favorably with refined design validation procedures.

Brown, K. W.↗
NASA NTRS1994

Climatic, eustatic, and tectnoic controls on Quarternary deposits and landforms, Red Sea coast, Egypt

The degree to which local climatic variations, eustatic sea level fluctuations, and tectonic uplift have influenced the development of Quaternary marine and fluvial landforms and deposits along the Red Sea coast, Eastern Desert, was investigated using a combination of remote sensing and field data, age determinations of corals, and numerical simulations. False color composites generated from Landsat Thematic Mapper and SPOT image data, digital elevation models derived from sterophotogrammetric analysis of SPOT data, and field observations document that a approximately 10-km wide swath inland from the coast is covered in many places with coalescing alluvial fans of Quaternary age. Wadis cutting through the fans exhibit several pairs of fluvial terraces, and wadi walls expose alluvium interbedded with corraline limestone deposits Further, three distinct coral terraces are evident along the coatline. Climatic, eustatic, and tectonic uplift controls on the overall system were simulated using a cellular automata algorithm with the following characteristics: (1) uplift as a function of position and time, as defined by the elevations and ages of corals; (2) climatic variations driven by insolation changes associated with Milankovitch cycles; (3) sea level fluctuations based on U/Th ages of coral terraces and eustatic data; and (4) parametrized fluvial erosion and deposition. Results imply that the fans and coralline limestones were generated in a setting in which the tectonic uplift rate decreased over the Quarternary to negligible values at present. Coralline limestones formed furing eustatic highstands when alluvium was trapped uspstream and wadis filled with debris. During lowstands, wadis cut into sedimentary deposits; coupled with continuing uplift, fans were dissected, leaving remnant surfaces, and wadi-related terraces were generated by down cutting. Only landforms from the past three to four eustatic sea level cycles (i.e., approximately 300 to 400 kyr) are likely to have survived erosion and deposition associated with fluvial processes.

Arvidson, Raymond↗
NASA NTRS1982

Stability analysis of flexible wind turbine blades using finite element method

Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.

Kamoulakos, A.↗
NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗
NASA NTRSDate not supplied

Comparison of Multivariate Time Series Prediction Techniques for Emulating Noah-LSM Soil Moisture Outputs

Land surface models are crucial tools for many earth science applications including numerical weather prediction, water resource and crop monitoring, and climatological analysis. Given a set of atmospheric forcings, seasonal data, and static parameters, models like Noah-LSM solve for land surface quantities including skin temperature, sensible heat flux, and soil moisture. While these calculations are theoretically robust, they are often computationally expensive. Since artificial neural networks (ANNs) are universal function approximators, they can learn to emulate the output of a deterministic numerical model given a time series of input forcings, with the learned ANN having substantially shorter execution time. The ANN could efficiently parameterize other models, generate ensembles, and provide first-guess inputs for retrievals. As such, with the goal of developing a model that efficiently mimics the output of Noah-LSM given NLDAS2 forcings on a region covering much of the central US, we examine and compare several neural network architectures for the multi-horizon multivariate time series forecasting problem. Recent literature includes a diverse set of approaches including autoregressive architectures like LSTM and GRU, parametric and non-parametric statistical predictors (ForecastNet and MQRNN), self-attention (LSTM-attention-LSTM), and temporal convovlution (DeepTCN). We implement several of these models for the Noah-LSM prediction task, highlighting the features and challenges for each and providing practical insight on the training process.

Mitchell Dodson↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
NASA NTRS1973

The use of the modified Cholesky decomposition in divergence and classification calculations

The use of the Cholesky decomposition technique is analyzed as applied to the feature selection and classification algorithms used in the analysis of remote sensing data (e.g. as in LARSYS). This technique is approximately 30% faster in classification and a factor of 2-3 faster in divergence, as compared with LARSYS. Also numerical stability and accuracy are slightly improved. Other methods necessary to deal with numerical stablity problems are briefly discussed.

Vanroony, D. L.↗
NASA NTRS1973

The use of the modified Cholesky decomposition in divergence and classification calculations

This report analyzes the use of the modified Cholesky decomposition technique as applied to the feature selection and classification algorithms used in the analysis of remote sensing data (e.g., as in LARSYS). This technique is approximately 30% faster in classification and a factor of 2-3 faster in divergence, as compared with LARSYS. Also numerical stability and accuracy are slightly improved. Other methods necessary to deal with numerical stability problems are briefly discussed.

Van Rooy, D. L.↗
NASA NTRS1980

Solution of Euler's Equations of Motion and Eulerian Angles for near symmetric rigid bodies subject to constant moments

Analytic expressions are found for Euler's Equations of Motion and for the Eulerian Angles for both symmetric and near symmetric rigid bodies under the influence of arbitrary constant body-fixed torques. These solutions provide the body-fixed angular velocities and the attitude of the body, respectively, as functions of time. They are of special interest in applications to spinning spacecraft (such as the Galileo Spacecraft to be launched in 1984) because they include the effect of time-varying spin rate. Thus they can be applied to spin-up and spin-down maneuvers as well as to error analysis for thruster misalignments. The solutions are given for arbitrary initial conditions in terms of Fresnel, Sine and Cosine Integrals. Numerical integration of the governing differential equations has verified that the approximate analytic solutions are very accurate in many physical situations of interest.

Longuski, J. M.↗
NASA NTRS1984

Error analysis and correction of discrete solutions from finite element codes

Many structures are an assembly of individual shell components. Therefore, results for stresses and deflections from finite element solutions for each shell component should agree with the equations of shell theory. This paper examines the problem of applying shell theory to the error analysis and the correction of finite element results. The general approach to error analysis and correction is discussed first. Relaxation methods are suggested as one approach to correcting finite element results for all or parts of shell structures. Next, the problem of error analysis of plate structures is examined in more detail. The method of successive approximations is adapted to take discrete finite element solutions and to generate continuous approximate solutions for postbuckled plates. Preliminary numerical results are included.

Thurston, G. A.↗
NASA NTRS1987

Model-size reduction for the non-linear dynamic analysis of quasi-symmetric structures

A numerical technique is developed to reduce the size of models describing the nonlinear dynamic response of quasi-symmetric structures (i.e., structures with unsymmetric geometry). The response vectors of the structure are approximated by a linear combination of the symmetric and antisymmetric vectors at each time step. The mathematical formulation and numerical implementation of the method are described in detail, and results for a shallow laminated anisotropic panel of quadrilateral planform are presented in graphs and normalized contour plots.

Noor, Ahmed K.↗
NASA NTRS1956

Theoretical analysis of incompressible flow through a radial-inlet centrifugal impeller at various weight flows

A method for the solution of the incompressible nonviscous flow through a centrifugal impeller, including the inlet region, is presented. Several numerical solutions are obtained for four weight flows through an impeller at one operating speed. These solutions are refined in the leading-edge region. The results are presented in a series of figures showing streamlines and relative velocity contours. A comparison is made with the results obtained by using a rapid approximate method of analysis.

Kramer, James J↗