Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational discretization”

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 37 records · Page 2

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An inexact semismooth Newton method with application to adaptive randomized sketching for dynamic optimization

In many applications, one can only access the inexact gradients and inexact hessian times vector products. Thus it is essential to consider algorithms that can handle such inexact quantities with a guaranteed convergence to solution. An inexact adaptive and provably convergent semismooth Newton method is considered to solve constrained optimization problems. In particular, dynamic optimization problems, which are known to be highly expensive, are the focus. A memory efficient semismooth Newton algorithm is introduced for these problems. The source of efficiency and inexactness is the randomized matrix sketching. Further, applications to optimization problems constrained by partial differential equations are also considered.

97 MATHEMATICS AND COMPUTING↗

Improving reproducibility in synchrotron tomography using implementation-adapted filters

For reconstructing large tomographic datasets fast, filtered backprojection-type or Fourier-based algorithms are still the method of choice, as they have been for decades. These robust and computationally efficient algorithms have been integrated in a broad range of software packages. The continuous mathematical formulas used for image reconstruction in such algorithms are unambiguous. However, variations in discretization and interpolation result in quantitative differences between reconstructed images, and corresponding segmentations, obtained from different software. This hinders reproducibility of experimental results, making it difficult to ensure that results and conclusions from experiments can be reproduced at different facilities or using different software. In this paper, a way to reduce such differences by optimizing the filter used in analytical algorithms is proposed. These filters can be computed using a wrapper routine around a black-box implementation of a reconstruction algorithm, and lead to quantitatively similar reconstructions. Use cases for this approach are demonstrated by computing implementation-adapted filters for several open-source implementations and applying them to simulated phantoms and real-world data acquired at the synchrotron. Our contribution to a reproducible reconstruction step forms a building block towards a fully reproducible synchrotron tomography data processing pipeline.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Time Series Dimension Reduction for Surrogate Models of Port Scanning Cyber Emulations

Surrogate model development is a key resource in the scientific modeling community for providing computational expedience when simulating complex systems without loss of great fidelity. The initial step to development of a surrogate model is identification of the primary governing components of the system. Principal component analysis (PCA) is a widely used data science technique that provides inspection of such driving factors, when the objective for modeling is to capture the greatest sources of variance inherent to a dataset. Although an efficient linear dimension reduction tool, PCA makes the fundamental assumption that the data is continuous and normally distributed. Thus, it provides ideal performance when these conditions are met. In the case for which cyber emulations provide realizations of a port scanning scenario, the data to be modeled follows a discrete time series function comprised of monotonically increasing piece-wise constant steps. The sources of variance are related to the timing and magnitude of these steps. Therefore, we consider using XPCA, an extension to PCA for continuous and discrete random variates. This report provides the documentation of the trade-offs between the PCA and XPCA linear dimension reduction algorithms, for the intended purpose to identify key components of greatest variance in our time series data. These components will ultimately provide the basis for future surrogate models of port scanning cyber emulations.

97 MATHEMATICS AND COMPUTING↗

The UCSD/MIT hard X-ray and low energy gamma ray experiment for HEAO-A

An HEAO-A spacecraft is described which is scheduled for launch in April 1977 to survey the sky in the 0.15 keV to 10 MeV range of X-rays with improved sensitivity and angular resolution. The spectral and time variations of discrete sources, and the spectrum and isotropy of the diffuse component will be measured. A hard X-ray and low energy gamma ray experiment is described which will use as an array of seven (NaI)(Tl)/CsI(Na) phoswich detectors and eight large CsI(Na) shields. Detector events are pulse-height analyzed with a selectable time resolution of 0.1, 0.005, or 0.00002 sec. The instrument background is expected to be comprised of approximately equal parts of diffuse gamma rays entering the aperture, radioactive decay of spallation products, and gamma ray leakage through the shields. Detailed background predictions are presented, and sensitivity estimates indicate that the Crab nebula should be detected up to approximately 5 MeV. Approximate 1% isotropy limits for the diffuse component of cosmic gamma rays should be achieved, and its spectrum measured with improved precision.

Matteson, J. L.↗

The UCSD/MIT hard X-ray and low energy gamma-ray experiment for HEAO-1 - Design and early results

Primary scientific objectives of the experiment are the measurement of the spectra and time variations of discrete sources and the spectrum and isotropy of the diffuse background, all in the 10 keV to 10 MeV energy range. The instrument consists of an array of 7 NaI(Tl)/CsI(Na) phoswich-type detectors and 8 large CsI(Na) shields. In order to achieve the maximum sensitivity allowed by weight, data rate and funding constraints of the mission, three detector and aperture geometries were selected. The instrument's total sensitive area is 494 sq cm. The fields of view vary from 1.7 deg x 20 deg FWHM at low energies to 43 deg FWHM at high energies. The basic instrument concepts, details of the mechanical and electrical design, and pre-launch calibrations are presented. In-orbit performance is described and background measurements are analyzed and used to estimate the experiment's ultimate sensitivity.

Matteson, J. L.↗

CryoScout: A Descent Through the Mars Polar Cap

CryoScout was proposed as a subsurface investigation of the stratigraphic climate record embedded in Mars North Polar cap. After landing on a gentle landscape in the midst of the mild summer season, CryoScout was to use the continuous polar sunlight to power the descent of a cryobot, a thermal probe, into the ice at a rate of about 1 m per day. CryoScout would probe deep enough into this time capsule to see the effects of planetary obliquity variations and discrete events such as dust storms or volcanic eruptions. By penetrating tens of meters of ice, the mission would explore at least one of the dominant "MOC layers" observed in exposed layered terrain.

Hecht, M. H.↗

Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle

Within this work, we establish a stochastic contact variational integrator and its discrete version via stochastic Herglotz variational principle for stochastic contact Hamiltonian systems. A general structure-preserving stochastic contact method is provided to seek the stochastic contact variational integrators. Numerical experiments are performed to verify the validity of this approach.

97 MATHEMATICS AND COMPUTING↗

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Total Variation of High-Order Semi-Discrete Central Schemes for Conservation Laws

We discuss a new fifth-order, semi-discrete, central-upwind scheme for solving one-dimensional systems of conservation laws. This scheme combines a fifth-order WENO reconstruction, a semi-discrete central-upwind numerical flux, and a strong stability preserving Runge-Kutta method. We test our method with various examples, and give particular attention to the evolution of the total variation of the approximations.

Bryson, Steve↗

A Three-Dimensional, Thermodynamically and Variationally Consistent, Fully Coupled, Electro-Chemo-Thermo-Mechanical Model of Solid-State Batteries

Here we present a theoretical framework for high fidelity modeling of three-dimensional electro-chemo-thermo-mechanical coupled systems such as solid-state batteries. The framework provides a useful means of utilizing molecular level information to predict the overall performance of the system at the continuum level. The central statement of the theoretical framework is a thermodynamically consistent inf-sup problem. The total potential of the coupled system consists of the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. All governing equations derive from variational principles. We also show that the resulting initial-boundary-value problem (IBVP) is consistent with existing, widely adopted equations for modeling batteries. Moreover, the IBVP can be naturally discretized using the (variational) finite element method to simulate practical solid-state batteries with arbitrary geometries. Several examples illustrate multiphysics coupling between pairs of the physical field equations, culminating in a comprehensive three-dimensional model of charge-discharge in a solid state battery.

25 ENERGY STORAGE↗

Accurate interlaminar stress recovery from finite element analysis

The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Tessler, Alexander↗

Inventory of forest and rangeland and detection of forest stress

The author has identified the following significant results. At the Atlanta site (226B) it was found that bulk color composites for October 15, 1972, and April 13, 1973, can be interpreted together to disclose the location of the perennial Kudzu vine (Pyeraria lobata). Land managers concerned with Kudzu eradication could use ERTS-1 to inventory locations over 200 meters (660 feet) square. Microdensitometer data collected on ERTS-1 Bulk photographic products for the Manitou test site (226C) have shown that the 15-step gray-scale tablets are not of systematic equal values corresponding to 1/14 the maximum radiant energy incident on the MSS sensor. The gray-scale values present a third-order polynomial function rather than a direct linear relationship. Although data collected on step tablets for precision photographic products appear more discrete, the density variation within blocks in almost as great as variations between blocks. These system errors will cause problems when attempting to analyze radiometric variances among vegetation and land use classes.

Heller, R. C.↗

Coherent substructure of turbulence near the stagnation zone of a bluff body

The evolution of freestream turbulence in crossflow about a circular cylinder was studied in order to identify the existence of a coherent substructure which is the outcome of the amplification of freesteam turbulence by the stretching mechanism in diverging flow about a bluff body. Visualization of the flow events revealed the selective stretching of cross-vortex tubes and the emergence of an organized turbulent flow pattern near the cylinder stagnation zone. Significant amplification of the total turbulent energy of the streamwise fluctuating velocity was consistently monitored. Realization of selective amplification at scales larger than the neutral scale of the stagnation flow was indicated by the variation of the discrete streamwise turbulent energy. A most amplified scale, characteristic of the energy containing eddies within the coherent substructure and commensurate with the boundary-layer thickness, was detected. Penetration of the amplified turbulence into the cylinder boundary layer led to the retardation of separation and to a concurrent decrease in the drag coefficient at subcritical cylinder-diameter Reynolds numbers.

Sadeh, W. Z.↗

Turbulence amplification in flow about an airfoil

An experimental investigation of the evolution of freestream turbulence in flow about an airfoil was conducted in order to ascertain its selective amplification induced by the stretching mechanism according to the vorticity-amplification theory. Significant amplification of the streamwise turbulent energy transpired even in the limiting flow situation studied of a symmetric airfoil at zero angle of attack where the stretching is the least. Substantiation of the stretching effect was provided by the almost 100 percent amplification of turbulence with respect to its background level in the absence of the airfoil. Realization of preferred amplification at scales larger than the neutral scale of the stagnation flow was clearly indicated by the variation of the discrete streamwise turbulent energy. Particularly important was the detection of a most amplified scale which is characteristic of the coherent substructure near the airfoil stagnation zone and, concurrently, commensurate with the boundary-layer thickness.

Sadeh, W. Z.↗

On the solution of problems involving impact type loading

The response of beam and plate members to pulse and impact loading is investigated, and numerical methods are analyzed and compared. The equations of motion, the finite difference and finite element methods, the time integration/Runge-Kutta techniques, and material modeling are discussed in detail. It is found that both the finite difference and the finite element methods could be accurately employed to discretize the spatial variation in the displacements. Central differences or a fourth order Runge-Kutta algorithm could be used for the time integration. The total energy of the system would give the stability and accuracy of the solution. Results showed that the finite element method provided better efficiency in obtaining accurate solution than the finite difference method when a scalar processor is used. However, the finite difference method was more efficient on the vector processor. Therefore, the most efficient method of solution depends on the type of computer present for the analysis.

Moyer, E. T., Jr.↗