SEARCH · Engineering Papers

Results for “object oriented modeling”

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

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Design of Object-Oriented Distributed Simulation Classes

Distributed simulation of aircraft engines as part of a computer aided design package being developed by NASA Lewis Research Center for the aircraft industry. The project is called NPSS, an acronym for "Numerical Propulsion Simulation System". NPSS is a flexible object-oriented simulation of aircraft engines requiring high computing speed. It is desirable to run the simulation on a distributed computer system with multiple processors executing portions of the simulation in parallel. The purpose of this research was to investigate object-oriented structures such that individual objects could be distributed. The set of classes used in the simulation must be designed to facilitate parallel computation. Since the portions of the simulation carried out in parallel are not independent of one another, there is the need for communication among the parallel executing processors which in turn implies need for their synchronization. Communication and synchronization can lead to decreased throughput as parallel processors wait for data or synchronization signals from other processors. As a result of this research, the following have been accomplished. The design and implementation of a set of simulation classes which result in a distributed simulation control program have been completed. The design is based upon MIT "Actor" model of a concurrent object and uses "connectors" to structure dynamic connections between simulation components. Connectors may be dynamically created according to the distribution of objects among machines at execution time without any programming changes. Measurements of the basic performance have been carried out with the result that communication overhead of the distributed design is swamped by the computation time of modules unless modules have very short execution times per iteration or time step. An analytical performance model based upon queuing network theory has been designed and implemented. Its application to realistic configurations has not been carried out.

Schoeffler, James D.↗

NASA NTRS2002

From following edges to pursuing objects

Primates can generate accurate, smooth eye-movement responses to moving target objects of arbitrary shape and size, even in the presence of complex backgrounds and/or the extraneous motion of non-target objects. Most previous studies of pursuit have simply used a spot moving over a featureless background as the target and have thus neglected critical issues associated with the general problem of recovering object motion. Visual psychophysicists and theoreticians have shown that, for arbitrary objects with multiple features at multiple orientations, object-motion estimation for perception is a complex, multi-staged, time-consuming process. To examine the temporal evolution of the motion signal driving pursuit, we recorded the tracking eye movements of human observers to moving line-figure diamonds. We found that pursuit is initially biased in the direction of the vector average of the motions of the diamond's line segments and gradually converges to the true object-motion direction with a time constant of approximately 90 ms. Furthermore, transient blanking of the target during steady-state pursuit induces a decrease in tracking speed, which, unlike pursuit initiation, is subsequently corrected without an initial direction bias. These results are inconsistent with current models in which pursuit is driven by retinal-slip error correction. They demonstrate that pursuit models must be revised to include a more complete visual afferent pathway, which computes, and to some extent latches on to, an accurate estimate of object direction over the first hundred milliseconds or so of motion.

NASA Discipline Neuroscience↗

NASA NTRS2005

Measuring Low-Order Aberrations in a Segmented Telescope

The in-focus PSF optimizer (IPO) is an algorithm for use in monitoring and controlling the alignment of the segments of a segmented-mirror astronomical telescope. IPO is so named because it computes wave-front aberrations of the telescope from digitized pointspread functions (PSFs) measured in infocus images. Inasmuch as distant astronomical objects that behave optically as point sources can typically be seen in almost any astronomical image, the main benefit afforded by IPO may be to enable maintenance of mirror-segment alignments without detracting from valuable scientific-observation time. IPO evolved from prescription-retrieval type algorithms. Prescription retrieval uses in-focus and out-of-focus PSFs to infer the state of an imaging optical system. The state, in this context, refers to the positions, orientations, and low-order figure errors of the optical elements in the system. Both prescription- retrieval and IPO use an iterative, nonlinear, least-squares optimizer to compute the optimal state parameters such that a digital computer-generated model image matches the digitized image acquired from the real system. The difference between IPO and prescription- retrieval algorithms is that IPO is specifically designed to utilize infocus images only. Although the restriction to in-focus images limits IPO to calculating only the lowest-order wave front aberrations, it also causes the resulting computation to take much less time because fewer degrees of freedom are included in the optimization process. In the prescription retrieval software developed at JPL, the model images are generated using the ray-trace/physical optics program, MACOS. IPO, on the other hand, uses a linear sensitivity matrix to compute the exit-pupil wave front from the system parameters; the wave front is then converted into a complex pupil field, which is then propagated to the image plane via a fast Fourier transform. This approach is computationally faster and requires less computer memory than is needed for prescription retrieval.

Ohara, Catherine↗

NASA NTRS2005

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗

NASA NTRS2007

Using EIGER for Antenna Design and Analysis

EIGER (Electromagnetic Interactions GenERalized) is a frequency-domain electromagnetics software package that is built upon a flexible framework, designed using object-oriented techniques. The analysis methods used include moment method solutions of integral equations, finite element solutions of partial differential equations, and combinations thereof. The framework design permits new analysis techniques (boundary conditions, Green#s functions, etc.) to be added to the software suite with a sensible effort. The code has been designed to execute (in serial or parallel) on a wide variety of platforms from Intel-based PCs and Unix-based workstations. Recently, new potential integration scheme s that avoid singularity extraction techniques have been added for integral equation analysis. These new integration schemes are required for facilitating the use of higher-order elements and basis functions. Higher-order elements are better able to model geometrical curvature using fewer elements than when using linear elements. Higher-order basis functions are beneficial for simulating structures with rapidly varying fields or currents. Results presented here will demonstrate curren t and future capabilities of EIGER with respect to analysis of installed antenna system performance in support of NASA#s mission of exploration. Examples include antenna coupling within an enclosed environment and antenna analysis on electrically large manned space vehicles.

Champagne, Nathan J.↗

NASA NTRS2015

Spherical Coordinate Systems for Streamlining Suited Mobility Analysis

Introduction: When describing human motion, biomechanists generally report joint angles in terms of Euler angle rotation sequences. However, there are known limitations in using this method to describe complex motions such as the shoulder joint during a baseball pitch. Euler angle notation uses a series of three rotations about an axis where each rotation is dependent upon the preceding rotation. As such, the Euler angles need to be regarded as a set to get accurate angle information. Unfortunately, it is often difficult to visualize and understand these complex motion representations. It has been shown that using a spherical coordinate system allows Anthropometry and Biomechanics Facility (ABF) personnel to increase their ability to transmit important human mobility data to engineers, in a format that is readily understandable and directly translatable to their design efforts. Objectives: The goal of this project was to use innovative analysis and visualization techniques to aid in the examination and comprehension of complex motions. Methods: This project consisted of a series of small sub‐projects, meant to validate and verify a new method before it was implemented in the ABF's data analysis practices. A mechanical test rig was built and tracked in 3D using an optical motion capture system. Its position and orientation were reported in both Euler and spherical reference systems. In the second phase of the project, the ABF estimated the error inherent in a spherical coordinate system, and evaluated how this error would vary within the reference frame. This stage also involved expanding a kinematic model of the shoulder to include the rest of the joints of the body. The third stage of the project involved creating visualization methods to assist in interpreting motion in a spherical frame. These visualization methods will be incorporated in a tool to evaluate a database of suited mobility data, which is currently in development. Results: Initial results demonstrated that a spherical coordinate system is helpful in describing and visualizing the motion of a space suit. The system is particularly useful in describing the motion of the shoulder, where multiple degrees of freedom can lead to very complex motion paths.

Benson, Elizabeth↗

NASA NTRS1990

Robot tracking system improvements and visual calibration of orbiter position for radiator inspection

The following separate topics are addressed: (1) improving a robotic tracking system; and (2) providing insights into orbiter position calibration for radiator inspection. The objective of the tracking system project was to provide the capability to track moving targets more accurately by adjusting parameters in the control system and implementing a predictive algorithm. A computer model was developed to emulate the tracking system. Using this model as a test bed, a self-tuning algorithm was developed to tune the system gains. The model yielded important findings concerning factors that affect the gains. The self-tuning algorithms will provide the concepts to write a program to automatically tune the gains in the real system. The section concerning orbiter position calibration provides a comparison to previous work that had been performed for plant growth. It provided the conceptualized routines required to visually determine the orbiter position and orientation. Furthermore, it identified the types of information which are required to flow between the robot controller and the vision system.

Tonkay, Gregory↗

NASA NTRS2006

Calculation of Zonal Winds using Accelerometer and Rate Data from Mars Global Surveyor

The Mars Global Surveyor spacecraft was initially placed into a high eccentricity, nearly polar orbit about Mars with a 45-hour period. To accomplish the science objectives of the mission, a 2-hour, circular orbit was required. Using a method known as aerobraking, numerous passes through the upper atmosphere slowed the spacecraft, thereby reducing the orbital period and eccentricity. To successfully perform aerobraking, the spacecraft was designed to be longitudinally, aerodynamically stable in pitch and yaw. Since the orbit is nearly polar, the yaw orientation of the spacecraft was sensitive to disturbances caused by the zonal components of wind (east-to-west or west-to-east) acting on the spacecraft at aerobraking altitudes. Zonal wind velocities were computed by equating the aerodynamic and inertia-related torques acting on the spacecraft. Comparisons of calculated zonal winds with those computed from the Mars Thermospheric Global Circulation Model are discussed.

Mars thermospheric winds↗

NASA NTRS1999

A Comprehensive Orbit Reconstruction for the Galileo Prime Mission in the J2000 System

The Galileo spacecraft arrived at Jupiter in December of 1995 to begin an orbital tour of the Jovian system. The objective of the tour was up close study of the planet, its satellites, and its magnetosphere. The spacecraft completed its 11 orbit prime mission in November of 1997 having had 16 successful close encounters with the Galilean satellites (including two prior to Jupiter orbit insertion). Galileo continues to operate and will have made an additional 10 orbits of Jupiter by the date of this Conference. Earlier papers discuss the determination of the spacecraft orbit in support of mission operations from arrival at Jupiter through the first 9 orbits. In this paper we re-examine those earlier orbits and extend the analysis through orbit 12, the first orbit of the Galileo Europa Mission (GEM). The objective of our work is the reconstruction of the spacecraft trajectory together with the development of a consistent set of ephemerides for the Galilean satellites. As a necessary byproduct of the reconstruction we determine improved values for the Jovian system gravitational parameters and for the Jupiter pole orientation angles. Our preliminary analyses have already led to many of the results reported in the scientific literature. Unlike the Galileo Navigation Team which operates in the EME-1950 coordinate system, we elected to work in the (J2000) International Celestial Reference Frame (ICRF), the reference frame of the current JPL planetary and satellite ephemerides as well as the standard frame of the international astronomical and planetary science community. Use of this frame permits more precise modelling of the spacecraft and satellite observations. Moreover, it is the frame of choice for all other operational JPL missions and will probably be the frame for future missions for some time. Consequently, our adoption of the ICRF will facilitate the combination of our results with any obtained from future missions (e.g. the proposed Europa Orbiter mission). In addition, our results may be used by the science community, without need of a reference frame conversion.

Jacobson, Robert A.↗

NASA NTRS1991

NASA 1990 Multisensor Airborne Campaigns (MACs) for ecosystem and watershed studies

The Multisensor Airborne Campaign (MAC) focus within NASA's former Land Processes research program was conceived to achieve the following objectives: to acquire relatively complete, multisensor data sets for well-studied field sites, to add a strong remote sensing science component to ecology-, hydrology-, and geology-oriented field projects, to create a research environment that promotes strong interactions among scientists within the program, and to more efficiently utilize and compete for the NASA fleet of remote sensing aircraft. Four new MAC's were conducted in 1990: the Oregon Transect Ecosystem Research (OTTER) project along an east-west transect through central Oregon, the Forest Ecosystem Dynamics (FED) project at the Northern Experimental Forest in Howland, Maine, the MACHYDRO project in the Mahantango Creek watershed in central Pennsylvania, and the Walnut Gulch project near Tombstone, Arizona. The OTTER project is testing a model that estimates the major fluxes of carbon, nitrogen, and water through temperate coniferous forest ecosystems. The focus in the project is on short time-scale (days-year) variations in ecosystem function. The FED project is concerned with modeling vegetation changes of forest ecosystems using remotely sensed observations to extract biophysical properties of forest canopies. The focus in this project is on long time-scale (decades to millenia) changes in ecosystem structure. The MACHYDRO project is studying the role of soil moisture and its regulating effects on hydrologic processes. The focus of the study is to delineate soil moisture differences within a basin and their changes with respect to evapotranspiration, rainfall, and streamflow. The Walnut Gulch project is focused on the effects of soil moisture in the energy and water balance of arid and semiarid ecosystems and their feedbacks to the atmosphere via thermal forcing.

Wickland, Diane E.↗

DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Deflection of Solar Energetic Protons by an Ion Thruster

While Galactic Cosmic Rays (GCR) represent a long-term cancer risk to deep space human missions, Solar Energetic Particles (SEP) represent a significant acute radiation hazard to human spaceflight. This work investigates wave-particle interactions caused by the incidence of a beam of energetic protons, representing a notional SEP event, onto a cloud of ions representing a notional ion thruster plume. The resulting electromagnetic ion/ion instability is developed to investigate the degree to which energetic protons can be scattered out of the initial beam. This work includes investigations of deflection rate dependences on incident SEP energy and flux, ion thruster power and orientation relative to SEP incidence, plume ion species, and background plasma. Also investigated are the spectrum and angular distribution of electromagnetic radiation resulting from the wave-particle interactions. The primary objective of this investigation is to develop a technique for active mitigation of SEP radiation threats to spacecraft and crew, by exploring the wave-particle interaction of a SEP event with the plasma plume expelled from an ion thruster. Since the cost in expended fuel restricts this technique to intermittent use during interplanetary flights, effective prediction of when a SEP event might impact the mission is still important, and the observation of conditions that might lead to SEP events, human-in-the-loop space weather analyses, and development of SEP prediction models remain essential to crew safety, spacecraft function, and mission success.

Interplanetary Space Flight↗

NASA NTRS1990

Photovoltaic array space power plus diagnostics experiment

The objective of the Photovoltaic Array Space Power Plus Diagnostics (PASP Plus) experiment is to measure the effects of the interaction of the low- to mid-altitude space environment on the performance of a diverse set of small solar-cell arrays (planar and concentrator, representative of present and future military technologies) under differing conditions of velocity-vector orientation and simulated (by biasing) high-voltage operation. Solar arrays to be tested include Si and GaAs planar arrays and several types of GaAs concentrator arrays. Diagnostics (a Langmuir probe and a pressure gauge) and a transient pulse monitor (to measure radiated and conducted EMI during arcing) will be used to determine the impact of the environment on array operation to help verify various interactions models. Results from a successful PASP Plus flight will furnish answers to important interactions questions and provide inputs for design and test standards for photovoltaic space-power subsystems.

Guidice, Donald A.↗

NASA NTRS1994

Measurement and characterization of force dynamics in high T(sub c) superconductors

Magnetic bearing implementations using more exotic superconducting phenomena have been proliferating in recent years because they have important advantages over conventional implementations. For example, the stable suspension of a six degrees-of-freedom object by superconducting means can be achieved without a control system and with the use of only a single superconductor. It follows that the construction becomes much simpler with decreased need for position sensors and stabilizers. However, it is recognized that the design of superconducting systems can be difficult because important characteristics relating to the 6 degree-of-freedom dynamics of an object suspended magnetically are not readily available and the underlying principles of superconducting phenomena are not yet completely understood. To eliminate some of the guesswork in the design process, this paper proposes a system which can resolve the mechanical properties of suspension by superconductivity and provide position and orientation dependent data about the system's damping, stiffness, and frequency response characteristics. This system employs an actively-controlled magnetically-suspended fine-motion device that can also be used as a six degree-of-freedom force sensor. By attaching the force sensor to a permanent magnet that is being levitated above a superconducting magnet, mechanical characteristics of the superconductor levitation can be extracted. Such information would prove useful for checking the validity of theoretical models and may even give insights into superconducting phenomena.

Higuchi, Toshiro↗

DOE OSTI2021

Modelling Nuclear Thermal Propulsion Reactor Startup Transients

The National Aeronautics and Space Administration (NASA) has set the goal of a manned mission to Mars by the year 2030 [1] and charged the national academy of sci- ences "to identify primary technical and programmatic chal- lenges, merits, and risks for maturing space nuclear propulsion technologies of interest to a future human Mars exploration mission" [2]. One relevant technology, nuclear thermal propul- sion (NTP), has notable advantages over traditional chemical rockets; most important among them is the ability to produce larger specific impulse on the order of 900s. The reduction of mission time is crucial for a manned mission to Mars to reduce the risk for the crew. Due to its higher specific impulse, NTP satisfies this need and is pursued as one technology to get humans to Mars [3, 4]. The construction of an NTP sys- tem has to negotiate several challenges laid out in Ref. [2]; one of these challenges is the need to startup the NTP sys- tem from essentially cold conditions to full power within one minute. This paper focuses on studying the neutronics and thermal-hydraulics behavior of a simplified NTP model dur- ing prescribed reactivity insertions and mass flow rate (MFR) ramps. It is the goal of this paper to investigate startup, peak material temperatures, and average specific impulse for a low enriched Uranium (LEU), ceramic and metal material (CER- MET) NTP system when varying reactivity insertion and MFR ramps.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1993

System identification of the Large-Angle Magnetic Suspension Test Facility (LAMSTF)

The Large-Angle Magnetic Suspension Test Facility (LAMSTF), a laboratory-scale research project to demonstrate the magnetic suspension of objects over wide ranges of attitudes, has been developed. This system represents a scaled model of a planned Large-Gap Magnetic Suspension System (LGMSS). The LAMSTF system consists of a planar array of five copper electromagnets which actively suspend a small cylindrical permanent magnet. The cylinder is a rigid body and can be controlled to move in five independent degrees of freedom. Five position variables are sensed indirectly by using infra-red light-emitting diodes and light-receiving phototransistors. The motion of the suspended cylinder is in general nonlinear and hence only the linear, time-invariant perturbed motion about an equilibrium state is considered. One of the main challenges in this project is the control of the suspended element over a wide range of orientations. An accurate dynamic model plans an essential role in controller design. The analytical model of the LAMSTF system includes highly unstable real poles (about 10 Hz) and low-frequency flexible modes (about 0.16 Hz). Projection filters are proposed to identify the state space model from closed-loop test data in time domain. A canonical transformation matrix is also derived to transform the identified state space model into the physical coordinate. The LAMSTF system is stabilized by using a linear quadratic regulator (LQR) feedback controller. The rate information is obtained by calculating the back difference of the sensed position signals. The reference inputs contain five uncorrelated random signals. This control input and the system reponse are recorded as input/output data to identify the system directly from the projection filters. The sampling time is 4 ms and the model is fairly accurate in predicting the step responses for different controllers while the analytical model has a deficiency in the pitch axis.

Huang, Jen-Kuang↗

DOE OSTI2025

High-Performance Computing Based EMT Simulation: Power Grid with IBRs

Electromagnetic transient (EMT) simulation of power grids with high-fidelity models of inverter-based resources (IBRs) is time-consuming and difficult to scale. The necessity for high-fidelity models of IBRs that incorporate the dynamics of individual inverters within IBRs has been showcased in recent studies. These studies focused on events with partial power reduction in each IBR during a transmission line fault in the power grid. These types of events have been documented in multiple North American Electric Reliability Council (NERC) reports in the past decade. It is imperative then to find solutions to speed-up EMT simulations and scale the size of the region with IBRs studied in EMT simulations. In this paper, a combination of numerical simulation algorithms with high-performance computing techniques are employed in discretization and linear solvers employed in the proposed RE-INTEGRATE EMT simulation platform for power grid with IBRs. For ease of scalability, modular and object-oriented programming is used as these techniques are implemented. Additionally, automation software is developed to convert legacy software codes to the proposed RE-INTEGRATE EMT simulation platform. Thereafter, this platform is evaluated on multi-core central processing units (CPUs). Finally, scale-up tests are performed to showcase the scalability that is possible.

Marthi, Phani Ratna Vanamali [ORNL] (ORCID:0000000↗

DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗