Results for “Structural modeling and simulation” — Page 34

NASA NTRS2011

Scientific Visualization to Study Flux Transfer Events at the Community Coordinated Modeling Center

In this paper we present results of modeling of reconnection at the dayside magnetopause with subsequent development of flux transfer event signatures. The tools used include new methods that have been added to the suite of visualization methods that are used at the Community Coordinated Modeling Center (CCMC). Flux transfer events result from localized reconnection that connect magnetosheath magnetic field and plasma with magnetospheric fields and plasma and results in flux rope structures that span the dayside magnetopause. The onset of flux rope formation and the three-dimensional structure of flux ropes are studied as they have been modeled by high-resolution magnetohydrodynamic simulations of the dayside magnetosphere of the Earth. We show that flux transfer events are complex three-dimensional structures that require modern visualization and analysis techniques. Two suites of visualization methods are presented and we demonstrate the usefulness of those methods through the CCMC web site to the general science user.

Rastatter, Lutz↗

NASA NTRS2016

Effects of Zoledronate and Mechanical Loading during Simulated Weightlessness on Bone Structure and Mechanical Properties

Space flight modulates bone remodeling to favor bone resorption. Current countermeasures include an anti-resorptive drug class, bisphosphonates (BP), and high-force loading regimens. Does the combination of anti-resorptives and high-force exercise during weightlessness have negative effects on the mechanical and structural properties of bone? In this study, we implemented an integrated model to mimic mechanical strain of exercise via cyclical loading (CL) in mice treated with the BP Zoledronate (ZOL) combined with hindlimb unloading (HU). Our working hypothesis is that CL combined with ZOL in the HU model induces additive structural and mechanical changes. Thirty-two C57BL6 mice (male,16 weeks old, n8group) were exposed to 3 weeks of either HU or normal ambulation (NA). Cohorts of mice received one subcutaneous injection of ZOL (45gkg), or saline vehicle, prior to experiment. The right tibia was axially loaded in vivo, 60xday to 9N in compression, repeated 3xweek during HU. During the application of compression, secant stiffness (SEC), a linear estimate of slope of the force displacement curve from rest (0.5N) to max load (9.0N), was calculated for each cycle once per week. Ex vivo CT was conducted on all subjects. For ex vivo mechanical properties, non-CL left femurs underwent 3-point bending. In the proximal tibial metaphysis, HU decreased, CL increased, and ZOL increased the cancellous bone volume to total volume ratio by -26, +21, and +33, respectively. Similar trends held for trabecular thickness and number. Ex vivo left femur mechanical properties revealed HU decreased stiffness (-37),and ZOL mitigated the HU stiffness losses (+78). Data on the ex vivo Ultimate Force followed similar trends. After 3 weeks, HU decreased in vivo SEC (-16). The combination of CL+HU appeared additive in bone structure and mechanical properties. However, when HU + CL + ZOL were combined, ZOL had no additional effect (p0.05) on in vivo SEC. Structural data followed this trend with ZOL not modulating trabecular thickness in CL + NAHU mice. In summary, our integrated model simulates the combination of weightlessness, exercise-induced mechanical strain, and anti-resorptive treatment that astronauts experience during space missions. Based on these results, we conclude that, at the structural and stiffness level, zoledronate treatment during simulated spaceflight does not impede the skeletal response to axial compression. In contrast to our hypothesis, our data show that zoledronate confers no additional mechanical or structural benefit beyond those gained from cyclical loading.

Mechanical Loading↗

NASA NTRS1983

A two-equation turbulence model for two-phase flows

A two-equation turbulence model has been developed for predicting two-phase flows. The two equations describe the conservation of turbulence kinetic energy and dissipation rate of that energy for the carrier fluid in a two-phase flow. They have been derived rigorously from the momentum equations of the carrier fluid. Closure of the time-mean equations is achieved by modeling the turbulent correlations up to third order. The new model eliminates the need to simulate in an ad hoc manner the effects of the dispersed phase on turbulence structure. Preliminary testing indicates that the model is successful in predicting the main features of a round gaseous jet laden with uniform-size solid particles.

Elghobashi, S. E.↗

NASA NTRS1983

The N/Rev phenomenon in simulating a blade-element rotor system

When a simulation model produces frequencies that are beyond the bandwidth of a discrete implementation, anomalous frequencies appear within the bandwidth. Such is the case with blade element models of rotor systems, which are used in the real time, man in the loop simulation environment. Steady state, high frequency harmonics generated by these models, whether aliased or not, obscure piloted helicopter simulation responses. Since these harmonics are attenuated in actual rotorcraft (e.g., because of structural damping), a faithful environment representation for handling qualities purposes may be created from the original model by using certain filtering techniques, as outlined here. These include harmonic consideration, conventional filtering, and decontamination. The process of decontamination is of special interest because frequencies of importance to simulation operation are not attenuated, whereas superimposed aliased harmonics are.

Mcfarland, R. E.↗

NASA NTRS1983

The N/rev phenomenon in simulating a blade-element rotor system

When a simulation model produces frequencies that are beyond the bandwidth of a discrete implementation, anomalous frequencies appear within the bandwidth. Such is the case with blade element models of rotor systems, which are used in the real time, man in the loop simulation environment. Steady state, high frequency harmonics generated by these models, whether aliased or not, obscure piloted helicopter simulation responses. Since these harmonics are attenuated in actual rotorcraft (e.g., because of structural damping), a faithful environment representation for handling qualities purposes may be created from the original model by using certain filtering techniques, as outlined here. These include harmonic consideration, conventional filtering, and decontamination. The process of decontamination is of special interest because frequencies of importance to simulation operation are not attenuated, whereas superimposed aliased harmonics are.

Mcfarland, R. E.↗

NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

DOE OSTI2024

The Role of Unit-Cell Topology in Modulating the Compaction Response of Additively Manufactured Cellular Materials using Simulations and Validation Experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗

NASA NTRS2022

Impact of Meteorological Factors on the Mesoscale Morphology of Cloud Streets during a Cold-Air Outbreak over the Western North Atlantic

Postfrontal clouds (PFC) are ubiquitous in the marine boundary layer, and their morphology is essential to estimating the radiation budget in weather and climate models. Here we examine the roles of sea surface temperature (SST) and meteorological factors in controlling the mesoscale morphology and evolution of shallow clouds associated with a cold-air outbreak that occurred on 1 March 2020 during phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). Our results show that the simulated PFC structure and ambient conditions by the Weather Research and Forecasting (WRF) Model are generally consistent with observations from GOES-16 and dropsonde measurements. We also examine the thermodynamical and dynamical influences in the cloud mesoscale morphology using WRF sensitivity experiments driven by two meteorological forcing datasets with different domain-mean SST and spatial gradients, which lead to dissimilar values of hydrometeor water path and cloud core fraction. The SST from ERA5 leads to weaker stability and higher inversion height than the SST from FNL does. In addition, the use of large-scale meteorological forcings from ERA5 yields a distinctive time evolution of wind direction shear in the inner domain, which favors the formation and persistence of longer cloud rolls. Both factors contribute to a change in the time evolution of domain-mean water path and cloud core fraction of cloud streets. Our study takes advantage of the simulation driven by the differences between two large-scale forcing datasets to illustrate the importance of SST and wind direction shear in the cloud street morphology in a realistic scenario

Jingyi Chen↗

DOE OSTI2024

The role of unit cell topology in modulating the compaction response of additively manufactured cellular materials using simulations and validation experiments

Additive manufacturing has enabled a transformational ability to create cellular structures (or foams) with tailored topology. Compared to their monolithic polymer counterparts, cellular structures are potentially suitable for systems requiring materials with high specific energy-absorbing capability to provide enhanced damping. In this work, we demonstrate the utility of controlling unit-cell topology with the intent of obtaining a desired stress–strain response and energy density. Using mesoscale simulations that resolve the unit-cell sub-structures, we validate the role of unit-cell topology in selectively activating a buckling mode and thereby modulating the characteristic stress–strain response. Simulations incorporate a linear viscoelastic constitutive model and a hyperelastic model for simulating large deformation of the polymer under both tension and compression. Simulated results for nine different cellular structures are compared with experimental data to gain insights into three different modes of buckling and the corresponding stress–strain response.

36 MATERIALS SCIENCE↗

NASA NTRS1986

Space Station dynamic modeling, control and simulation

Dynamic disturbances of many orders of magnitude greater than those of conventional spacecraft will be routine for the Space Station. Accurate knowledge of inflight structural dynamics and disturbances will be lacking. System identification will reduce uncertainties. To deal with the remaining model errors and time-varying elements, adaptive control may be required. Dynamic models for two Space Station configurations are derived. A direct model reference adaptive control algorithm is synthesized and evaluated with respect to the station model parameter errors and plant dynamics truncations. Both the rigid body and the flexible modes are treated. Simulation results show that convergence of the adaptive algorithm can be achieved in 100 to about 300 seconds with reasonable performance even during Shuttle hard docking operations in which station mass and inertia are instantaneously changed by more than 100 percent.

Ih, C.-H. C.↗

NASA NTRS1987

Advances and trends in computational structural mechanics

The development status and applicational range of techniques in computational structural mechanics (CSM) are evaluated with a view to advances in computational models for material behavior, discrete-element technology, quality assessment, the control of numerical simulations of structural response, hybrid analysis techniques, techniques for large-scale optimization, and the impact of new computing systems on CSM. Primary pacers of CSM development encompass prediction and analysis of novel materials for structural components, computational strategies for large-scale structural calculations, and the assessment of response prediction reliability together with its adaptive improvement.

Noor, Ahmed K.↗

NASA NTRS2012

Model-Based Verification and Validation of Spacecraft Avionics

Verification and Validation (V&V) at JPL is traditionally performed on flight or flight-like hardware running flight software. For some time, the complexity of avionics has increased exponentially while the time allocated for system integration and associated V&V testing has remained fixed. There is an increasing need to perform comprehensive system level V&V using modeling and simulation, and to use scarce hardware testing time to validate models; the norm for thermal and structural V&V for some time. Our approach extends model-based V&V to electronics and software through functional and structural models implemented in SysML. We develop component models of electronics and software that are validated by comparison with test results from actual equipment. The models are then simulated enabling a more complete set of test cases than possible on flight hardware. SysML simulations provide access and control of internal nodes that may not be available in physical systems. This is particularly helpful in testing fault protection behaviors when injecting faults is either not possible or potentially damaging to the hardware. We can also model both hardware and software behaviors in SysML, which allows us to simulate hardware and software interactions. With an integrated model and simulation capability we can evaluate the hardware and software interactions and identify problems sooner. The primary missing piece is validating SysML model correctness against hardware; this experiment demonstrated such an approach is possible.

MBV&V↗

NASA NTRS1991

Computation of rotor aerodynamic loads with a constant vorticity contour free wake model

An analytical method is presented which facilitates the study of isolated rotors with an improved approach to wake simulation. Vortex filaments are simulated along contours of constant sheet strength for the sheet of vorticity resulting from each rotor blade. Curved vortex elements comprise the filaments which can be distorted by the local velocity field. Called the Constant Vorticity Contour wake model, the approach permits the simulation of the blades' wakes corresponding to the full span of the rotor blade. The discretization of the wake of the rotor blade produces spacing and structure that are consistent with the spatial and temporal variations in the loading. A vortex-lattice aerodynamic model of the blade is also included which introduces a finite-element structural model of the blade and consideration of the force and moment trim analysis. Results of the present version of the simulation, called RotorCRAFT, are found to correlate well with H-34 flight-test data.

Quackenbush, Todd R.↗

NASA NTRS1989

Six-degree-of-freedom aircraft simulation with mixed-data structure using the applied dynamics simulation language, ADSIM

A realistic simulation of an aircraft in the flight using the AD 100 digital computer is presented. The implementation of three model features is specifically discussed: (1) a large aerodynamic data base (130,00 function values) which is evaluated using function interpolation to obtain the aerodynamic coefficients; (2) an option to trim the aircraft in longitudinal flight; and (3) a flight control system which includes a digital controller. Since the model includes a digital controller the simulation implements not only continuous time equations but also discrete time equations, thus the model has a mixed-data structure.

Savaglio, Clare↗

NASA NTRS1981

Electromagnetic deep-probing (100-1000 KMS) of the Earth's interior from artificial satellites: Constraints on the regional emplacement of crustal resources

An algorithm was developed to address the problem of electromagnetic coupling of ionospheric current systems to both a homogeneous Earth having finite conductivity, and to an Earth having gross lateral variations in its conductivity structure, e.g., the ocean-land interface. Typical results from the model simulation for ionospheric currents flowing parallel to a representative geologic discontinuity are shown. Although the total magnetic field component at the satellite altitude is an order of magnitude smaller than at the Earth's surface (because of cancellation effects from the source current), the anomalous behavior of the satellite observations as the vehicle passes over the geologic contact is relatively more important pronounced. The results discriminate among gross lithospheric structures because of difference in electrical conductivity.

Hermance, J. F.↗

DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗

NASA NTRS2018

Molecular Dynamics Simulations of Liquid and Polymer Electrolytes for Energy Storage Devices

Advancing beyond current lithium-ion technology is necessary in order to enable energy storage devices for electric airplanes. Electrolyte stability is a key limiting factor, yet the design of improved electrolytes remains a formidable challenge. Molecular dynamics (MD) simulations are a powerful tool for studying electrolytes, since they can be used to evaluate structural, thermodynamic, and transport properties, and can provide molecular-level detail often inaccessible to experimental techniques. Our computational materials groups at the NASA Ames Research Center has developed models and methods to accurately simulate both liquid and polymer electrolytes.We report the results from atomistic MD simulations of several electrolyte materials, with lithium salts dissolved in ionic liquids, dimethoxyethane (DME), and polyethylene oxide (PEO). For improved accuracy, we employ polarizable models, where each atom is given an environment-dependent atomic dipole. The simulations accurately predict bulk transport properties, including viscosity, diffusion, and ionic conductivity, in quantitative agreement with available experimental data. Moreover, the simulations provide important insights into the solvation structure of the lithium ions.We also report the results from coarse-grained MD simulations of polyanion electrolytes. In order to more efficiently capture the longer length- and time-scales of these systems, we employ a generic bead-spring model. These simulations provide important insight into how the polymer chain architecture and ionic interaction strengths affect the ionic aggregation behavior and cation dynamics. Despite the simplicity of the model, the simulations yield qualitative agreement with experimental data for similar systems.

Abbott, Lauren J.↗