Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Analysis of Auroral Morphology: Substorm Precursor and Onset on January 10, 1997

The solar wind interaction with the geomagnetic field is studied using global auroral images obtained by the Ultraviolet Imager (UVI) on Polar. We study the dynam,cs of the poleward and equatorward boundaries of the auroral oval in response to the solar wind IMF on January 10, 1997 using a neural network algorithm to perform an automated morphological analysis. Poleward and equatorward boundaries identified by the algorithm demonstrate a clear growth motion with the southward turning of the IMF and growth and poleward expansion at substorm onset. The area poleward of the oval (polar cap) is found to increase in size coincident with the'southward turning of the IMF Bz component at 0220 UT and peaks at substorm onset at 0334 UT. The area of the oval, however, decreases continuously throughout the period of the polar cap area increase with a slight recovery observed during the substorm onset. These observations are consistent with the concept that magnetospheric dynamics are directly driven by the solar wind-geomagnetic field interactions.

Germany, G. A.↗
DOE OSTI2020

Microsecond Molecular Dynamics Simulations of Proteins Using a Quasi-Equilibrium Solvation Shell Model

In this work, we describe the development and implementation of a quasi-equilibrium hydration shell model of biomolecular solvation with adaptive boundaries. Applying the model to microsecond-long molecular dynamics simulations of several protein systems of varying complexity, we find that the model simulation results are of comparable quality to those obtained from simulations of fully solvated systems, but at a reduced computational cost. We discuss the dominant sources of error in the model and outline directions for future improvements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗
DOE OSTI2020

Dislocation dynamics in polycrystalline materials

An extension of the dislocation dynamics method to polycrystalline materials is presented with emphasis on handling interactions between dislocations and grain boundaries. The key new features of the method include: (1) a grain boundary detection algorithm, (2) new dislocation emission criteria based on power dissipation and area growth, (3) mobility of dislocations at grain boundaries, and (4) extension of topological operations e.g. mesh adaption, collision, and dissociation to handle dislocations interacting with grain boundaries. Furthermore, verification and validation examples are executed to compare our calculations with existing analytic methods, molecular dynamics simulations, and experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗
DOE OSTI2026

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗
DOE OSTI2023

The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere

Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Refining Planetary Boundary Layer Height Retrievals From Micropulse‐Lidar at Multiple ARM Sites Around the World

Abstract Knowledge of the planetary boundary layer height (PBLH) is crucial for various applications in atmospheric and environmental sciences. Lidar measurements are frequently used to monitor the evolution of the PBLH, providing more frequent observations than traditional radiosonde‐based methods. However, lidar‐derived PBLH estimates have substantial uncertainties, contingent upon the retrieval algorithm used. In addressing this, we applied the Different Thermo‐Dynamic Stabilities (DTDS) algorithm to establish a PBLH data set at five separate Department of Energy's Atmospheric Radiation Measurement sites across the globe. Both the PBLH methodology and the products are subject to rigorous assessments in terms of their uncertainties and constraints, juxtaposing them with other products. The DTDS‐derived product consistently aligns with radiosonde PBLH estimates, with correlation coefficients exceeding 0.77 across all sites. This study delves into a detailed examination of the strengths and limitations of PBLH data sets with respect to both radiosonde‐derived and other lidar‐based estimates of the PBLH by exploring their respective errors and uncertainties. It is found that varying techniques and definitions can lead to diverse PBLH retrievals due to the inherent intricacy and variability of the boundary layer. Our DTDS‐derived PBLH data set outperforms existing products derived from ceilometer data, offering a more precise representation of the PBLH. This extensive data set paves the way for advanced studies and an improved understanding of boundary‐layer dynamics, with valuable applications in weather forecasting, climate modeling, and environmental studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Particle Dynamics in a Diblock-Copolymer-Based Dodecagonal Quasicrystal and Its Periodic Approximant by X-Ray Photon Correlation Spectroscopy

Temperature-dependent X-ray photon correlation spectroscopy (XPCS) measurements are reported for a binary diblock copolymer blend that self-assembles into an aperiodic dodecagonal quasicrystal (DDQC) and a periodic Frank-Kasper σ phase approximant. The measured structural relaxation times are Bragg scattering wavevector-independent and are five times faster in the DDQC than the σ phase, with minimal temperature dependence. Further, the underlying dynamical relaxations are ascribed to differences in particle motion at the grain boundaries within each of these tetrahedrally close-packed assemblies. These results identify unprecedented particle dynamics measurements of tetrahedrally-coordinated micellar block polymers, thus expanding the application of XPCS to ordered soft materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

The Auroral Oval Boundaries on January 10, 1997: A Comparison of Global Magnetospheric Simulations with UVI Images

We present the results of a global magnetospheric simulation of the initial period of the January 10 - 11, 1997 magnetic cloud event. Distinct magnetospheric boundaries derived from the model are mapped down to the ionosphere and compared to UVI images of the auroral oval from 1:00 to 4:30 UT. The convection reversal boundary, which coincides with the maximum Region I currents, tends to generally match the UVI equatorward boundary. This boundary is almost always poleward of the boundary separating Region 1 and Region 2 currents. The separatrix between open and closed magnetic field lines as mapped in the model matches the poleward boundary of the UVI images very well during quiet periods. During dynamic periods, however, when the separatrix can move several degrees in latitude in some sectors, the poleward boundary of the Region 1 currents matches the UVI images better.

Elsen, R.↗
DOE OSTI2023

Local structural ordering determines the mechanical damage tolerance of amorphous grain boundary complexions

Amorphous grain boundary complexions act as toughening features within a microstructure because they can absorb dislocations more efficiently than traditional grain boundaries. This toughening effect should be a strong function of the local internal structure of the complexion, which has recently been shown to be determined by grain boundary crystallography. To test this hypothesis, molecular dynamics are used here to simulate dislocation absorption and damage nucleation for complexions with different distributions of structural short-range order. Here, the complexion with a more disordered structure away from the dislocation absorption site is actually found to better resist crack nucleation, as damage tolerance requires delocalized deformation and the operation of shear-transformation zones through the complexion thickness. The more damage tolerant complexion accommodates plastic strain efficiently within the entire complexion, providing the key mechanistic insight that local patterning and asymmetry of structural short-range order controls the toughening effect of amorphous complexions.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Analytical studies of hypersonic viscous dissociated flows

This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Inger, George R.↗
NASA NTRS2004

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗
NASA NTRS2020

1- by 1-Foot Supersonic Wind Tunnel Sidewall Flowfield Compression Plate Study

A method for locally increasing the dynamic pressure in the test section of a supersonic wind tunnel is investigated from a theoretical, computational and experimental perspective, culminating in a reduced-scale, proof-of-concept test at Mach 3.44. Symmetric, variable-angle compression plates along the test section sidewalls generate oblique shockwaves to elevate the test section dynamic pressure (q). The experimental effort utilizes boundary-layer pitot rakes to sample the incoming flow near the floor, pitot and static rakes to measure the high-q core flow and numerous surface taps to map out the static pressure distribution along the floor, ceiling and sidewalls. A reliable and robust starting procedure is developed to establish supersonic flow in the test section and behind the compression plates. A characterization of the incoming floor boundary layer reveals previously documented features (such as a bulge toward the centerline) and an unexpected (and as-yet-unexplained) thickening of the boundary layer with increasing Reynolds number around total pressures of 40–60 psia. The high-q core flow is characterized for eight compression plate angles (4.0–8.0°), eight tunnel total pressures (20–139 psia) and three axial stations (24.95–27.35 inches). The peak dynamic pressure rise exceeds expectations. At conditions relevant to full-scale testing (6.0° and 139 psia), the dynamic pressure increases by a factor of 1.73, the Mach number becomes 2.85 and the total pressure drops imperceptibly. Basic distortion metrics show transverse non-uniformity of less than 5 %. Limited spatial resolution in the static rake measurements warrants additional, higher fidelity testing or computational analysis. Surface-mounted, fast-response instrumentation suggests minimal unsteadiness in the flow, less than a few percent of the mean (comparable to the sensor accuracy).

Wind tunnel testing, test techniques/methods, dyna↗
NASA NTRS1988

Thermospheric dynamics during November 21-22, 1981 - Dynamics Explorer measurements and thermospheric general circulation model predictions

Time-dependent aurora and magnetospheric convection parameterizations have been derived from solar wind and aurora particle data for November 21-22, 1981, and are used to drive the auroral and magnetospheric convection models that are embedded in the National Center for Atmospheric Research thermospheric general circulation model (TGCM). Neutral wind speeds and transition boundaries between the midlatitude solar-driven circulation and the high-latitude magnetospheric convection-driven circulation are examined on an orbit-by-orbit basis. The results show that TGCM-calculated winds and reversal boundary locations are in generally good agreement with Dynamics Explorer 2 measurements for the orbits studied. This suggests that, at least for this particular period of relatively moderate geomagnetic activity, the TGCM parameterizations on the eveningside of the auroral oval and polar cap are adequate.

Roble, R. G.↗