On the Development Of Preconditioning Methods for High-Order Multiphysics Problems
Explore the source record for details and available documents.
STMD↗
SEARCH · Search NASA
Results for “Multiphysics methods”
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.
Oxidation of Silicon Particles Suspended in Mullite Matrix With Dual-Mode Vacancy Diffusion of Oxygen
Mullite bond coats filled with silicon particles have recently been proposed for the next generation of environmental barrier coatings (EBCs). The sacrificial oxidation of the silicon particles eliminates some of the environmental durability issues associated with traditional EBC designs. To investigate this novel design approach, a mathematical formulation and resulting numerical method have been developed for simulating the oxidation of silicon particles embedded in a mullite matrix. The focus is on the oxidation of the silicon in a dry oxygen environment where the oxidant is presumably oxygen. Oxygen transport from the exterior environment through the mullite is assumed to occur via two lattice defect diffusion mechanisms. As a result, the rate of oxygen transport is dictated by the self-diffusivity of the lattice defects. The numerical method is implemented using the COMSOL Multiphysics® program (COMSOL, Inc.). The method is applied to simulate the oxidation of silicon particles in a mullite pellet in recent oxidation experiments. Numerical results compare favorably with the oxidation behavior observed in the oxidation experiments. The numerical method is also used to investigate the effects of various features of the mullite-silicon particle system on the oxidation behavior.
silicon oxidation↗
Pragmatic Stress Prediction on Additively Manufactured Coupons
Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.
Additive Manufacturing↗
Pragmatic Stress Prediction on Additively Manufactured Coupons
Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.
pragmatic method↗
The Role of Multiphysics Simulation in Multidisciplinary Analysis
This article describes the applications of the Spectrum(Tm) Solver in Multidisciplinary Analysis (MDA). Spectrum, a multiphysics simulation software based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling as well as the interaction between these disciplines. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena. Interaction constraints are enforced in a fully-coupled manner using the augmented-Lagrangian method. Within the multiphysics framework, the finite element treatment of fluids is based on Galerkin-Least-Squares (GLS) method with discontinuity capturing operators. The arbitrary-Lagrangian-Eulerian method is utilized to account for deformable fluid domains. The finite element treatment of solids and structures is based on the Hu-Washizu variational principle. The multiphysics architecture lends itself naturally to high-performance parallel computing. Aeroelastic, propulsion, thermal management and manufacturing applications are presented.
Rifai, Steven M.↗
Extracting Damping Ratio from Dynamic Data and Numerical Solutions
There are many ways to extract damping parameters from data or models. This Technical Memorandum provides a quick reference for some of the more common approaches used in dynamics analysis. Described are six methods of extracting damping from data: the half-power method, logarithmic decrement (decay rate) method, an autocorrelation/power spectral density fitting method, a frequency response fitting method, a random decrement fitting method, and a newly developed half-quadratic gain method. Additionally, state-space models and finite element method modeling tools, such as COMSOL Multiphysics (COMSOL), provide a theoretical damping via complex frequency. Each method has its advantages which are briefly noted. There are also likely many other advanced techniques in extracting damping within the operational modal analysis discipline, where an input excitation is unknown; however, these approaches discussed here are objective, direct, and can be implemented in a consistent manner.
Casiano, M. J.↗
Computationally Guided Development of Components for High Energy Density Solid-State Lithium-Sulfur Batteries
All electric vertical take-off and landing vehicles (eVTOL) for urban air mobility (UAM) concepts face numerous challenging technical barriers before their introduction into the consumer marketplace. The most challenging of these technical barriers to overcome is developing an energy storage system capable of meeting the rigorous aerospace safety and performance criteria1. The performance metrics for eVTOL craft, such as specific energy, specific power, and safety, exceed those of electric automobiles by a factor of two to four. Current state-of-the-art (SOA) lithium-ion batteries are incapable of meeting the key performance criteria of energy and safety for eVTOL. Therefore, next generation advanced chemistries and designs must be developed to meet required performance metrics for electric aviation2. Beyond lithium-ion chemistries, such as lithium-sulfur, show promise in their high energy, while limitations exist in their power and cyclability due to low electrical conductivity and high intermediate solubility in organic liquid electrolytes. Several strategies to overcome the low electrical conductivity involve the use of selenium as a dopant in the active sulfur material, along with the incorporation of 2-dimensional electron-conducting holey-graphene to improve the composite cathodes electronic conductivity. Furthermore, combining this chemistry with a solid-electrolyte avoids the components’ dissolution issues3. However, the development of composite solid-state cathodes is non-trivial as several components must be intimately mixed so that the active component has sufficient access to both electrons and lithium ions to undergo full electrochemical conversion. Mathematical modeling of battery components can assist experimental design through a robust and rigorous combination of computational modeling techniques covering multiple length scales. The objective is to leverage modern computational materials methods combined with battery multiphysics tools to develop radically advanced compatible cathode and electrolyte materials, build and test solid state lithium-sulfur cells and packs. A NASA-based cross-organizational team of high-powered experts combined integrated computational predictive modeling, fundamental chemistry analysis, advanced material science, and battery cell development to tackle this very challenging, multidisciplinary problem. This presentation will show a multiscale computational modeling approach that has produced a novel particle dynamics method called Solid Electrolyte Sphere Approximation Model (SESAM). SESAM modeling targets the 1-10 µm scale structures and provides electromechanical and grain interactions for predictive design guidelines for the manufacturing of solid-state components. Parameters such as particle size and volume fraction of the constituent materials were modeled and experimentally fabricated to optimize electrochemical performance through improved microstructure design. Experimental feedback was provided through ionic and electronic conductivity assessment and structural analysis of developed materials and cell components.
battery↗
Integrated Computational-Experimental Development of Lithium-Air Batteries for Electric Aircraft
The primary obstacle to enable NASA's vision of Green Aviation is the extraordinary energy storage requirements for electric aircraft. Significant advances in high energy, rechargeable, safe batteries are required to enable electric aviation. Boeing's SUGAR and NASA studies have identified 400 Wh/kg as the threshold energy density for general aviation and 750 Wh/kg for commercial regional air service. State of the Art Lithium Ion Battery (LIB) technology currently has a density of 200 Wh/kg and is expected to plateau at 300 Wh/kg due to fundamental chemistry limitations making it unsuitable for future electric aircraft. Additional demanding requirements include high power, rechargeability, and high safety. Such battery technology does not currently exist. The recent considerable activity in battery research (DOE, Tesla Gigafactory, etc) overwhelmingly has been geared towards reducing cost and improving safety of LIB technology in order to promote the adoption of electric automobiles; and thus it is expected to have little impact on electric aviation development. New battery materials will be needed for the "Beyond Li Ion" (BLI) technologies required for high energy, safe electric aviation. Li-Air batteries have the highest known theoretical energy density (3400 Wh/kg) and therefore and if realized promises to transform the global transportation system. These high energy batteries have the potential to meet the energy storage challenges of current and future NASA aeronautics and space missions in addition to many terrestrial transportation applications as well. However, this technology requires significant components development and integration, as it is currently unable to achieve aircraft requirements. The objective of this project is to leverage modern computational materials methods combined with battery multiphysics tools to develop radically advanced compatible cathode and electrolyte materials, build several Li-Air cells, and flight-demonstrate the corresponding Li-Air battery packs. A significant problem for current Lithium-Air batteries is large scale decomposition of the battery electrolyte during operation leading to battery failure after a handful of charge/discharge cycles. Therefore, development of large scale, ultra-high energy, rechargeable, and safe Lithium-Air batteries require highly stable electrolytes that are resistant to decomposition under operating conditions. A NASA-based cross-organizational "dream team" of high-powered experts combined integrated supercomputer modeling, fundamental chemistry analysis, advanced material science, and battery cell development to tackle this very challenging, multidisciplinary problem. The ultimate goal for the team is to develop an integrated experimental/computational infrastructure to produce a reliable predictive capability for the selection of optimal components, their fabrication parameters, and "design rules" of novel cell components for advanced ultra-high energy batteries that can meet energy storage challenges of NASA missions and many terrestrial transportation applications.
Li-air battery↗
Coupled Radiative Thermal and Nonlinear Stress Analysis for Thermal Deformation in Large Space Structures
Large space structures are capable of large thermal deformations in the space environment. A case of large-scale thermal deformation was observed in the analysis of the Near Earth Asteroid Scout solar sail, with predicted tip displacements of more than one meter in seven-meter booms. Experimental data supports the broad conclusions of the analysis, but shows poor agreement on the details of the thermal deformation. Prediction that is precise enough to drive engineering decisions will require coupled thermal-stress analysis with features that are not found in current multiphysics codes. This paper describes a simple method for stepwise coupling between commercial nonlinear stress analysis software and radiative thermal analysis software. Results are presented for a round stainless steel tube, which is a common case in existing literature.
Stohlman, Olive R.↗
Multiphysics Simulation of Active Hypersonic Lip Cooling
This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.
Melis, Matthew E.↗
Comparing Acoustic Prediction Methods for Additively Manufactured Porous Strutures
While macroscale methods for predicting the acoustic properties of porous structures have been popular in the past, they often require time-consuming manufacturing and testing workflows. Meanwhile, microscale approaches allow the prediction of transport parameters based exclusively on a periodic structure’s unit cell geometry. Here, we compare these methods to predict the characteristic impedance of additively manufactured porous structures. We use the microscale approach to estimate the geometry’s transport parameters, then predict the characteristic properties using the Johnson-Champoux-Allard (JCA) model. We measure the acoustic properties of the printed structures using a normal incidence impedance tube and estimate the transport parameters using an inverse characterization approach. We use the two-thickness method as a macroscale approach to predict the characteristic properties from the measured surface impedances of two sample thicknesses. Finally, we compare these characteristic prediction methods. Our results show that the inverse characterization and two-thickness methods offer the closest match to the measured values at low frequencies.
impedance↗
Analytical Design and Performance Estimation Methods for Aircraft Permanent Magnet Synchronous Machines
The design of an electric motor drivetrain is a complex multiphysics problem. Low fidelity motor drivetrain sizing can be a key tool in the design cycle of an electric motor drivetrain and for system level studies of aircraft configurations. However, low fidelity sizing can lead to misleading results if all the physics involved in a motor design are not properly accounted for. This paper provides details on modeling approaches for initial design and sizing of permanent magnet synchronous electric machines. The goal of this paper is to provide the reader an understanding of the key principles of motor design and some modeling approaches to perform initial sizing of an electric motor and its inverter.
Thomas F. Tallerico↗
Piezoelectric Vibration Damping Study for Rotating Composite Fan Blades
Resonant vibrations of aircraft engine blades cause blade fatigue problems in engines, which can lead to thicker and aerodynamically lower performing blade designs, increasing engine weight, fuel burn, and maintenance costs. In order to mitigate undesirable blade vibration levels, active piezoelectric vibration control has been investigated, potentially enabling thinner blade designs for higher performing blades and minimizing blade fatigue problems. While the piezoelectric damping idea has been investigated by other researchers over the years, very little study has been done including rotational effects. The present study attempts to fill this void. The particular objectives of this study were: (a) to develop and analyze a multiphysics piezoelectric finite element composite blade model for harmonic forced vibration response analysis coupled with a tuned RLC circuit for rotating engine blade conditions, (b) to validate a numerical model with experimental test data, and (c) to achieve a cost-effective numerical modeling capability which enables simulation of rotating blades within the NASA Glenn Research Center (GRC) Dynamic Spin Rig Facility. A numerical and experimental study for rotating piezoelectric composite subscale fan blades was performed. It was also proved that the proposed numerical method is feasible and effective when applied to the rotating blade base excitation model. The experimental test and multiphysics finite element modeling technique described in this paper show that piezoelectric vibration damping can significantly reduce vibrations of aircraft engine composite fan blades.
Min, James B.↗
Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials
The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.
Liou, Frank↗
Microwave Processing of Planetary Surfaces for the Extraction of Volatiles
In-Situ Resource Utilization will be necessary for sustained exploration of space. Volatiles are present in planetary soils, but water by far has the most potential for effective utilization. The presence of water at the lunar poles, Mars, and possibly on Phobos opens the possibility of producing LOX for propellant. Water is also a useful radiation shielding material , and valuable to replenish expendables (water and oxygen) required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating soil effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on soil dielectric properties. Methods for complex electric permittivity and magnetic permeability measurement are being developed and used for measurements of lunar soil simulants. A new method for delivery of microwaves deep into a planetary surface is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. We are planning to set up a planetary testbed in a large vacuum chamber in the coming year. Recent results are discussed.
Ethridge, Edwin C.↗
Microwave Processing of Planetary Surfaces for Volatile Extraction
In-Situ Resource Utilization will be necessary for sustained exploration of space. Volatiles are present in planetary soils, but water by far has the strongest potential for effective utilization. The presence of water at the lunar poles, Mars, and possibly on Phobos opens the possibility of producing LOX for propellant. Water is also a useful radiation shielding material and water (and oxygen) are expendables that are also required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating soil effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on soil dielectric properties. Methods for measuring the complex electric permittivity and magnetic permeability are being developed and have been measured for some lunar soil simulants at 0.5, 2.45, and 10 GHz from room temperature down to liquid nitrogen temperature. A new method for delivery of microwaves deep into a planetary surface is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. We have plans to set up a planetary testbed in a large vacuum chamber in the coming year. Recent results will be presented.
Ethridge, Edwin C.↗
Finite Element Analysis of Three Methods for Microwave Heating of Planetary Surfaces
In-Situ Resource Utilization will be Ground Breaking technology for sustained exploration of space. Volatiles are present in planetary regolith, but water by far has the most potential for effective utilization. The presence of water at the lunar poles and Mars opens the possibility of using the hydrogen for propellant on missions beyond Earth orbit. Likewise, the oxygen could be used for in-space propulsion for lunar ascent/descent and for space tugs from low lunar orbit to low Earth orbit. Water is also an effective radiation shielding material as well as a valuable expendable (water and oxygen) required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating regolith effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within, much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on regolith dielectric properties. New methods for delivery of microwaves into lunar and planetary surfaces is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. Recent results are discussed.
Ethridge, Edwin↗
Analysis of Arc Welding Process in Space
This work is motivated by NASA plans to conduct welding experiments on the ISS. It is expected that deployment of welding and additive manufacturing technologies in the space environment has the potential to revolutionise how orbiting platforms are designed, manufactured, and assembled. However, the structure, composition and quality of a weld is extremely dependent on the environment and can be difficult to control in space. Shielding gases would also be tough to manage as gases behave differently in zero gravity and airless environments. Additional points of concern are related to the spatter and sparks dynamics in space. Therefore, there is a need for a more basic understanding of welding processes by computational modelling. To provide such an insight we developed state of art models of the ARC torch, droplet detachment and transfer, and the melt pool build up using magnetohydrodynamics approximation and level set method for two-phase liquid metal/-gas flow modelling. Two finite element models were built using 2D axisymmetric geometry in COMSOL Multiphysics®: (i) stationary model of the ARC torch and (ii) dynamical model of droplet detachment and transfer. Both models demonstrate reasonable agreement with earlier experimental observations and high sensitivity to the temperature dependence of the thermophysical parameters on the system materials. The models were used to provide physical insight into ARC welding in various environments and geometries. It was shown, in particular, that there is a significant probability of gas bubbles trapping in the meltpool and a possibility of unbounded wondering of sparks in the welding chamber in zero gravity. In addition, we estimated metal evaporation rate that may be hazardous in the confined environment of the ISS. These issues raise concern of quality of the weld and safety for ISS applications.
ARC welding↗