Lower Length Scale Informed Continuum Modeling of Femtosecond Laser Processing
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Despite the availability of computational resources and advancements in numerical computing capabilities, the multiscale models core to understanding, predicting the behaviors of, and designing energy and environmental systems involving porous media are still 1.) developed through by-hand derivations and 2.) limited by many methodological assumptions employed during model derivation. As a result, the advancement of effective media models for engineering DOE mission-critical systems (e.g., batteries, flow batteries, electrolyzers, geothermal systems, subsurface chemical storage systems, etc.) is slow (i.e., it takes years for models to traverse from stages of “development” to “practical utilization”), hindering our ability to effectively optimize such systems and stay at the cutting-edge of the energy frontier. In this work, we aimed to address these limitations by 1.) automating and accelerating multiscale model derivation via symbolic computing and 2.) develop a novel multiscale modeling methodology for flow and transport through porous media that avoids the typical assumptions hindering previous models. As a result of our efforts, we 1.) developed a hybrid symbolic-numeric code called Fouriera for fully-automating the implementation of multiphysical and phase-field models via the Fourier spectral method for materials science research, and 2.) advanced a multiscale modeling methodology called The Method of Finite Averages that rigorously predicts the behaviors of flow and transport through heterogeneous porous media under the influence of non-local effects and strong advection. Ultimately, these deliverables provide strong foundations from which further efforts can advance multiscale modeling tools and capabilities that do not intrinsically rely on 1.) the speed and mathematical capabilities of humans, nor 2.) the methodological assumptions limiting current models.
The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.
A mathematical and computational analysis capability has been developed for calculating the effective mechanical properties of three-dimensional periodic truss-like structures. Two models are studied in detail. The first, called the octetruss model, is a three-dimensional extension of a two-dimensional model, and the second is a cubic model. Symmetry considerations are employed as a first step to show that the specific octetruss model has four independent constants and that the cubic model has two. The actual values of these constants are determined by averaging the contributions of each rod element to the overall structure stiffness. The individual rod member contribution to the overall stiffness is obtained by a three-dimensional coordinate transformation. The analysis shows that the effective three-dimensional elastic properties of both models are relatively close to each other.
An approach for modeling dynamic equations of motion of a plate attached with rigid bodies is presented. The equations of motion are developed using the principle of virtual work. Lagrange multipliers are used as interaction forces and/or moments to maintain prescribed constraints which is the basis of the interconnection between the plate and rigid bodies. The overall approach is unique in the sense that a continuous model described by a family of partial differential equations is established. An approximate formulation by using variational method is established yielding a solution compatible with the assumed degree of approximation. The formulation is useful particularly when parametric study of dynamic response for a satellite power system is desired. As an example, an approximate governing equation of algebraic eigenvalue problem is given for a dual microwave power transmission system. Controller design is discussed.
A methodology is proposed for the integrated design of structures and their controls. The integration of design is accomplished through simultaneous selection, via optimization, of designer-specified parameters governing the structure and its controller. The cost function is chosen to address the primary structural and control goals. Inequality constraints are added to insure that additional design requirements are met. The methodology is illustrated through development of an example involving a cantilever beam. A distributed transfer function is derived from the Bernoulli-Euler beam equation. A controller parameterization is then chosen, and performance requirements specified. A controller is then obtained via optimization, which is shown to give good performance.
A new method for obtaining the maximum frequency design of a beam-like repetitive lattice structure is presented. Using existing techniques, the lattice is first modeled as an equivalent anisotropic Timoshenko beam. The computation of the stiffness and inertial properties of the beam, determined by matching the strain and kinetic energies of the beam with those of the lattice, is facilitated by the repetitive nature of the lattice. The optimum design is obtained by maximizing Rayleigh's quotient using methods of variational calculus. For the problem selected, results show excellent agreement with those obtained by traditional finite-element methods. Moreover, unlike FE methods, cpu time is relatively unaffected by the size of the truss.
Results of a direct simulation Monte Carlo method for a hypersonic flow about a flat plate at a 40 deg angle of attack have been compared with corresponding results from a theory for fully viscous shock layers (FVSLs). Using the 13-moment equations for a Maxwell gas, it is demonstrated that nonequilibrium and equilibrium FVSL flows can be correlated. With the exception of the pressure density, most of the flow properties along a streamline (including heat flux and shear and normal stresses) are correctly predicted to leading order by the Navier-Stokes model.
A two dimensional hydrocode based on 2-D SALE was modified to include strength effects and fragmentation equations for fracture resulting from tensile stress in one dimension. Output from this code includes a complete fragmentation summary for each cell of the modeled object: fragment size (mass) distribution, vector velocities of particles, peak values of pressure and tensile stress, and peak strain rates associated with fragmentation. Contour plots showing pressure and temperature at given times within the object are also produced. By invoking axial symmetry, three dimensional events can be modeled such as zero impact parameter collisions between asteroids. The code was tested against the one dimensional model and the analytical solution for a linearly increasing tensile stress under constant strain rate.
A continuous beam model is developed for approximate analysis of a large, slender, beam-like truss. The model is incorporated in a multi-level optimization scheme for the weight minimization of such trusses. This scheme is tested against traditional optimization procedures for savings in computational cost. Results from both optimization methods are presented for comparison.
It is clear that computer software is needed to assist in the generation of the equations of motion for complex, flexible spacecraft. Daniel Poelaert of ESTEC has developed the software DISTEL with which he has modeled the structural dynamics for different satellites. He is interested in expanding the capabilities of DISTEL to include structural damping and control systems. Unfortunately, the software has not been released. The author has developed similar software, PDEMOD, which has been used to model the Spacecraft control Laboratory Experiment (SCOLE), the Solar Array Flight Experiment (SAFE), the Mini-MAST truss, and the LACE satellite. PDEMOD has been used also for optimal parameter estimation and integrated control-structures design. PDEMOD is also being extended to include structural damping and control systems which are imbedded into the same equations for the structural dynamics. This paper will address the formulation of the equations for the structural dynamics of spacecraft structures which are constructed of a 3-dimensional arrangement of rigid bodies and flexible beam elements. Control system dynamics are imbedded into the same equations so that model order reduction approximations are not necessary. The input data consists of the physical data of the elements and the topological information describing how the elements are connected. PDEMOD accomplishes the following: (1) automatically assembles the equations of motion for the entire structural model; (2) calculates the modal frequencies; (3) calculates the mode shapes; (4) generates perspective views of the mode shapes; and (5) forms selected transfer functions. The software PDEMOD continues to be developed to provide additional features to assist in analyzing and synthesizing control and structural systems for flexible spacecraft.
Explore the source record for details and available documents.
The Space Shuttle External Tank is covered with rigid polymeric closed-cell foam insulation to prevent ice formation, protect the metallic tank from aerodynamic heating, and control the breakup of the tank during re-entry. The cryogenic state of the tank, as well as the ascent into a vacuum environment, places this foam under significant stress. Because the loss of the foam during ascent poses a critical risk to the shuttle orbiter, there is much interest in understanding the stress state in the foam insulation and how it may contribute to fracture and debris loss. Several foam applications on the external tank have been analyzed using finite element methods. This presentation describes the approach used to model the foam material behavior and compares analytical results to experiments.
This document is a viewgraph presentation reporting on work in modeling the foam insulation of the Space Shuttle External Tank. An analytical understanding of foam mechanics is required to design against structural failure. The Space Shuttle External Tank is covered primarily with closed cell foam to: Prevent ice, Protect structure from ascent aerodynamic and engine plume heating, and Delay break-up during re-entry. It is important that the foam does not shed unacceptable debris during ascent environment. Therefore a modeling of the foam insulation was undertaken.
The behavior of arrays of coupled oscillators has been previously studied by computatinal solution of a set of non-linear differential equations describing the time dependence of each oscillator in the presence of signals coupled from neighboring oscillators.
A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.