SEARCH · Search NASA

Results for “variational 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 1,747 records · Page 97

All sources NASA NTRS DOE OSTI

NASA NTRS1960

The Numerical Calculation of Flow Past Conical Bodies Supporting Elliptic Conical Shock Waves at Finite Angles of Incidence

The inverse method, with the shock wave prescribed to be an elliptic cone at a finite angle of incidence, is applied to calculate numerically the supersonic perfect-gas flow past conical bodies not having axial symmetry. Two formulations of the problem are employed, one using a pair of stream functions and the other involving entropy and components of velocity. A number of solutions are presented, illustrating the numerical methods employed, and showing the effects of moderate variation of the initial parameters.

Briggs, Benjamin R.↗

NASA NTRS1995

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

NASA NTRS2001

Sensitivity Studies for Space-Based Global Measurements of Atmospheric Carbon Dioxide

Carbon dioxide (CO2) is well known as the primary forcing agent of global warming. Although the climate forcing due to CO2 is well known, the sources and sinks of CO2 are not well understood. Currently the lack of global atmospheric CO2 observations limits our ability to diagnose the global carbon budget (e.g., finding the so-called "missing sink") and thus limits our ability to understand past climate change and predict future climate response. Space-based techniques are being developed to make high-resolution and high-precision global column CO2 measurements. One of the proposed techniques utilizes the passive remote sensing of Earth's reflected solar radiation at the weaker vibration-rotation band of CO2 in the near infrared (approx. 1.57 micron). We use a line-by-line radiative transfer model to explore the potential of this method. Results of sensitivity studies for CO2 concentration variation and geophysical conditions (i.e., atmospheric temperature, surface reflectivity, solar zenith angle, aerosol, and cirrus cloud) will be presented. We will also present sensitivity results for an O2 A-band (approx. 0.76 micron) sensor that will be needed along with CO2 to make surface pressure and cloud height measurements.

Mao, Jian-Ping↗

NASA NTRS1999

Fundamental Study of Material Flow in Friction Stir Welds

The presented research project consists of two major parts. First, the material flow in solid-state, friction stir, butt-welds as been investigated using a marker insert technique. Changes in material flow due to welding parameter as well as tool geometry variations have been examined for different materials. The method provides a semi-quantitative, three-dimensional view of the material transport in the welded zone. Second, a FSW process model has been developed. The fully coupled model is based on fluid mechanics; the solid-state material transport during welding is treated as a laminar, viscous flow of a non-Newtonian fluid past a rotating circular cylinder. The heat necessary for the material softening is generated by deformation of the material. As a first step, a two-dimensional model, which contains only the pin of the FSW tool, has been created to test the suitability of the modeling approach and to perform parametric studies of the boundary conditions. The material flow visualization experiments agree very well with the predicted flow field. Accordingly, material within the pin diameter is transported only in the rotation direction around the pin. Due to the simplifying assumptions inherent in the 2-D model, other experimental data such as forces on the pin, torque, and weld energy cannot be directly used for validation. However, the 2-D model predicts the same trends as shown in the experiments. The model also predicts a deviation from the "normal" material flow at certain combinations of welding parameters, suggesting a possible mechanism for the occurrence of some typical FSW defects. The next step has been the development of a three-dimensional process model. The simplified FSW tool has been designed as a flat shoulder rotating on the top of the workpiece and a rotating, cylindrical pin, which extends throughout the total height of the flow domain. The thermal boundary conditions at the tool and at the contact area to the backing plate have been varied to fit experimental data such as temperature profiles, torque and tool forces. General aspects of the experimentally visualized material flow pattern are confirmed by the 3-D model.

Reynolds, Anthony P.↗

NASA NTRS2002

Effects of Liner Geometry on Acoustic Impedance

Current aircraft engine nacelles typically contain acoustic liners consisting of perforated sheets bonded onto honeycomb cavities. Numerous models have been developed to predict the acoustic impedance of these liners in the presence of grazing flow, and to use that information with aeroacoustic propagation codes to assess nacelle liner noise suppression. Recent efforts have provided advances in impedance education methodologies that offer more accurate determinations of acoustic liner properties in the presence of grazing flow. The current report provides the results of a parametric study, in which a finite element method was used to assess the effects of variations of the following geometric parameters on liner impedance, with and without the presence of grazing flow: percent open area, sheet thickness, sheet thickness-to-hole diameter ratio and cavity depth. Normal incidence acoustic impedances were determined for eight acoustic liners, consisting of punched aluminum facesheets bonded to hexcell honeycomb cavities. Similar liners were tested in the NASA Langley Research Center grazing incidence tube to determine their response in the presence of grazing flow. The resultant data provide a quantitative assessment of the effects of these perforate, single-layer liner parameters on the acoustic impedance of the liner.

Jones, Michael G.↗

NASA NTRS2003

Meter-Scale 3-D Models of the Martian Surface from Combining MOC and MOLA Data

We have extended our previous efforts to derive through controlled photoclinometry, accurate, calibrated, high-resolution topographic models of the martian surface. The process involves combining MGS MOLA topographic profiles and MGS MOC Narrow Angle images. The earlier work utilized, along with a particular MOC NA image, the MOLA topographic profile that was acquired simultaneously, in order to derive photometric and scattering properties of the surface and atmosphere so as to force the low spatial frequencies of a one-dimensional MOC photoclinometric model to match the MOLA profile. Both that work and the new results reported here depend heavily on successful efforts to: 1) refine the radiometric calibration of MOC NA; 2) register the MOC to MOLA coordinate systems and refine the pointing; and 3) provide the ability to project into a common coordinate system, simultaneously acquired MOC and MOLA with a single set of SPICE kernels utilizing the USGS ISIS cartographic image processing tools. The approach described in this paper extends the MOC-MOLA integration and cross-calibration procedures from one-dimensional profiles to full two-dimensional photoclinometry and image simulations. Included are methods to account for low-frequency albedo variations within the scene.

Soderblom, Laurence A.↗

NASA NTRS2012

Investigation of Effects of Material Architecture on the Elastic Response of a Woven Ceramic Matrix Composite

To develop methods for quantifying the effects of the microstructural variations of woven ceramic matrix composites on the effective properties and response of the material, a research program has been undertaken which is described in this paper. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, CVI SiC/SiC, composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents and collect relevant statistics such as within ply tow spacing. This information was then used to build two dimensional finite element models that approximated the observed section geometry. With the aid of geometrical models generated by the microstructural characterization process, finite element models were generated and analyses were performed to quantify the effects of the microstructure and its variation on the effective stiffness and areas of stress concentration of the material. The results indicated that the geometry and distribution of the porosity appear to have significant effects on the through-thickness modulus. Similarly, stress concentrations on the outer surface of the composite appear to correlate to regions where the transverse tows are separated by a critical amount.

Goldberg, Robert K.↗

NASA NTRS2015

Effect of Test Specimen Shape and Size on Interlaminar Tensile Properties of Advanced Carbon-Carbon Composites

The interlaminar tensile strength of 1000-tow T-300 fiber ACC-6 carbon-carbon composites was measured using the method of bonding the coupons to adherends at room temperature. The size, 0.70 to 1.963 inches maximum width or radius, and shape, round or square, of the test coupons were varied to determine if the test method was sensitive to these variables. Sixteen total variations were investigated and the results modeled.

Vaughn, Wallace L.↗

NASA NTRSDate not supplied

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗

NASA NTRS2018

An Algorithm for Trajectory Generation in Redundant Manipulators with Joint Transmission Accommodation

Trajectory generation for manipulators involves generating incremental updates of joint variables to achieve a desired end–eﬀector motion. The Jacobian matrix maps incremental joint motion to incremental end–eﬀector motion in a linear fashion and is typically used in manipulator trajectory generation algorithms. In the case of a redundant manipulator, the Jacobian matrix is not square or invertible and therefore algorithms based on pseudoinverses and their variations are commonly used for trajectory generation. These methods either are computationally not eﬃcient or do not utilize all the joints in motion generation and therefore do not completely exploit the redundancy of the manipulator. The method presented in this paper is a simple method that maximizes transmission of all joint variables onto a desired end–eﬀector motion trajectory. The method is based on aligning the null–space of an augmented Jacobian matrix with the path of the desired end–eﬀector motion, from which a linear combination of joints that projects fully onto the desired end–eﬀector trajectory is obtained. In this manner, all joints of the redundant manipulator are used to generate the end–eﬀector trajectory accommodating the ability of each joint in terms of its motion transmission.

Ravani, Bahram↗

DOE OSTI2018

On Concavity and Utilization Analytics of Wide-Area Network Transport Protocols

Wide-area data transfers in high-performance computing infrastructures and big data applications are increasingly carried over dynamically provisioned dedicated network connections, which provide high capacities with no competing traffic. Throughput of TCP and UDT over wide-area connections depends, often in a complex non-linear way, on the congestion control mechanism, buffer sizes, and the number of parallel streams. In addition, our extensive TCP and UDT throughput measurements and time traces over a suite of physical and emulated 10 Gbps connections with 0-366 ms round-trip times (RTTs) show significant statistical and temporal variations. Consequently, parameter selection and optimization for transport methods require data analytics to complement analytical throughput models. We present analytics based on the concavity-convexity of throughput profiles which provide insights into peak throughput and superior/inferior trend compared to linear interpolations based on RTT. In particular, we propose the utilization-concavity coefficient, a scalar metric that incorporates both utilization and concavity profile to characterize the overall performance of a transport protocol. These measurements-based anaytics enable us to select a transport protocol and its parameters for a given host and connection configuration to achieve high throughput with statistical guarantees.

Rao, Nageswara↗

DOE OSTI2022

Using distributed power electronics-based devices to improve the voltage and frequency stability of distribution systems

Methods include, in response to a line frequency variation of a power grid, adjusting a voltage setpoint of a voltage regulator coupled to the power grid at a grid edge to maintain a voltage at the grid edge, wherein the adjusting the regulated voltage setpoint is configured to reduce the line frequency variation to stabilize the line frequency of the power grid. Apparatus include a voltage regulator configured to couple to a power grid at a grid edge and to maintain a voltage at the grid edge, wherein the voltage regulator is further configured to adjust a voltage setpoint of a voltage regulator in response to a line frequency variation of the power grid to reduce the line frequency variation and stabilize the line frequency of the power grid.

Du, Wei↗

DOE OSTI2016

fuel cycle cost analysis of Th/U MOX fuel in a block-type HTR - 18619

With a rapid development of nuclear energy in these years, the supply of nuclear fuel is becoming more and more important. Thorium regains attention because of its abundant reserves and excellent physical properties. As a next generation reactor, the high temperature reactor (HTR) is a promising choice for thorium utilization because of its deeper burnup, harder spectrum and better inherent safety. Therefore, a 600 MWth 108-column modular high temperature reactor (MHR) core model with 20%UO{sub 2}-80%ThO{sub 2}; MOX fuel is adapted. A traditional two-step calculation scheme, lattice calculation then core calculation, is used to get the neutronic performance of equilibrium cycle including uranium enrichment and mass of fuel. Using these data, the fuel cycle cost of HTR can be calculated by the levelized lifetime cost method. In order to take into account the variation of unit price, a well-established, standard statistical analysis technique is used by Monte Carlo simulation, (thus, the probability distribution of fuel cycle cost can be obtained. Furthermore, sensitivity calculations are made to analyze the impact on the fuel cycle cost of variations in the unit prices for each fuel cycle component. The result shows that the most influential parameters are unit price of natural uranium, enrichment and discount rate. Finally a comparison of fuel cycle cost between 100%UO{sub 2} and 20%UO{sub 2}-80%ThO{sub 2} MOX fuel is made. The latter is 11% cheaper than the former. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Economies of Numbers Formulations for Optimal Process Family Design of Carbon Capture Systems

Enabling rapid, widespread deployment of many process variants across a range of design requirements is key to ensuring the widespread availability of critical decentralized technologies such as carbon capture systems. Unfortunately, traditional process design methods are not efficient for designing many similar variations of a process rapidly across a wide range of system requirements. Rather than considering each process variant as a separate design task, we have proposed viewing the entire set of variants as a “family.” In our work, we present rigorous optimization formulations that simultaneously design each variant within the family from a platform of available unit modules, largely inspired by product family and platform design.

Stinchfield, Georgia↗

DOE OSTI2024

Accelerating Application Bulk Synchronous Writes in HPC Environments

High-bandwidth storage tiers are becoming more common for their capability to absorb high-rate, bursty I/Os. Notably, the designs of these fast storage tiers differ from system to system. The variation of these layers and non-uniform methods of access can pose chal- lenges for applications seeking to run at multiple HPC facilities. Therefore, in this work, we present Spectral, a rapid-output ab- straction library to accelerate application, bulk-synchronous writes on HPC systems. We design Spectral to enable applications to use high-bandwidth storage, such as node-local storage and dis- tributed, write-caches (e.g., burst buffers) transparently without requiring modifications to the application or file system source code. The key idea is to allow applications to spend most of the time performing productive work and to not require any source code changes for maximum portability on different HPC archi- tectures. Spectral internally re-routes write-only files through available, high-performance I/O resources before ultimately mi- grating them to the shared global parallel file system. For instance, on Summit, Spectral transparently places application outputs on node-local storage and then utilizes asynchronous migration to the center-wide GPFS file system. We evaluate Spectral on the Summit HPC system (1024 nodes) using the IOR benchmark and real scientific applications. Spectral shows linear performance scaling, improving application write performance by over an order of magnitude when compared to GPFS.

Khan, Awais↗

DOE OSTI2026

Renormalon subtracted nonrelativistic QCD for heavy hadron systems

We present a renormalon-subtracted formulation of potential nonrelativistic QCD (pNRQCD) for precision spectroscopy of heavy hadron systems, combining variational and Green's function Monte Carlo (VMC/GFMC) methods with NNLO static two- and three-body potentials. Minimal renormalon subtraction (MRS) systematically sums leading factorially growing terms, thereby stabilizing perturbative convergence and reducing renormalization-scale dependence. We tune charm and bottom quark masses to spin-averaged $1S$ quarkonium states and predict $Ω_{ccc}$, $Ω_{ccb}$, $Ω_{cbb}$, and $Ω_{bbb}$ baryon masses, as well as QCD-stable baryons containing top quarks. NNLO MRS results undershoot lattice QCD by 125--175~MeV, with fractional differences decreasing as $\sim 1/m_Q$, consistent with neglected $\mathcal{O}(1/m_Q)$ corrections. Applying these methods to unequal-mass fully-heavy tetraquarks, we determine the critical heavy-to-light mass ratio for binding and compute binding energies across the mass-ratio landscape.

Assi, Benoît [Cincinnati U.; Fermilab]↗

NASA NTRS1994

Computer modeling of the sensitivity of a laser water vapor sensor to variations in temperature and air speed

Currently, there is disagreement among existing methods of determining atmospheric water vapor concentration at dew-points below -40 C. A major source of error is wall effects which result from the necessity of bringing samples into the instruments. All of these instruments also have response times on the order of seconds. NASA Langley is developing a water vapor sensor which utilizes the absorption of the infrared radiation produced by a diode laser to estimate water vapor concentration. The laser beam is directed through an aircraft window to a retroreflector located on an engine. The reflected beam is detected by an infrared detector located near the laser. To maximize signal to noise, derivative signals are analyzed. By measuring the 2f/DC signal and correcting for ambient temperature, atmospheric pressure and air speed (which results in a Doppler shifting of the laser beam), the water vapor concentration can be retrieved. Since this is an in situ measurement there are no wall effects and measurements can be made at a rate of more than 20 per second. This allows small spatial variations of water vapor to be studied. In order to study the sensitivity of the instrument to variations in temperature and air speed, a computer program which generated the 2f, 3f, 4f, DC and 2f/DC signals of the instrument as a function of temperature, pressure and air speed was written. This model was used to determine the effect of errors in measurement of the temperature and air speed on the measured water vapor concentration. Future studies will quantify the effect of pressure measurement errors, which are expected to be very small. As a result of these studied, a retrieval algorithm has been formulated, and will be applied to data taken during the PEM-West atmospheric science field mission. Spectroscopic studies of the water vapor line used by the instrument will be used to refine this algorithm. To prepare for these studies, several lasers have been studied to determine their output frequency range and power.

Tucker, George F.↗

NASA NTRS1984

Maximum likelihood estimation of turbulence spectrum parameters

Estimation of the integral scale and intensity of a generic turbulence record is treated as a statistical problem of parameter estimation. Properties of parameter estimators and the method of maximum likelihood are reviewed. Likelihood equations are derived for estimation of the integral scale and intensity applicable to a general class of turbulence spectra that includes the von Karman and Dryden transverse and longitudinal spectra as special cases. The method is extended to include the Bullen transverse and longitudinal spectra. Coefficients of variation are given for maximum likelihood estimates of the integral scale and intensity of the von Karman spectra. Application of the method is illustrated by estimating the integral scale and intensity of an atmospheric turbulence vertical velocity record assumed to be governed by the von Karman transverse spectrum.

Mark, W. D.↗