Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiscale dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Mechanical Properties of Graphene Nanoplatelet/Carbon Fiber/Epoxy Hybrid Composites: Multiscale Modeling and Experiments

Because of the relatively high specific mechanical properties of carbon fiber/epoxy composite materials, they are often used as structural components in aerospace applications. Graphene nanoplatelets (GNPs) can be added to the epoxy matrix to improve the overall mechanical properties of the composite. The resulting GNP/carbon fiber/epoxy hybrid composites have been studied using multiscale modeling to determine the influence of GNP volume fraction, epoxy crosslink density, and GNP dispersion on the mechanical performance. The hierarchical multiscale modeling approach developed herein includes Molecular Dynamics (MD) and micromechanical modeling, and it is validated with experimental testing of the same hybrid composite material system. The results indicate that the multiscale modeling approach is accurate and provides physical insight into the composite mechanical behavior. Also, the results quantify the substantial impact of GNP volume fraction and dispersion on the transverse mechanical properties of the hybrid composite while the effect on the axial properties is shown to be insignificant.

nanocomposites↗

Mechanical Properties of Graphene Nanoplatelet Carbon Fiber Epoxy Hybrid Composites: Multiscale Modeling and Experiments

Because of the relatively high specific mechanical properties of carbon fiber/epoxy composite materials, they are often used as structural components in aerospace applications. Graphene nanoplatelets (GNPs) can be added to the epoxy matrix to improve the overall mechanical properties of the composite. The resulting GNP/carbon fiber/epoxy hybrid composites have been studied using multiscale modeling to determine the influence of GNP volume fraction, epoxy crosslink density, and GNP dispersion on the mechanical performance. The hierarchical multiscale modeling approach developed herein includes Molecular Dynamics (MD) and micromechanical modeling, and it is validated with experimental testing of the same hybrid composite material system. The results indicate that the multiscale modeling approach is accurate and provides physical insight into the composite mechanical behavior. Also, the results quantify the substantial impact of GNP volume fraction and dispersion on the transverse mechanical properties of the hybrid composite, while the effect on the axial properties is shown to be insignificant.

computational mechanics↗

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Epidermal Homeostasis and Radiation Responses in a Multiscale Tissue Modeling Framework

The surface of skin is lined with several thin layers of epithelial cells that are maintained throughout life time by a small population of stem cells. High dose radiation exposures could injure and deplete the underlying proliferative cells and induce cutaneous radiation syndrome. In this work we propose a multiscale computational model for skin epidermal dynamics that links phenomena occurring at the subcellular, cellular, and tissue levels of organization, to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis. Incorporating experimentally measured histological and cell kinetic parameters, we obtain results of population kinetics and proliferation indexes comparable to observations in unirradiated and acutely irradiated swine experiments. At the sub-cellular level, several recently published Wnt signaling controlled cell-cycle models are applied and the roles of key components and parameters are analyzed. Based on our simulation results, we demonstrate that a moderate increase of proliferation rate for the survival proliferative cells is sufficient to fully repopulate the area denuded by high dose radiation, as long as the integrity of underlying basement membrane is maintained. Our work highlights the importance of considering proliferation kinetics as well as the spatial organization of tissues when conducting in vivo investigations of radiation responses. This integrated model allow us to test the validity of several basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and enhance our understanding of the pathophysiological effects of ionizing radiation on skin.

Hu, Shaowen↗

Newton-Euler Dynamic Equations of Motion for a Multi-body Spacecraft

The Magnetospheric MultiScale (MMS) mission employs a formation of spinning spacecraft with several flexible appendages and thruster-based control. To understand the complex dynamic interaction of thruster actuation, appendage motion, and spin dynamics, each spacecraft is modeled as a tree of rigid bodies connected by spherical or gimballed joints. The method presented facilitates assembling by inspection the exact, nonlinear dynamic equations of motion for a multibody spacecraft suitable for solution by numerical integration. The building block equations are derived by applying Newton's and Euler's equations of motion to an "element" consisting of two bodies and one joint (spherical and gimballed joints are considered separately). Patterns in the "mass" and L'force" matrices guide assembly by inspection of a general N-body tree-topology system. Straightforward linear algebra operations are employed to eliminate extraneous constraint equations, resulting in a minimum-dimension system of equations to solve. This method thus combines a straightforward, easily-extendable, easily-mechanized formulation with an efficient computer implementation.

Stoneking, Eric↗

SynopFrame: Multiscale time-dependent visual abstraction framework for analyzing DNA nanotechnology simulations

We present an open-source framework, SynopFrame, that allows DNA nanotechnology (DNA-nano) experts to analyze and understand molecular dynamics simulation trajectories of their designs. We use a multiscale multi-dimensional abstraction space, connect the representations to a projected conformational space plot of the structure’s temporal sequence, and thus enable experts to analyze the dynamics of their structural designs and, specifically, failure cases of the assembly. In addition, our time-dependent abstraction representation allows the biologists, for the first time in a smooth and structurally clear way, to identify and observe temporal transitions of a DNA-nano design from one configuration to another, and to highlight important periods of the simulation for further analysis. We realize SynopFrame as a dashboard of the different synchronized 3D spatial and 2D schematic visual representations, with a color overlay to show essential properties such as the status of hydrogen bonds. The linking of the spatial, schematic, and abstract views ensures that users can effectively analyze the high-frequency motion. We also categorize the status of the hydrogen bonds into a new format to allow us to color-encode it and overlay it on the representations. To demonstrate the utility of SynopFrame, we describe example usage scenarios and report user feedback.

Abstraction space↗

Multiscale Modeling of Thermal Conductivity of Polymer/Carbon Nanocomposites

Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between nanoparticles and amorphous and crystalline polymer matrices. Bulk thermal conductivities of the nanocomposites were then estimated using an established effective medium approach. To study functionalization, oligomeric ethylene-vinyl alcohol copolymers were chemically bonded to a single wall carbon nanotube. The results, in a poly(ethylene-vinyl acetate) matrix, are similar to those obtained previously for grafted linear hydrocarbon chains. To study the effect of noncovalent functionalization, two types of polyethylene matrices. -- aligned (extended-chain crystalline) vs. amorphous (random coils) were modeled. Both matrices produced the same interfacial thermal resistance values. Finally, functionalization of edges and faces of plate-like graphite nanoparticles was found to be only modestly effective in reducing the interfacial thermal resistance and improving the composite thermal conductivity

Clancy, Thomas C.↗

Towards a multiscale approach for understanding irradiation induced swelling and creep in 316 stainless steels - A coupled cluster dynamics and crystal plasticity approach

Structural materials undergo mechanical degradation, in part due to irradiation-induced swelling and creep, under nuclear reactor conditions. While swelling results from the migration and clustering of irradiation-induced atomic-scale mobile defects, the interaction of mesoscale dislocations with these defects causes creep deformation. A coupled crystal plasticity (CP) and mean-field cluster dynamics (CD) approach is presented to investigate the effect of irradiation on the long-term mechanical behavior of 316 stainless steel, which are under consideration for use in nuclear reactors. The temporal evolution of Frenkel pairs and extended defect population, under a chosen irradiation flux and temperature, is predicted using the CD model. The impact of the irradiation defects on the stress state, and the resulting dislocation-mediated inelastic deformation, is modeled concurrently with the CP model. The inelastic deformation is irradiation flux dependent, and early-stage defect evolution determines the later-stage mechanical behavior in 316 stainless steel.

36 - MATERIALS SCIENCE↗

MultiSector Dynamics: 2023 Inaugural Workshop Report

Preface The MultiSector Dynamics (MSD) Community of Practice (CoP) hosted an inaugural workshop on October 3-5, 2023 at the University of California, Davis, to bring together members of the MSD community of practice to advance understanding of the co-evolution of human and natural systems, and to build the next generation of tools that bridge sectors, scales, and systems to realize a more resilient and equitable future. The theme of the workshop was "Advancing Complex Adaptive Human-Earth Systems Science in a World of Interconnected Risks". This document outlines the motivation for the workshop, its goals and objectives, the application process, the agenda, overviews of the training sessions offered to the workshop participants and a summary of each breakout session. The MSD workshop report further discusses the feedback from workshop participants and presents some reflections and next steps. The MSD Workshop organizers thank the DOE Office of Science, Earth and Environmental System Modeling, MultiSector Dynamics program area for financial support of its activities through the Integrated Multisector Multiscale Modeling (IM3) project. For more information related to the broader DOE MultiSector Dynamics Program please see https://climatemodeling.science.energy.gov/program-area/multisector-dynamics. D.L.M. and C.M.B. acknowledge support from the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory (ORNL), managed by UT-Battelle, LLC, for the US Department of Energy (DOE). Disclaimer This report was prepared as an account of work sponsored by an agency of the United States Government. Neither theUnited States Government nor any agency thereof, nor Battelle Memorial Institute, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, complete- ness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial products, process, or service by trade name,trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Pacific Northwest National Laboratory operated by Battelle for the United States Department of Energy Available from:Office of Scientific and Technical Information http://www.OSTI.gov multisectordynamics.org This work is made available under the terms of the Creative Commons Attribution- NonCommercial 4.0 International (CC BY-NC 4.0) https://creativecommons.org/licenses/by-nc/4 Suggested citation: Monier, E., Reed, P.M., Vernon, C.R., Hadjimichael, A., Brelsford, C.M., Burleyson, C.B., Dyreson, A.R., Fletcher, S.M., Giang, A., Gupta, R.S., Jackson, N.D., Jones, A.D., Lamontagne, J.R., McCollum, D.L., Morris, J.F., Moss, R.H., Peng, W., Saari, R.K., Srikrishnan, V., Szinai, J.K., Yoon, J. (2024) MultiSector Dynamics: 2023 Inaugural Workshop Report. MSD-LIVE Data Repository. doi:10.57931/2371710.

Monier, Erwan↗

Potential quantum advantage for simulation of fluid dynamics

Numerical simulation of turbulent fluid dynamics needs to either parametrize turbulence—which introduces large uncertainties—or explicitly resolve the smallest scales—which is prohibitively expensive. Here, we provide evidence through analytic bounds and numerical studies that a potential quantum speedup can be achieved to simulate fluid dynamics using quantum computing. Specifically, we provide a lattice Boltzmann formulation of fluid dynamics for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems. This is achieved via a combination of reformulating the Navier-Stokes nonlinearity (u·$\triangledown$u) to lattice-Boltzmann nonlinearity (u 2 ) and accurately linearizing the dynamical equations, which effectively trades nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this, we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size compared with polynomial scaling in the best known classical algorithms. In this paper, we suggest that a quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.

42 ENGINEERING↗

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗

Center of Mass Estimation for a Spinning Spacecraft Using Doppler Shift of the GPS Carrier Frequency

A sequential filter is presented for estimating the center of mass (CM) of a spinning spacecraft using Doppler shift data from a set of onboard Global Positioning System (GPS) receivers. The advantage of the proposed method is that it is passive and can be run continuously in the background without using commanded thruster firings to excite spacecraft dynamical motion for observability. The NASA Magnetospheric Multiscale (MMS) mission is used as a test case for the CM estimator. The four MMS spacecraft carry star cameras for accurate attitude and spin rate estimation. The angle between the spacecraft nominal spin axis (for MMS this is the geometric body Z-axis) and the major principal axis of inertia is called the coning angle. The transverse components of the estimated rate provide a direct measure of the coning angle. The coning angle has been seen to shift slightly after every orbit and attitude maneuver. This change is attributed to a small asymmetry in the fuel distribution that changes with each burn. This paper shows a correlation between the apparent mass asymmetry deduced from the variations in the coning angle and the CM estimates made using the GPS Doppler data. The consistency between the changes in the coning angle and the CM provides validation of the proposed GPS Doppler method for estimation of the CM on spinning spacecraft.

Spacecraft↗

Center of Mass Estimation for a Spinning Spacecraft Using Doppler Shift of the GPS Carrier Frequency

A sequential filter is presented for estimating the center of mass (CM) of a spinning spacecraft using Doppler shift data from a set of onboard Global Positioning System (GPS) receivers. The advantage of the proposed method is that it is passive and can be run continuously in the background without using commanded thruster firings to excite spacecraft dynamical motion for observability. The NASA Magnetospheric Multiscale (MMS) mission is used as a test case for the CM estimator. The four MMS spacecraft carry star cameras for accurate attitude and spin rate estimation. The angle between the spacecraft nominal spin axis (for MMS this is the geometric body Z-axis) and the major principal axis of inertia is called the coning angle. The transverse components of the estimated rate provide a direct measure of the coning angle. The coning angle has been seen to shift slightly after every orbit and attitude maneuver. This change is attributed to a small asymmetry in the fuel distribution that changes with each burn. This paper shows a correlation between the apparent mass asymmetry deduced from the variations in the coning angle and the CM estimates made using the GPS Doppler data. The consistency between the changes in the coning angle and the CM provides validation of the proposed GPS Doppler method for estimation of the CM on spinning spacecraft.

Estimation↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale repeating unit cell model of a single spherulite containing four disparate length scales was developed to predict the thermoelastic behavior of semicrystalline thermoplastic materials for composite aerospace applications. The continuum level scales were fully coupled and modeled using the generalized method of cells and the high fidelity generalized method of cells micromechanics theories. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents in the multiscale continuum models. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, showing good agreement with the available experimental data from the open literature. Moreover, it is shown that predicted properties are fairly insensitive to the fidelity of the micromechanics model used at the highest continuum scale or the assumed shape of the spherulite.

thermoplastics↗