Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Numerical Modeling”

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.

974 records · Page 3

Surface Tension Driven Convection Experiment (STDCE)

Results are reported of the Surface Tension Driven Convection Experiment (STDCE) aboard the USML-1 (first United States Microgravity Laboratory) Spacelab which was launched on June 25, 1992. In the experiment 10 cSt silicone oil was placed in an open circular container which was 10 cm wide by 5 cm deep. The fluid was heated either by a cylindrical heater (1.11 cm dia.) located along the container centerline or by a CO 2 laser beam to induce thermocapillary flow. The flow field was studied by flow visualization. Several thermistor probes were placed in the fluid to measure the temperature distribution. The temperature distribution along the liquid free surface was measured by an infrared imager. Tests were conducted over a range of heating powers, laser beam diameters, and free surface shapes. In conjunction with the experiments an extensive numerical modeling of the flow was conducted. In this paper some results of the velocity and temperature measurements with flat and curved free surfaces are presented and they are shown to agree well with the numerical predictions.

S Ostrach

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials

Measure sea level air pressure from space to improve knowledge and forecasting of the atmospheric state

Modern numerical weather prediction (NWP) and analysis models require globally-observed meteorological data including sea-level pressure (SLP) for accurate operations. Up until now, SLP has only been measured by in situ instruments from ships, buoys, and ocean platforms. These measurements are sparse with large gaps, leaving models starved of this critical information to constrain the atmospheric state. Recent advancements in differential absorption radar (DAR) provide a path to close this critical observation gap through spaceborne observations in the coming decade, improving the analysis models relied upon for atmospheric research and the weather forecasts depended upon daily for public safety and commerce.

Matthew L Walker McLinden

Behaviors of Lunar Regolith Simulant Under Varying Gravitational Conditions

Understanding the behavior of regolith in varying gravity conditions is critical for space exploration and future missions. In this work, the gravity-driven hopper flow of lunar regolith simulant in different gravitational accelerations (terrestrial, lunar) is first observed experimentally. Numerical simulations (DEM) are then developed to understand the role which cohesive inter-particle forces play in such gravity-driven flow, using the theoretical framework of granular Bond number. Qualitative comparison between a terrestrial experiment and numerical simulation validated this framework. Following that, we numerically studied the dynamic behavior under varying gravitational conditions (from terrestrial to lunar to asteroid gravitational accelerations). We find that this behavior is extremely sensitive to the interplay of the gravity conditions and the attractive/cohesive forces among particles. The numerical and experimental results show that the complex interaction of these forces can drastically change the dynamics of the material producing effects relevant for variable gravity applications.

Soft Matter

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters

Expansion of Check-Cases for 6DOF Simulation

This is the Appendix containing a description of the solution for Case 1 in the assessment, “Expansion of Check-Cases for 6DOF Simulation”. For cases of spherical gravity, it is possible to provide a two-body solution without recourse to numerical integration and thus it is accurate to machine precision. Python code for a Keplerian Propagator (propagate.py) which produced a reference trajectory for Case 1 is provided in this appendix. There is also code for generating test cases which was used as an independent verification of the propagator. This is a high-level description of the algorithm employed. The documentation of each function includes implementation details, including equations for each task.

Modeling

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff

1998 NASA Seal/Secondary Air System Workshop

The 1998 NASA Seal/Secondary Air System Workshop was divided into three major areas with limited materials presented in Volume II: (1) overviews of the (NASAs high speed research (HSR) and DOE's advanced turbine engine systems (ATS)) gas turbine programs and the general aviation program (GAP) with emphasis on sealing methods and results; (2) sealing concepts and methods and results including experimental facilities and numerical predictions; and (3) reviews of the numerical engine simulation and aerospace vehicles and concepts (Trailblazer, Bantam, and X-38).

Bruce Steinetz

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig

Improvement in the Thermal-to-Structural Model Mapping Process for Integrated Modeling for the Roman Space Telescope

Integrated Modeling has been a key component of verifying optical requirements for the Nancy Grace Roman Space Telescope (RST) that are either impossible or impractical to verify exclusively through ground testing. Two major areas for integrated Modeling are Jitter and Thermal Distortion that require the exchanges of model performance predictions across disciplines. In both cases, distortions are impressed on optical models to evaluate the impact on boresight alignment and wave front error. In the case of Jitter, the disturbances are driven by reactions to motions most often from actuators; however, in the case of thermal distortion, the motions are driven by thermal expansion or contraction as a result of changing temperatures. This then requires a link further upstream to the thermal model, which is used to predict the thermal performance and temperature gradients and stability. The process for mapping temperatures from a thermal model to a corresponding structural model has been performed numerous times through the RST project lifecycle, with improvements in the accuracy, verification, and effort sought throughout. This paper describes some of the recent improvements to the process, including: capture of the visualization parameters, automatic generation of the mapped images for both the thermal and structural model groupings, and reduction in the effort to assemble the full set of mapped temperatures. These upgrades have greatly reduced the manual effort associated with thermal mapping and allowed for faster turn-around of Integrated Modeling predictions.

Thermal Mapping

Impact Testing of 410 Stainless Steel for Material Impact Model Development

A project is underway to develop a consistent set of material properties, impact test data, and failure analysis for a variety of metallic aircraft materials that can be used to develop improved impact failure and deformation models. This project is jointly funded by the NASA Glenn Research Center and the Federal Aviation Administration William J. Hughes Technical Center. Particular features of this set of data are that all material property and impact test data are obtained using traceable material, the test methods and procedures are extensively documented, and all the raw data are available. Four parallel efforts are currently underway: measurement of material deformation and failure response over a wide range of strain rates and temperatures, failure analysis of material property specimens and impact test articles, development of improved numerical modeling techniques for deformation and failure, and impact testing of flat panels and substructures for model validation. This report describes impact testing performed on 410 stainless steel sheet and plate samples of different thicknesses with two different types of projectiles, one a regular cylinder and one with a more complex geometry incorporating features representative of a jet engine fan blade. Data from this testing will be used in validating material models developed under this program. The material tests and the material models developed in this program will be published in separate reports.

410 Stainless Steel

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi

In-Space Demonstration of Spray-on Solar White Paint

Discs coated with spray-on solar white and two similar thermal control coatings were flown in low earth orbit for roughly five hours. While in orbit, the coated discs were exposed to direct sunlight and to deep space for short intervals. Measurements of disc temperatures showed that the spray-on solar white coating had the lowest absorption of sunlight. Numerical modeling was used to infer that spray-on solar white had a solar absorptance between 0.02 and 0.04 and a thermal emittance between 0.93 and 0.99, both of which are consistent with terrestrial laboratory measurements.

Solar White

Active learning using hybrid surrogate tool life modeling for machining process optimization

Here, this paper describes an active learning approach for part-to-part iterative machining process optimization using a hybrid surrogate tool life model. A probabilistic interpolating tool life model is developed by combining the empirical Taylor-type tool life equation and the model fit error. The probabilistic tool life model is then used to calculate the machining cost per part distribution. The optimal machining parameters are selected using an expected improvement in machining cost per part criterion. The method is validated numerically using experimental results; the results show a median convergence error of 2.2% after three tests over 400 simulations. The method is validated experimentally on two industrial applications for Ti-6Al-4V roughing resulting in a cost per part reduction greater than 23% after two tests. The described method is a robust solution for rapid convergence to optimal machining parameters in an industrial production environment.

Active learning

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics