SEARCH · Engineering Papers

Results for “interior point 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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1981

A new solution-adaptive grid generation method for transonic airfoil flow calculations

The clustering algorithm is controlled by a second-order, ordinary differential equation which uses the airfoil surface density gradient as a forcing function. The solution to this differential equation produces a surface grid distribution which is automatically clustered in regions with large gradients. The interior grid points are established from this surface distribution by using an interpolation scheme which is fast and retains the desirable properties of the original grid generated from the standard elliptic equation approach.

Nakamura, S.↗

NASA NTRS2001

High Order Numerical Simulation of Sound Generated by the Kirchhoff Vortex

An improved high order finite difference method for low Mach number computational aeroacoustics (CAA) is described. The improvements involve the conditioning of the Euler equations in perturbation form to minimize numerical cancellation error, and the use of a stable non-dissipative sixth-order central spatial differencing for the interior points and third-order at the boundary points. The spatial difference operator satisfies the summation-by-parts property to guarantee strict stability for linear hyperbolic systems. Spurious high frequency oscillations are damped by a third-order characteristic-based filter. The objective of this paper is to apply these improvements in the simulation of sound generated by the Kirchhoff vortex.

Mueller, Bernhard↗

DOE OSTI2022

Computational Methods to Characterize Panel Loading Conditions for Accelerated Testing

Panel cracking and degradation due to wind loading are known to have a detrimental effect on power output. Recreating these damaging conditions in a controlled experimental setting requires an understanding of the loads being generated at the panel surface as a function of both wind speed, wind direction, and panel orientation. To better understand these relationships, a computational fluid dynamics (CFD) simulation package was constructed using an open-source Python library for solving partial differential equations with the finite element method. This CFD package allows the simulated wind speed and panel orientation of a multi-panel array to be easily changed and provides the capability to measure the traction forces at discrete points along the sun-facing and ground-facing surface of an interior panel residing in the wake of one or more upstream panels. A parameter exploration was performed in which the panel angle was varied in even increments from -70 deg to 70 deg and the wind speed was varied from 2 to 30 m/s. During post processing, the measured traction along the panel surface was averaged spatially and interpreted as a time-varying signal, where further processing of this signal yielded the root mean square amplitude and a characteristic frequency associated with the loading. This study found that higher wind speeds are generally associated with increased amplitude of loading and that the panel orientation angle can significantly exacerbate or mitigate this loading. These outcomes are presented along with current work on higher-fidelity verification simulations and recommendations for performing accelerated experimental testing.

loading↗

NASA NTRS1982

Numerical generation of composite three dimensional grids by quasilinear elliptic systems

A three-dimensional grid generation technique designed to numerically construct a boundary-conforming grid within a three-dimensional region bounded by a closed surface is described. The boundary values are generated numerically by a modified elliptic system and are used to compute grid control parameters that are contained in the elliptic systems. The interior grid distribution is governed by the distribution of points on the boundary as well as by the boundary's geometric shape. The composite three-dimensional grid remains both continuous and smooth across the surface of juncture between any two adjoining subregions. The details of the method and its implementation are presented, displaying numerical results for three-dimensional composite grid for a wing-body combination and surface grids.

Thomas, P. D.↗

DOE OSTI2025

Homotopy Solver

This software implements parallel versions of an interior-point solver, based on the publicly available ipopt solver. Here we have full control over the linear solver and our algorithm is fully parallel thus enabling scalability to large-scale optimization problems. This package also has a parallel implementation of a homotopy solver developed under the scalable methods for contact LDRD project 23-ERD-017. This solver is an mfem-based implementation of algorithm described in ``A filter trust-region Newton continuation method for nonlinear complementarity problems''. Cosmin G. Petra, Nai-Yuan Chiang, Jingyi Wang, Tucker Hartland, and Michael Puso (submitted), LLNL-JRNL-869761.

Hartland, Tucker [Lawrence Livermore National Labo↗

NASA NTRS1975

A uniformly valid solution for motion about the interior libration point of the perturbed elliptic-restricted problem

Bounded motion about Lagrange collinear libration points is considered for a perturbed elliptic-restricted problem. A practical application is the motion of a satellite near a libration point collinear with the sun and the earth-moon barycenter. Such a system is treated here as an earth-sun-satellite elliptic restricted problem with lunar perturbations. The method of dual time scales is utilized to develop a uniformly valid three-dimensional analytical solution to the satellite's equations of motion. The analytical development applies somewhat generally to that class of four-body problems where the second primary mass is much greater than the first, and the third primary mass is much greater than the second.

Richardson, D. L.↗

DOE OSTI2022

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗

NASA NTRS1998

A New Approach to Overcoming Spatial Aliasing in Structural Damage Detection

Aircraft, reusable launch vehicles,unmanned aircraft, and other advanced structures are being built using lightweight composite materials/metals with design safety factors as low as 1.25. These advanced structures operate in uncertain and severe environments and are susceptible to damage such as delamination, fiber/matrix damage, hydrothermal strain in composite materials, and fatigue and cracking in metals. To ensure human safety and load-bearing integrity these structures must be inspected to detect and locate often invisible damage and faults before they become catastrophic. Conventional methods of non-destructive evaluation sometimes miss significant damage and are time consuming and expensive to perform. In contrast, vibrometry or vibration signature techniques are a global method of structural integrity monitoring that potentially can efficiently detect damage on large structures, including damage that is away from sensor locations, and in the interior of structures. However, a barrier problem in damage detection using vibration measurements is the need to measure the vibration response at a large number of points on the structure. Typically, model reduction or expansion procedures such as Guyan reduction or dynamic expansion are attempted to overcome the problem of insufficient measurements or spatial aliasing. These approaches depend on using information from the healthy model, and thus put error into the reduction used to represent the damaged structure. In this paper, a Frequency Response Function technique is used to detect damage to a fixed-free beam. The technique uses measured frequency response functions from the healthy structure as reference data, and then monitors vibration measurements during the life of the structure to detect damage. In an analytical simulation using a finite-element model of a beam, damage was located using only sparse measurements because the technique uses both rotation and translation measurements from the damaged structure.

Schulz, Mark J.↗

DOE OSTI2023

Accelerated Aging of PV Cables - The Development of Methods Towards Combined-Accelerated Stress Testing

PV cables allow the collection and distribution of electricity from modules to the energy grid. Durable cabling allows continuous operation of PV installations. Cables with a life less than that of the modules must be replaced - reducing electricity generation and adding to the system operating expense. We have examined the durability of a variety of representative cable samples using artificial UV weathering (IEC TS 62788-7-2) and combined accelerated stress testing (C-AST). Jackets in this study include PV-rated (polyolefin), building wire (PVC/PA), and multi-conductor (rubber, as in PV trackers), with some examples of black and red cables. Characterizations at read points include appearance (camera and optical microscope), surface roughness (linear profilometry), instrumented indentation, surface morphology (SEM), chemical composition (EDS), thermal decomposition (TGA), phase transition (DSC), crystalline structure (WAXS), and polymer chemical structure (FTIR). Complimentary examinations confirm damage to the external and interior layer(s) of less durable samples, including changes in morphology (e.g., cracking), chemistry (e.g., oxidation), structure, and mechanical performance (e.g., through-thickness modulus and hardness). The results may be generally compared to other examples from the field (PV installations) or after accelerated testing (IEC 61215). The methods in this study are presently being applied in a follow-on C-AST study, incorporating both dynamic and kinematic fixtures.

balance of systems↗

NASA NTRS1974

Solution of an eigenvalue problem for the Laplace operator on a spherical surface

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (also referred to as the membrane eigenvalue problem for the vibration equation) on the unit spherical surface are developed. Two specific types of spherical surface domains are considered: (1) the interior of a spherical triangle, i.e., the region bounded by arcs of three great circles, and (2) the exterior of a great circle arc extending for less than pi radians on the sphere (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is derived. The symmetric (generally five-point) finite difference equations that develop are written in matrix form and then solved by the iterative method of point successive overrelaxation. Upon convergence of this iterative method, the fundamental eigenvalue is approximated by iteration utilizing the power method as applied to the finite Rayleigh quotient.

Walden, H.↗

DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2004

Discontinuous Spectral Difference Method for Conservation Laws on Unstructured Grids

A new, high-order, conservative, and efficient discontinuous spectral finite difference (SD) method for conservation laws on unstructured grids is developed. The concept of discontinuous and high-order local representations to achieve conservation and high accuracy is utilized in a manner similar to the Discontinuous Galerkin (DG) and the Spectral Volume (SV) methods, but while these methods are based on the integrated forms of the equations, the new method is based on the differential form to attain a simpler formulation and higher efficiency. Conventional unstructured finite-difference and finite-volume methods require data reconstruction based on the least-squares formulation using neighboring point or cell data. Since each unknown employs a different stencil, one must repeat the least-squares inversion for every point or cell at each time step, or to store the inversion coefficients. In a high-order, three-dimensional computation, the former would involve impractically large CPU time, while for the latter the memory requirement becomes prohibitive. In addition, the finite-difference method does not satisfy the integral conservation in general. By contrast, the DG and SV methods employ a local, universal reconstruction of a given order of accuracy in each cell in terms of internally defined conservative unknowns. Since the solution is discontinuous across cell boundaries, a Riemann solver is necessary to evaluate boundary flux terms and maintain conservation. In the DG method, a Galerkin finite-element method is employed to update the nodal unknowns within each cell. This requires the inversion of a mass matrix, and the use of quadratures of twice the order of accuracy of the reconstruction to evaluate the surface integrals and additional volume integrals for nonlinear flux functions. In the SV method, the integral conservation law is used to update volume averages over subcells defined by a geometrically similar partition of each grid cell. As the order of accuracy increases, the partitioning for 3D requires the introduction of a large number of parameters, whose optimization to achieve convergence becomes increasingly more difficult. Also, the number of interior facets required to subdivide non-planar faces, and the additional increase in the number of quadrature points for each facet, increases the computational cost greatly.

Liu, Yen↗

DOE OSTI2024

Shifts in belowground processes along a temperate forest edge

Abstract Context Forests are increasingly fragmented, and as a result most forests in the United States are within one km of an edge. Edges change environmental conditions of the forest—especially radiation, roughness, temperature, and moisture—that can have consequences for plant productivity and ecosystem functions. However, edge effects on aboveground characteristics of plants and the environment are better understood relative to plant roots and soil in the belowground environment. Objectives Our main objectives were to determine if soil C pools and fluxes are higher at the edge relative to other landscape positions, and to understand how specific belowground processes contribute to bulk differences in pools and fluxes. Methods We measured environmental conditions, live and dead fine root traits, soil chemistry, and soil respiration along a 75 m transect from interior forest to meadow in Gaithersburg, MD. Results We observed differences in the soil chemical, biological and hydrological environment between the forest interior, edge and adjacent meadow. In some cases, the forest edge represented a mid-point in environmental or belowground characteristics between the forest interior and meadow ( e.g. , pH, C-to-N ratio [C:N], live fine root biomass, heterotrophic respiration), likely reflecting the change in litter type and quality associated with the transition from grass to woody species. In other cases, neighboring landscape positions were different from the forest edge, which was drier and had higher dead fine root biomass. Although soil C contents were not significantly different across landscape positions, there was a tendency towards higher average soil C content at the edge relative to other landscape positions, suggesting that increased C loss related to root decay and greater soil respiration at the edge relative to the forest interior may have been offset by increased C gain from high plant productivity and subsequent inputs to soil. Conclusions This research provides insight into how forest edge environments may differ from the interior and how concurrent processes above- and belowground may contribute to those differences.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

NASA NTRS1985

Prediction of Acoustical Response of Three-dimensional Cavities Using an Indirect Boundary Element Method

A technique which characterizes the acoustics of generalized cavities with the minimum model possible is developed. All boundary element methods have two advantages over finite element methods: (1) the models are smaller; and (2) the assumed variable behavior, inherent in the method to allow discretization, is harmonic rather than polynomial. Further, IBEM often requires one rather than two numerical boundary integrals as required by DBEM. Thus, a quadratic, isoparametric IBEM program was developed. The source distribution in this solution is continuous and quadratically variable rather than continuous and constant. The program was also formulated to include the additional capability of interior point sources and impedance boundary conditions. To test the quadratic, isoparametric IBEM program, several simple cavity enclosure problems where studied. Results are shown.

Bernhard, R. J.↗

DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

Microscopic stardust grains, found in trace amounts in primitive meteorites, afford a unique opportunity to study stellar nucleosynthesis in the laboratory by measuring their heavy element isotopic compositions. Applying LLNL’s state-of-the art material characterization capabilities to these micrometer-sized grains yields large new datasets with unprecedented precision. However, the astrophysical models necessary to interpret the isotopic record require improved nuclear physics data to constrain stellar nucleosynthesis conditions. For example, the branch point in the s-process path at 95 Zr controls the relative amounts of heavier nuclides but is poorly understood. We conducted a coordinated effort across disciplines to combine a new determination of the neutron capture cross section of 95 Zr with the laboratory analysis of hundreds of stardust grains to yield a view into stellar interiors. We developed new methods to quantify multiple isotopic systems simultaneously (e.g., Zr, Ba, and W) in stardust grains and improved the determination of the neutron capture cross section of unstable 95 Zr via the surrogate reaction method. This interdisciplinary approach is allowing us to fill gaps in our understanding of the environments in which elements are produced while enhancing LLNL’s experimental and theoretical capabilities relevant to nuclear threat reduction missions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2023

Deciphering Fingerprints of Stellar Nucleosynthesis Through Nuclear Reaction Rate Measurements and Isotopic Analyses of Stardust

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

NASA NTRS1992

Experimental verification of the asymtotic modal analysis method as applied to a rectangular acoustic cavity excited by structural vibration

An experiment was performed on a rigid wall rectangular acoustic cavity driven by a flexible plate mounted in a quarter of one end wall and excited by white noise. The experiment was designed so that the assumptions of Asymptotic Modal Analysis (AMA) were satisfied for certain bandwidths and center frequencies. Measurements of sound pressure levels at points along the boundaries and incrementally into tbe interior were taken. These were compared with the theoretical results predicted with AMA, and found to be in good agreement, particularly for moderate (1/3 octave) bandwidths and sufficiently high center frequencies. Sound pressure level measurements were also taken well into the cavity interior at various points along the 5 totally rigid walls. The AMA theory, including boundary intensification effects, was shown to be accurate provided the assumption of large number of acoustic modes is satisfied, and variables such as power spectra of the wall acceleration, frequency, and damping are slowly varying in the frequency of bandwidth.

Peretti, L. F.↗