Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable coefficients”

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 91 records · Page 5

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗

Developing sensitivity coefficients for a state variable model of a flexible spacecraft with modular attitude control

The numerical value for the sensitivity coefficient of each scalar state variable with respect to each parameter of a five body model of the rotational dynamics of a flexible spacecraft with modular attitude control is generated. The maximum magnitudes of these sensitivity coefficients are arranged in descending order in one table while the final time (steady state) values of these sensitivity coefficients are arranged in order of descending magnitude in a second table.

Chichester, F. D.↗

Usefulness of Derived Frank Lead Parameters in Screening for Coronary Artery Disease and Cardiomyopathy

We investigated the accuracy of several known as well as newly-introduced derived Frank-lead ECG parameters in differentiating healthy individuals from patients with obstructive coronary artery disease (CAD) and cardiomyopathy (CM). Advanced high-fidelity 12-lead ECG tests (approx. 5-min supine) were first performed on a "training set" of 99 individuals: 33 with ischemic or dilated CM and low ejection fraction (EF less than 40%); 33 with catheterization-proven obstructive CAD but normal EF; and 33 age-/gender-matched healthy controls. The following derived Frank lead parameters were studied for their accuracy in detecting CAD and CM: the spatial ventricular gradient (VG), including its beat-to-beat coefficient of variability (VG CV); the spatial mean QRS (SM-QRS) and T-waves (SM-T) and their beat-to-beat coefficients of variability; the spatial ventricular activation time (VAT); the mean and maximum spatial QRS-T angles; and standard late potentials parameters (RMS40, fQRSD and LAS). Several of these parameters were accurate in discriminating between the control group and both diseased groups at p less than 0.0001. For example the fQRSD, VG CV, mean spatial QRS-T angle and VG minus SM-QRS (which is similar to the SM-T) had retrospective areas under the ROC curve of 0.78, 0.78, 0.80, and 0.84 (CAD vs. controls) and 0.93, 0.88, 0.98 and 0.99 (CM vs. controls), respectively. The single most effective parameter in discriminating between the CAD and CM groups was the spatial VAT (44 plus or minus 5.8 vs. 53 plus or minus 9.9 ms, p less than 0.0001), with an area under the ROC curve of 0.80. Since subsequent prospective analyses using new groups of patients and healthy subjects have yielded only slightly less accurate results, we conclude that derived Frank-lead parameters show great promise for potentially contributing to the development of a rapid and inexpensive resting ECG-based screening test for heart disease.

DePalma, J. L.↗

Quasi-steady flight to quasi-steady flight transition in a windshear - Trajectory optimization and guidance

The near-optimum guidance of an aircraft from quasi-steady flight to quasi-steady flight in a windshear is studied. The take-off problem is considered with reference to flight in a vertical plane; allowance is made for the presence of a downdraft as well as horizontal shear. It is assumed that the power setting is held at the maximum value and that the aircraft is controlled through the angle of attack. While the shear guidance and the initial aftershear guidance use constant gain coefficients, the final aftershear guidance employs a variable gain coefficient. The results show that the guidance scheme for quasi-steady flight recovery yields a transition from quasi-steady flight to quasi-steady flight which is close to that of the optimal trajectory; it guarantees the restoration of the initial quasi-steady state and has good stability properties.

Miele, A.↗

Further Progress in Noise Source Identification in High Speed Jets via Causality Principle

To locate noise sources in high-speed jets, the sound pressure fluctuations p/, measured at far field locations, were correlated with each of density p, axial velocity u, radial velocity v, puu and pvv fluctuations measured from various points in fully expanded, unheated plumes of Mach number 0.95, 1.4 and 1.8. The velocity and density fluctuations were measured simultaneously using a recently developed, non-intrusive, point measurement technique based on molecular Rayleigh scattering (Seasholtz, Panda, and Elam, AIAA Paper 2002-0827). The technique uses a continuous wave, narrow line-width laser, Fabry-Perot interferometer and photon counting electronics. The far field sound pressure fluctuations at 30 to the jet axis provided the highest correlation coefficients with all flow variables. The correlation coefficients decreased sharply with increased microphone polar angle, and beyond about 60 all correlation mostly fell below the experimental noise floor. Among all correlations < puu; p/> showed the highest values. Interestingly, , in all respects, were very similar to . The and correlations with 90deg microphone fell below the noise floor. By moving the laser probe at various locations in the jet it was found that the strongest noise source lies downstream of the end of the potential core and extends many diameters beyond. Correlation measurement from the lip shear layer showed a Mach number dependency. While significant correlations were measured in Mach 1.8 jet, values were mostly below the noise floor for subsonic Mach 0.95 jet. Various additional analyses showed that fluctuations from large organized structures mostly contributed to the measured correlation, while that from small scale structures fell below the noise floor.

Panda, J.↗

An improved multiple linear regression and data analysis computer program package

NEWRAP, an improved version of a previous multiple linear regression program called RAPIER, CREDUC, and CRSPLT, allows for a complete regression analysis including cross plots of the independent and dependent variables, correlation coefficients, regression coefficients, analysis of variance tables, t-statistics and their probability levels, rejection of independent variables, plots of residuals against the independent and dependent variables, and a canonical reduction of quadratic response functions useful in optimum seeking experimentation. A major improvement over RAPIER is that all regression calculations are done in double precision arithmetic.

Sidik, S. M.↗

Metal/Silicate Partitioning of P, Ga, and W at High Pressures and Temperatures: Dependence on Silicate Melt Composition

The distinctive pattern of element concentrations in the upper mantle provides essential evidence in our attempts to understand the accretion and differentiation of the Earth (e.g., Drake and Righter, 2002; Jones and Drake, 1986; Righter et al., 1997; Wanke 1981). Core formation is best investigated through use of metal/silicate partition coefficients for siderophile elements. The variables influencing partition coefficients are temperature, pressure, the major element compositions of the silicate and metal phases, and oxygen fugacity. Examples of studies investigating the effects of these variables on partitioning behavior are: composition of the metal phase by Capobianco et al. (1999) and Righter et al. (1997); silicate melt composition by Watson (1976), Walter and Thibault (1995), Hillgren et al. (1996), Jana and Walker (1997), and Jaeger and Drake (2000); and oxygen fugacity by Capobianco et al. (1999), and Walter and Thibault (1995). Here we address the relative influences of silicate melt composition, pressure and temperature.

Bailey, Edward↗

Spherical Harmonics Analysis of the ECMWF Global Wind Fields at the 10-Meter Height Level During 1985: A Collection of Figures Illustrating Results

Half-daily global wind speeds in the east-west (u) and north-south (v) directions at the 10-meter height level were obtained from the European Centre for Medium Range Weather Forecasts (ECMWF) data set of global analyses. The data set covered the period 1985 January to 1995 January. A spherical harmonic expansion to degree and order 50 was used to perform harmonic analysis of the east-west (u) and north-south (v) velocity field components. The resulting wind field is displayed, as well as the residual of the fit, at a particular time. The contribution of particular coefficients is shown. The time variability of the coefficients up to degree and order 3 is presented. Corresponding power spectrum plots are given. Time series analyses were applied also to the power associated with degrees 0-10; the results are included.

Sanchez, Braulio V.↗

Updated Chemical Kinetics and Sensitivity Analysis Code

An updated version of the General Chemical Kinetics and Sensitivity Analysis (LSENS) computer code has become available. A prior version of LSENS was described in "Program Helps to Determine Chemical-Reaction Mechanisms" (LEW-15758), NASA Tech Briefs, Vol. 19, No. 5 (May 1995), page 66. To recapitulate: LSENS solves complex, homogeneous, gas-phase, chemical-kinetics problems (e.g., combustion of fuels) that are represented by sets of many coupled, nonlinear, first-order ordinary differential equations. LSENS has been designed for flexibility, convenience, and computational efficiency. The present version of LSENS incorporates mathematical models for (1) a static system; (2) steady, one-dimensional inviscid flow; (3) reaction behind an incident shock wave, including boundary layer correction; (4) a perfectly stirred reactor; and (5) a perfectly stirred reactor followed by a plug-flow reactor. In addition, LSENS can compute equilibrium properties for the following assigned states: enthalpy and pressure, temperature and pressure, internal energy and volume, and temperature and volume. For static and one-dimensional-flow problems, including those behind an incident shock wave and following a perfectly stirred reactor calculation, LSENS can compute sensitivity coefficients of dependent variables and their derivatives, with respect to the initial values of dependent variables and/or the rate-coefficient parameters of the chemical reactions.

Radhakrishnan, Krishnan↗

Aerodynamic Effects Caused by Icing of an Unswept NACA 65A004 Airfoil

The effects of ice formations on the section lift, drag, and pitching-moment coefficients of an unswept NACA 65A004 airfoil section of 6-foot chord were studied.. The magnitude of the aerodynamic penalties was primarily a function of the shape and size of the ice formation near the leading edge of the airfoil. The exact size and shape of the ice formations were determined photographically and found to be complex functions of the operating and icing conditions. In general, icing of the airfoil at angles of attack less than 40 caused large increases in section drag coefficients (as much as 350 percent in 8 minutes of heavy glaze icing), reductions in section lift coefficients (up to 13 percent), and changes in the pitching-moment coefficient from diving toward climbing moments. At angles of attack greater than 40 the aerodynamic characteristics depended mainly on the ice type. The section drag coefficients generally were reduced by the addition of rime ice (by as much as 45 percent in 8 minutes of icing). In glaze icing, however, the drag increased at these angles of attack. The section lift coefficients were variably affected by rime-ice formations; however, in glaze icing, lift increases at high angles of attack amounted to as much as 9 percent for an icing time of 8 minutes. Pitching-moment-coefficient changes in icing conditions were somewhat erratic and depended on the icing condition. Rotation of the iced airfoil to angles of attack other than that at which icing occurred caused sufficiently large changes in the pitching-moment coefficient that, in flight, rapid corrections in trim might be required in order to avoid a hazardous situation.

Gray, Vernon H.↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LSDYNA (Livermore Software Technology Corporation), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic yield function with a nonassociative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Polymer Matrix Composites↗

Theoretical Development of an Orthotropic Elasto-Plastic Generalized Composite Material Model

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within LS-DYNA (Registered), there are several features that have been identified that could improve the predictive capability of a composite model. To address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed and is being implemented into LS-DYNA as MAT_213. A key feature of the improved material model is the use of tabulated stress-strain data in a variety of coordinate directions to fully define the stress-strain response of the material. To date, the model development efforts have focused on creating the plasticity portion of the model. The Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic material model with a non-associative flow rule. The coefficients of the yield function, and the stresses to be used in both the yield function and the flow rule, are computed based on the input stress-strain curves using the effective plastic strain as the tracking variable. The coefficients in the flow rule are computed based on the obtained stress-strain data. The developed material model is suitable for implementation within LS-DYNA for use in analyzing the nonlinear response of polymer composites.

Plasticity↗

Effect of variable viscosity and thermal conductivity on high-speed slip flow between concentric cylinders

Schamberg was the first to solve the differential equations of slip flow, including the Burnett terms, for concentric circular cylinders assuming constant coefficients of viscosity and thermal conductivity. The problem is solved for variable coefficients of viscosity and thermal conductivity in this paper by applying a transformation which leads to an iteration method. Starting with the solution for constant coefficients, this method enables one to approximate the solution for variable coefficients very closely after one or two steps. Satisfactory results are shown to follow from Schamberg's solution by using his values of constant coefficients multiplied by a constant factor, leading to what are denoted as the effective coefficients of viscosity and thermal conductivity.

Lin, T C↗

Future Climate Change Significantly Alters Interannual Wheat Yield Variability Over Half of Harvested Areas

Climate change affects the spatial and temporal distribution of crop yields, which can critically impair food security across scales. A number of previous studies have assessed the impact of climate change on mean crop yield and future food availability, but much less is known about potential future changes in interannual yield variability. Here, we evaluate future changes in relative interannual global wheat yield variability (the coefficient of variation (CV)) at 0.25° spatial resolution for two representative concentration pathways (RCP4.5 and RCP8.5). A multi-model ensemble of crop model emulators based on global process-based models is used to evaluate responses to changes in temperature, precipitation, and CO2. The results indicate that over 60% of harvested areas could experience significant changes in interannual yield variability under a high-emission scenario by the end of the 21st century (2066–2095). About 31% and 44% of harvested areas are projected to undergo significant reductions of relative yield variability under RCP4.5 and RCP8.5, respectively. In turn, wheat yield is projected to become more unstable across 23% (RCP4.5) and 18% (RCP8.5) of global harvested areas—mostly in hot or low fertilizer input regions, including some of the major breadbasket countries. The major driver of increasing yield CV change is the increase in yield standard deviation, whereas declining yield CV is mostly caused by stronger increases in mean yield than in the standard deviation. Changes in temperature are the dominant cause of change in wheat yield CVs, having a greater influence than changes in precipitation in 53% and 72% of global harvested areas by the end of the century under RCP4.5 and RCP8.5, respectively. This research highlights the potential challenges posed by increased yield variability and the need for tailored regional adaptation strategies.

Climate change↗

Rotorcraft derivative identification from analytical models and flight test data

A general procedure is presented for systematic development of rotorcraft models for use in systems identification, which includes fuselage and rotor degrees of freedom (DOF). Formulations for rigid blade flap and lag as well as the normal mode representation of an elastic blade are developed for hingeless and articulated rotor systems. The method of multiblade coordinates is used to obtain linear constant coefficient state variable models of various levels of approximation. Two of the approximate models, a 6 DOF, are identified from a nonlinear articulated helicopter computer simulation. The results demonstrate the accuracy attainable for each model. Advanced results outline the status of rotorcraft modeling and systems identification and indicate areas that require further investigation.

Molusis, J. A.↗