Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cauchy”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Effects of hindlimb unweighting on the mechanical and structure properties of the rat abdominal aorta

Previous studies have shown that hindlimb unweighting of rats, a model of microgravity, reduces evoked contractile tension of peripheral conduit arteries. It has been hypothesized that this diminished contractile tension is the result of alterations in the mechanical properties of these arteries (e.g., active and passive mechanics). Therefore, the purpose of this study was to determine whether the reduced contractile force of the abdominal aorta from 2-wk hindlimb-unweighted (HU) rats results from a mechanical function deficit resulting from structural vascular alterations or material property changes. Aortas were isolated from control (C) and HU rats, and vasoconstrictor responses to norepinephrine (10(-9)-10(-4) M) and AVP (10(-9)-10(-5) M) were tested in vitro. In a second series of tests, the active and passive Cauchy stress-stretch relations were determined by incrementally increasing the uniaxial displacement of the aortic rings. Maximal Cauchy stress in response to norepinephrine and AVP were less in aortic rings from HU rats. The active Cauchy stress-stretch response indicated that, although maximum stress was lower in aortas from HU rats (C, 8.1 +/- 0.2 kPa; HU, 7.0 +/- 0.4 kPa), it was achieved at a similar hoop stretch. There were also no differences in the passive Cauchy stress-stretch response or the gross vascular morphology (e.g., medial cross-sectional area: C, 0.30 +/- 0.02 mm(2); HU, 0.32 +/- 0.01 mm(2)) between groups and no differences in resting or basal vascular tone at the displacement that elicits peak developed tension between groups (resting tension: C, 1.71 +/- 0.06 g; HU, 1.78 +/- 0.14 g). These results indicate that HU does not alter the functional mechanical properties of conduit arteries. However, the significantly lower active Cauchy stress of aortas from HU rats demonstrates a true contractile deficit in these arteries.

Non-NASA Center↗
NASA NTRS2006

Oscillation Amplitude Growth for a Decelerating Object with Constant Pitch Damping

The equations governing the deceleration and oscillation of a blunt body moving along a planar trajectory are re-expressed in the form of the Euler-Cauchy equation. An analytic solution of this equation describes the oscillation amplitude growth and frequency dilation with time for a statically stable decelerating body with constant pitch damping. The oscillation histories for several constant pitch damping values, predicted by the solution of the Euler-Cauchy equation are compared to POST six degree-of-freedom (6-DoF) trajectory simulations. The simulations use simplified aerodynamic coefficients matching the Euler-Cauchy approximations. Agreement between the model predictions and simulation results are excellent. Euler-Cauchy curves are also fit through nonlinear 6-DoF simulations and ballistic range data to identify static stability and pitch damping coefficients. The model os shown to closely fit through the data points and capture the behavior of the blunt body observed in simulation and experiment. The extracted coefficients are in reasonable agreement with higher fidelity, nonlinear parameter identification results. Finally, a nondimensional version of the Euler-Cauchy equation is presented and shown to be a simple and effective tool for designing dynamically scaled experiments for decelerating blunt capsule flight.

Schoenenberger, Mark↗
DOE OSTI2025

Informed unsupervised machine learning analysis of dislocation microstructure from high-resolution differential aperture X-ray structural microscopy data

This study leverages high-resolution differential-aperture X-ray structural microscopy (DAXM) to probe the local dislocation structure in deformed 304L-stainless steel at small strain, by measuring the lattice rotation and deviatoric elastic strain with a sub-micron resolution. For a single grain in a polycrystalline specimen, the measured lattice rotation field over the measured volume exhibited a multimodal distribution while the deviatoric elastic strain showed a single-mode distribution. An unsupervised Cauchy mixture machine learning model was developed to resolve the multimodal distribution of the lattice rotation. By mapping the lattice rotation data associated with each Cauchy peak in the model back onto the measured volume, we identify contiguous regions of the crystal rotated near the average values corresponding to the peaks of the overall rotation distribution. These regions represent the grain subdivision in the microstructure. Finally, the dislocation density tensor was also computed and its norm was laid over the rotation field to detect the subgrain boundaries. This step provided a validation of the Cauchy mixture model for the analysis of the lattice rotation distribution. The current study highlights the integration of advanced X-ray microscopy techniques with data-driven analysis methods to uncover detailed microstructure scales in deformed crystals.

Machine learning; Lattice rotation; High-energy X-↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
NASA NTRS1976

Arrangement of vortex lattices

A new method is developed for solving the lifting-surface equation for thin wings. The solution requires the downwash equation to be in the form of Cauchy integrals which can be interpreted as a vortex lattice with the positions of the vortices and control points dictated by the finite sum used to approximate the integrals involved. Lan's continuous loading method is employed for the chordwise integral since it properly accounts for the leading-edge singularity, Cauchy singularity, and Kutta condition. Unlike Lan, the spanwise loading is also continuous and the Cauchy singularity in the spanwise integral is also properly accounted for by using the midpoint trapezoidal rule and the theory of Chebyshev polynomials. This technique yields the exact classical solution to Prandtl's lifting-line equation.The solution to the lifting-surface equation for rectangular wings was found to compare well with other continuous loading methods, but with much smaller computational times, and to converge faster than other vortex lattice methods.

Dejarnette, F. R.↗
NASA NTRS1976

Lifting surface theory for rectangular wings

A new incompressible lifting-surface theory is developed for thin rectangular wings. The solution requires the downwash equation to be in the form of Cauchy-type integrals. Lan's method is employed for the chordwise integrals since it properly accounts for the leading-edge singularity, Cauchy singularity and Kutta condition. The Cauchy singularity in the spanwise integral is also accounted for by using the midpoint trapezoidal rule and theory of Chebychev polynomials. The resulting matrix equation, formed by satisfying the boundary condition at control points, is simpler and quicker to compute than other lifting surface theories. Solutions were found to converge with only a small number of control points and to compare favorably with results from other methods.

Dejarnette, F. R.↗
DOE OSTI2024

Constrained HRT Surfaces and their Entropic Interpretation

Abstract Consider two boundary subregionsAandBthat lie in a common boundary Cauchy surface, and consider also the associated HRT surfaceγ B forB. In that context, the constrained HRT surfaceγ A:B can be defined as the codimension-2 bulk surface anchored toAthat is obtained by a maximin construction restricted to Cauchy slices containingγ B . As a result,γ A:B is the union of two pieces,$$ {\gamma}_{A:B}^B $$ γ A : B B and$$ {\gamma}_{A:B}^{\overline{B}} $$ γ A : B B ¯ lying respectively in the entanglement wedges ofBand its complement$$ \overline{B} $$ B ¯ . Unlike the area$$ \mathcal{A}\left({\gamma}_A\right) $$ A γ A of the HRT surfaceγ A , at least in the semiclassical limit, the area$$ \mathcal{A}\left({\gamma}_{A:B}\right) $$ A γ A : B ofγ A:B commutes with the area$$ \mathcal{A}\left({\gamma}_B\right) $$ A γ B ofγ B . To study the entropic interpretation of$$ \mathcal{A}\left({\gamma}_{A:B}\right) $$ A γ A : B , we analyze the Rényi entropies of subregionAin a fixed-area state of subregionB. We use the gravitational path integral to show that then ≈1 Rényi entropies are then computed by minimizing$$ \mathcal{A}\left({\gamma}_A\right) $$ A γ A over spacetimes defined by a boost angle conjugate to$$ \mathcal{A}\left({\gamma}_B\right) $$ A γ B . In the case where the pieces$$ {\gamma}_{A:B}^B $$ γ A : B B and$$ {\gamma}_{A:B}^{\overline{B}} $$ γ A : B B ¯ intersect at a constant boost angle, a geometric argument shows that then ≈1 Rényi entropy is then given by$$ \frac{\mathcal{A}\left({\gamma}_{A:B}\right)}{4G} $$ A γ A : B 4 G . We discuss how then ≈1 Rényi entropy differs from the von Neumann entropy due to a lack of commutativity of then→ 1 andG →0 limits. We also discuss how the behaviour changes as a function of the width of the fixed-area state. Our results are relevant to some of the issues associated with attempts to use standard random tensor networks to describe time dependent geometries.

Physics↗
DOE OSTI2020

Bond-Based Peridynamics: a Tale of Two Poisson’s Ratios

This paper explores the restrictions imposed by bond-based peridynamics, particularly with respect to plane strain and plane stress models. We begin with a review of the derivations in Gerstle et al. (2005) wherein for isotropic materials a Poisson’s ratio restriction of $\frac{1}{4}$ for plane strain and $\frac{1}{3}$ for plane stress is deduced. Next, we show Cauchy’s relations are an intrinsic limitation of bond-based peridynamics and specialize these restrictions to plane strain and plane stress models. This generalizes the results from Gerstle et al. (2005) and demonstrates that the Poisson’s ratio restrictions described in Gerstle et al. (2005) are merely a consequence of Cauchy’s relations for isotropic materials. We conclude with a discussion of the validity of peridynamic plane strain and plane stress models formulated from two-dimensional bond-based peridynamic models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Application of the locally self-consistent embedding approach to the Anderson model with non-uniform random distributions

Highlights: • Typical Medium Theory (TMT) for the Anderson Localization. • Locally Self-Consistent Multiple Scattering Method (LSMS) for Random Disordered Systems. • Linear Scaling Computational Method for Random Disordered Systems. We apply the recently developed embedding scheme for the locally self-consistent method to random disorder electrons systems. The method is based on the locally self-consistent multiple scattering theory and the typical medium theory. The locally self-consistent multiple scattering theory divides a system into many small designated local interaction zones. The subsystem within each local interaction zone is embedded in a self-consistent field from the typical medium theory. This approximation allows the study of random systems with large numbers of sites. We present results for the three dimensional Anderson model with different random disorder potential distributions. Using the typical density of states as an indicator of Anderson localization, we find that the method can capture the localization for commonly studied disorder potentials. These include the uniform distribution, the Gaussian distribution, and even the unbounded Cauchy distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Robust Online Sequential RVFLNs for Data Modeling of Dynamic Time-Varying Systems with Application of an Ironmaking Blast Furnace

In a world where the increasing complexity of modern industrial processes brings difficulties for accurate mathematical modeling, taking advantage of data has become an efficient solution to complex dynamic process modeling issue. In this paper, we develop a novel robust online sequential version of random vector functional-link networks (RVFLNs) for data-driven modeling of dynamic time-varying system and applied it in a blast furnace (BF) ironmaking process. First, to overcome the time-varying dynamics of process and to enable the RVFLNs to learn online with avoiding data saturation, an improved online sequential version of RVFLNs (OS-RFVLNs) is first presented by online sequential learning with forgetting factor. This improved OS-RVFLNs algorithm is not only suitable for the real-time and large data transfer situation, but also can adjust the sensitivity of the algorithm to different samples with the help of the introduced forgetting factor. Second, since the output weights of the improved OS-RVFLNs as well as other RVFLNs algorithms are obtained by the least squares approach, a robustness problem may occur when the training dataset is contaminated with various outliers. To solve this problem, a Cauchy distribution weighted M-estimator is introduced to improve the robustness of the improved OS- RVFLNs. For this proposed robust OS-RVFLNs (R-OS- RVFLNs), since the weights of different outlier data are properly determined by the Cauchy distribution function, their corresponding contribution on modeling can be properly distinguished. Thus robust and better modeling results can be achieved. Experiments using actual industrial data of BF ironmaking process and comparative studies have demonstrated that the proposed method produces a better estimation accuracy and stronger robustness than other methods.

Blast furnace (BF), Modelling, Dynamic systems↗
DOE OSTI2024

Diffeomorphism invariance and quantum mechanical paradoxes

Paradoxes in gravitational physics, such as grandfather paradoxes with closed timelike curves or the AMPS paradox, are often constructed in a weak gravity regime but upon further examination engender strong gravitational responses such as unstable Cauchy horizons or firewalls. In contrast, some proposed paradoxes in nonrelativistic quantum mechanics ignore gravity completely. Such nonrelativistic proposed paradoxes are often not gauge invariant — they use local operators which should be forbidden by diffeomorphism-invariant quantum gravity. Ignoring this complication is reasonable if there is a consistent weak gravity regime where the paradox can be formulated with gauge invariant observables up to some order in Newton’s constant. Here, we show that this approach remains inconsistent due to a lack of analyticity of the necessary solutions. As a consequence, there is no consistent weak gravity, diffeomorphism invariant embedding of such paradoxes—just as in other gravitational paradoxes they generate a strong gravitational response and are in principle sensitive to quantum gravity.

quantum gravity↗
DOE OSTI2022

Exponential time differencing scheme for mass transport and depletion in molten salt reactors

This work extends the capability previously shown for addressing the problem of computing depletion and mass transport calculations in molten salt reactors (MSRs) by calculating matrix exponentials. Additional algorithms are implemented to compute the matrix exponential and the action of the matrix exponential on a matrix. These algorithms include two methods based on the Pade approximation, a Taylor series method, and three methods based on Cauchy's integral formula. In addition to the added matrix exponential solvers, a variable-order total variation diminishing scheme is applied to the convective flux approximation to provide enhanced accuracy. Finally, a simplified MSR problem is shown for each of the exponential time differencing solvers along with classical backwards differencing integrators. The results show excellent convergence for exponential time differencing methods. Computation time is a key element for selecting the optimal solver in these problems, and this work shows that Pade and Cauchy-based solvers may provided the fastest and most accurate solutions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Pressure dependence of Poisson’s ratio of glassy Baltic amber studied by Brillouin scattering spectroscopy

Amber is a typical natural glass with very long aging time. Elastic properties of amber at high pressures have been studied using Brillouin scattering and a diamond anvil cell. Both longitudinal and transverse acoustic modes have been observed up to 12 GPa. The pressure dependences of longitudinal, shear, Young’s, and bulk moduli, compressibility, and Poisson’s ratio were determined. The longitudinal, shear, Young’s, bulk moduli show the remarkable increase, and compressibility shows a marked decrease with increasing pressure. However, it is found that the pressure dependence of Poisson’s ratio is very small. The mechanism of this small pressure dependence was discussed. Finally, the Cauchy type relation between longitudinal and shear moduli was examined. Its coefficients show the deviation from the Cauchy relation owing to the coexistence of different intermolecular interactions in amber.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS1972

Subsonic, transonic, and supersonic nozzle flow by the inverse technique.

The inverse technique is used to obtain a mathematically and physically consistent solution of the flowfield in a nozzle from the mass generation surface through the supersonic region. The inverse method employs an assumed centerline function which is of the Cauchy type in that the values and the derivatives of the function are known. Since the Cauchy boundary conditions can give rise to numerical instabilities, the governing gasdynamic equations for rotational steady flow were transformed into a form which puts the geometry into a rectangular shape, and which spaces the network of interior points more finely in regions of the greatest gradients of the dependent variables. For arbitrarily specified centerline data, the solution of the governing flow equations may not exist, and if it does it may not depend continuously on the data.

Norton, D. J.↗
NASA NTRS1976

Stress-intensity factors for a thick-walled cylinder containing an annular imbedded or external or internal surface crack

The elastostatic axisymmetric problem for a long thick-walled cylinder containing a ring-shaped internal or edge crack is considered. Using the standard transform technique the problem is formulated in terms of an integral equation which has a simple Cauchy kernel for the internal crack and a generalized Cauchy kernel for the edge crack as the dominant part. As examples the uniform axial load and the steady-state thermal stress problems have been solved and the related stress intensity factors have been calculated. Among other findings the results show that in the cylinder under uniform axial stress containing an internal crack the stress intensity factor at the inner tip is always greater than that at the outer tip for equal net ligament thicknesses and in the cylinder with an edge crack which is under a state of thermal stress the stress intensity factor is a decreasing function of the crack depth, tending to zero as the crack depth approaches the wall thickness.

Erdol, R.↗
NASA NTRS1978

A thick-walled cylinder with an axisymmetric internal or edge crack

The paper considers the elastostatic axisymmetric problem for a long thick-walled cylinder containing a ring-shaped internal or edge crack. Using the standard transform technique the problem is formulated in terms of an integral equation which has a simple Cauchy kernel for the internal crack and a generalized Cauchy kernel for the edge crack as the dominant part. As examples the uniform axial load and the steady-state thermal stress problems have been solved and the related stress-intensity factors have been calculated. Among other findings the results show that in the cylinder under uniform axial stress containing an internal crack the stress-intensity factor at the inner tip is always greater than that at the outer tip for equal net ligament thicknesses and in the cylinder with an edge crack which is under a state of thermal stress the stress-intensity factor is a decreasing function of the crack depth, tending to zero as the crack depth approaches the wall thickness.

Erdol, R.↗
NASA NTRS1982

Analysis of airfoil leading edge separation bubbles

A local inviscid-viscous interaction technique was developed for the analysis of low speed airfoil leading edge transitional separation bubbles. In this analysis an inverse boundary layer finite difference analysis is solved iteratively with a Cauchy integral representation of the inviscid flow which is assumed to be a linear perturbation to a known global viscous airfoil analysis. Favorable comparisons with data indicate the overall validity of the present localized interaction approach. In addition numerical tests were performed to test the sensitivity of the computed results to the mesh size, limits on the Cauchy integral, and the location of the transition region.

Carter, J. E.↗