Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Dynamic Programming”

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 253 records · Page 14

Experimental Study of Beam Dynamics in the PIP-II MEBT Prototype

The Proton Improvement Plan, Stage Two (PIP-II) [1] is a program of upgrades proposed for the Fermilab injection complex, which central part is an 800 MeV, 2 mA CW SRF linac. A prototype of the PIP-II linac front end called PIP-II Injector Test (PIP2IT) is being built at Fermilab. As of now, a 15 mA DC, 30-keV H- ion source, a 2 m-long Low Energy Beam Transport (LEBT), a 2.1 MeV CW RFQ, followed by a 10 m Medium Energy Beam Transport (MEBT) have been assembled and commissioned. The MEBT bunch-by-bunch chopping system and the requirement of a low uncontrolled beam loss put stringent limitations on the beam envelope and its variation. Measurements of transverse and longitudinal beam dynamics in the MEBT were performed in the range of 1- 10 mA of the RFQ beam current. Almost all measurements are made with 10 μs beam pulses in order to avoid damage to the beam line. This report presents measurements of the transverse optics with differential trajectories, reconstruction of the beam envelope with scrapers and an Allison emittance scanner, as well as bunch length measurements with a Fast Faraday Cup.

43 PARTICLE ACCELERATORS↗

Chemical Imaging and Dynamical Studies of Reactivity and Emergent Behavior in Complex Interfacial Systems. Final Technical Report

This research program explored the efficacy of using molecular-level manipulation, imaging and scanning tunneling spectroscopy in conjunction with supersonic molecular beam gas-surface scattering to significantly enhance our understanding of chemical processes occurring on well-characterized interfaces. One program focus was on the spatially-resolved emergent behavior of complex reaction systems as a function of the local geometry and density of adsorbate-substrate systems under reaction conditions. Another focus was on elucidating the emergent electronic and related reactivity characteristics of intentionally constructed single and multicomponent atom- and nanoparticle-based materials. We also examined emergent chirality and self-organization in adsorbed molecular systems where collective interactions between adsorbates and the supporting interface lead to spatial symmetry breaking. In many of these studies we combined the advantages of scanning tunneling (STM) and atomic force (AFM) imaging, scanning tunneling local electronic spectroscopy (STS), and reactive supersonic molecular beams to elucidate precise details of interfacial reactivity that had not been observed by more traditional surface science methods. Using these methods, it was possible to examine, for example, the differential reactivity of molecules adsorbed at different bonding sites in conjunction with how reactivity is modified by the local configuration of nearby adsorbates. At the core of this effort was the goal of significantly extending our understanding of interfacial atomic-scale interactions to create, with intent, molecular assemblies and materials with advanced chemical and physical properties. This ambitious program addressed several key topics in DOE Grand Challenge Science, including emergent chemical and physical properties in condensed phase systems, novel uses of chemical imaging, and the development of advanced reactivity concepts in combustion and catalysis including carbon management. These activities directly benefitted national science objectives in the areas of chemical energy production and advanced materials development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Generic Biogeochemical Module for Earth System Models: Next Generation BioGeoChemical Module (NGBGC), Version 1.0

Physical and biogeochemical processes regulate soil carbon dynamics and CO2 flux to and from atmosphere, influencing global climate changes. Integration of these processes into earth system models (e.g., community land models (CLM)), however, currently faces three major challenges: 1) extensive efforts are required to modify modeling structures and to rewrite computer programs to incorporate new or updated processes as new knowledge is being generated, 2) computational cost is prohibitively expensive to simulate biogeochemical processes in land models due to large variations in the rates of biogeochemical processes, and 3) various mathematical representations of biogeochemical processes exist to incorporate different aspects of fundamental mechanisms, but systematic evaluation of the different mathematical representations is difficult, if not possible. To address these challenges, we propose a new computational framework to easily incorporate physical and biogeochemical processes into land models. The new framework consists of a new biogeochemical module with a generic algorithm and reaction database so that new and updated processes can be incorporated into land models without the need to manually set up the ordinary differential equations to be solved numerically. The reaction database consists of processes of nutrient flow through the terrestrial ecosystems in plants, litter and soil. This framework facilitates effective comparison studies of biogeochemical cycles in an ecosystem using different conceptual models under the same land modeling framework. The approach was first implemented in CLM and benchmarked against simulations from the original CLM-CN code. A case study was then provided to demonstrate the advantages of using the new approach to incorporate a phosphorus cycle into the CLM model. To our knowledge, the phosphorus-incorporated CLM is a new model that can be used to simulate phosphorus limitation on the productivity of terrestrial ecosystems.

Fang, Yilin↗

YAP/TAZ Related BioMechano Signal Transduction and Cancer Metastasis

Mechanoreciprocity refers to a cell’s ability to maintain tensional homeostasis in response to various types of forces. Physical forces are continually being exerted upon cells of various tissue types, even those considered static, such as the brain. Through mechanoreceptors, cells sense and subsequently respond to these stimuli. These forces and their respective cellular responses are prevalent in regulating everything from embryogenic tissue-specific differentiation, programmed cell death, and disease progression, the last of which being the subject of extensive attention. Abnormal mechanical remodeling of cells can provide clues as to the pathological status of tissues. This becomes particularly important in cancer cells, where cellular stiffness has been recently accepted as a novel biomarker for cancer metastasis. Several studies have also elucidated the importance of cell stiffness in cancer metastasis, with data highlighting that a reversal of tumor stiffness has the capacity to revert the metastatic properties of cancer. In this review, we summarize our current understanding of extracellular matrix (ECM) homeostasis, which plays a prominent role in tissue mechanics. We also describe pathological disruption of the ECM, and the subsequent implications toward cancer and cancer metastasis. In addition, we highlight the most novel approaches toward understanding the mechanisms which generate pathogenic cell stiffness and provide potential new strategies which have the capacity to advance our understanding of one of human-kinds’ most clinically significant medical pathologies. These new strategies include video-based techniques for structural dynamics, which have shown great potential for identifying full-field, high-resolution modal properties, in this case, as a novel application.

60 APPLIED LIFE SCIENCES↗

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A feedback linearization approach to spacecraft control using momentum exchange devices

Recent developments in the area of nonlinear control theory have shown how coordiante changes in the state and input spaces can be used with nonlinear feedback to transform certain nonlinear ordinary differential equations into equivalent linear equations. These feedback linearization techniques are applied to resolve two problems arising in the control of spacecraft equipped with control moment gyroscopes (CMGs). The first application involves the computation of rate commands for the gimbals that rotate the individual gyroscopes to produce commanded torques on the spacecraft. The second application is to the long-term management of stored momentum in the system of control moment gyroscopes using environmental torques acting on the vehicle. An approach to distributing control effort among a group of redundant actuators is described that uses feedback linearization techniques to parameterize sets of controls which influence a specified subsystem in a desired way. The approach is adapted for use in spacecraft control with double-gimballed gyroscopes to produce an algorithm that avoids problematic gimbal configurations by approximating sets of gimbal rates that drive CMG rotors into desirable configurations. The momentum management problem is stated as a trajectory optimization problem with a nonlinear dynamical constraint. Feedback linearization and collocation are used to transform this problem into an unconstrainted nonlinear program. The approach to trajectory optimization is fast and robust. A number of examples are presented showing applications to the proposed NASA space station.

Dzielski, John Edward↗

Artificial Neural Network Representation of Dynamic Criticality Excursion Experiment Data

This paper outlines time dependent considerations in applying the artificial neural network (ANN) paradigm to dynamic nuclear criticality excursion experiment data. The ANN amongst a variety of uses, can apply to function fitting of multiple input-multiple output data sets, pattern recognition, cluster analysis, and time series prediction. The ANN is a component of Machine Learning (ML) data analysis. The practical application of ML to criticality safety has been funded through the Y-12 Plant Directed Research and Demonstration (PDRD) program to support criticality safety and emergency planning efforts. The initial effort considered function fitting of criticality excursion data. Recent efforts successfully extended the ANN to criticality safety code validation. The initial ANN work applied to criticality excursions centered around multivariate function fitting of static (tabular) data taken from the French Consequences Radiologiques d'un Accident de Criticite (CRAC) and SILENE reactor campaigns. The former report (LeCorche) summarized important tabular data for the first 44 CRAC experiments. The initial ANN was trained to specific time-independent inputs, then applied to predict peak fission rate and total fissions for a wide variety of fissile solutions process conditions. This work extends that effort to apply the ANN paradigm to dynamic prediction of time-dependent criticality excursion experiment data. A dynamic criticality excursion in fissile solution represents a complex multiphysics phenomenon. Factors include among other things, coupled neutronic, multiphase fluid flow, temperature, pressure, acoustic shock, boiling heat transfer, reactivity feedback, and radiolytic gas considerations. There have been many efforts to model (over time) as coupled ordinary and (space) as partial differential equations (ODE, PDE). The ANN representation noted here does not rely on solving explicit ODE or PDEs. Rather the intent is to demonstrate that as a time series, the criticality excursion can be predicted by training specific ANN architectures to experimental input and output time series exclusively. The rationale for extending this work to the time-dependent regime is important in developing a 'criticality excursion emulator' solely on the basis of learning characteristics from experiment dynamic data. An emulator is different than a computer code simulator driven by the solution of ODEs and PDEs given select 'tuning parameters' that match outputs to experiments. A criticality excursion emulator, derived exclusively from an ANN or ML, could predict most likely excursion characteristics, trained from historic dynamic experiment data, presented with input values within the range of those experiments. The work noted shows that generic time-series ANN architectures can be satisfactorily applied to the prediction of dynamic, extremely chaotic criticality excursion power given an initial time dependent input series. In reality, an excursion emulator could be driven by user-supplied initial (t=0), time independent or constant values to generate time dependent outputs. Thus ongoing work includes developing a custom ANN such that initial dynamic conditions would be predicted. Those results then serve as inputs to a customized ANN architecture currently in development. (author)

42 ENGINEERING↗

Towards performance-portability of the Albany Finite Element analysis code using the Kokkos library of Trilinos

Performance portability on heterogeneous high-performance computing (HPC) systems is a major challenge faced today by code developers: parallel code needs to execute correctly as well as with high performance on machines with different architectures, operating systems, and software libraries. The Finite Element Method (FEM) is a popular and flexible method for discretizing partial differential equations arising in a wide variety of scientific, engineering, and industry applications that require HPC. This paper presents some preliminary results pertaining to our development of a performance portable implementation of the FEM-based Albany code. Performance portability is achieved using the Kokkos library of Trilinos. We present performance results for two different physics simulations modules in Albany: the Aeras global atmosphere dynamical code and the FELIX land-ice solver. As a result, numerical experiments show that our single code implementation gives reasonable performance across two multi-core/many-core architectures: NVIDIA GPUs and multi-core CPUs.

97 MATHEMATICS AND COMPUTING↗

Differential Absorption Lidar (DIAL) Measurements from Air and Space

Differential absorption lidar (DIAL) systems have been used for the measurement of ozone, water vapor, and aerosols from aircraft platforms for over 18 years, yielding new insights into atmospheric chemistry, composition, and dynamics in large-scale field experiments conducted all over the world. The successful deployment of the lidar in-space technology experiment (LITE) in September 1994 demonstrated that space-based lidars can also collect valuable information on the global atmosphere. This paper reviews some of the contributions of the NASA Langley Research Center's airborne ozone and water vapor DIAL systems and space-based LITE system to the understanding of the atmosphere and discusses the feasibility and advantages of putting DIAL systems in space for routine atmospheric measurements of ozone and/or water vapor and aerosols and clouds. The technology and applications of the differential absorption lidar (DIAL) technique have progressed significantly since the first DIAL measurements of Schotland, and airborne DIAL measurements of ozone and water vapor are frequently being made in a wide range of field experiments. In addition, plans are underway to develop DIAL systems for use on satellites for continuous global measurements. This paper will highlight the history of airborne lidar and DIAL systems, summarize the major accomplishments of the NASA Langley DIAL program, and discuss specifications and goals for DIAL systems in space.

Browell, E. V.↗

X ray microscope/telescope test and alignment

The tasks performed by the Center for Applied Optics (CAO) in support of the Normal Incidence Multilayer X-Ray Optics Program are detailed. The Multi-Spectral Solar Telescope Array (MSSTA) was launched on a Terrier-boosted Black Brant sounding rocket from White Sands Missile Range on 13 May 1991. High resolution images of the sun in the soft x ray to extreme ultraviolet (EUV) regime were obtained with normal-incidence Cassegrain, Ritchey-Chretien, and Herschelian telescopes mounted in the sounding rocket. MSSTA represents the first use of multilayer optics to study a very broad range of x ray and EUV solar emissions. Energy-selective properties of multilayer-coated optics allow distinct groups of emission lines to be isolated in the solar corona and transition region. Features of the near and far coronal structures including magnetic loops of plasmas, coronal plumes, coronal holes, faint structures, and cool prominences are visible in these images. MSSTA successfully obtained unprecedented information regarding the structure and dynamics of the solar atmosphere in the temperature range of 10(exp 4)-10(exp 7) K. The performance of the MSSTA has demonstrated a unique combination of ultra-high spatial resolution and spectral differentiation by use of multilayer optics.

Walker, Arthur B. C.↗

UAVSAR: A New NASA Airborne SAR System for Science and Technology Research

NASA's Jet Propulsion Laboratory is currently building a reconfigurable, polarimetric L-band synthetic aperture radar (SAR), specifically designed to acquire airborne repeat track SAR data for differential interferometric measurements. Differentian interferometry can provide key deformation measurements, important for studies of earthquakes, volcanoes and other dynamically changing phenomena. Using precision real-time GPS and a sensor controlled flight management system, the system will be able to fly predefined paths with great precision. The expected performance of the flight control system will constrain the flight path to be within a 10 m diameter tube about the desired flight track. The radar will be designed to be operable on a UAV (Unpiloted Aria1 Vehicle) but will initially be demonstrated on a NASA Gulfstream III. The radar will be fully polarimetric, with a range bandwidth of 80 MHz (2 m range resolution), and will support a 16 km range swath. The antenna will be electronically steered along track to assure that the antenna beam can be directed independently, regardless of the wind direction and speed. Other features supported by the antenna include elevation monopulse and pulse-to-pulse re-steering capabilities that will enable some novel modes of operation. The system will nominally operate at 45,000 ft (13800 m). The program began as an Instrument Incubator Project (IIP) funded by NASA Earth Science and Technology Office (ESTO).

synthetic aperture radar (SAR)↗

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗

Design Optimization Study on the Single Tank Packed-bed Thermal Energy Storage System

Idaho National Laboratory (INL) is establishing the Dynamic Energy Transport and Integration Lab (DETAIL) as part of its commitment to research on nuclear-renewable hybrid energy systems and associated advanced reactor technologies. DETAIL is designed to allow several different energy systems to work in unison. Current plans include a PWR simulator, high-temperature steam electrolysis (HSTE) unit, and a renewable energy system (e.g., photovoltaics). DETAIL will provide the real world basis for studies on the systems integration and system configurations to be completed. Encompassed in the DETAIL program is the Thermal Energy Storage (TES) system. The TES system is one of the key components in DETAIL because it allows for the delayed release of energy and can be used to simulate storage capabilities currently being considered by utility providers. Of the various TES concepts, this paper deals with the single-tank packed bed TES system. Particular attention is given to the transient thermal behavior of fluid and solid particles within the packed bed thermocline tank and heat storage efficiency influenced by various design parameters. The effects of tank geometry (height-to-diameter ratio), filler size, filler packing ratio and operating temperature differentials are investigated. Based on the parametric study and cost analysis, the optimal TES tank design for DETAIL is discussed.

25 ENERGY STORAGE↗

Developments and Challenges in the Design of the ITER DRGA

We report the ITER Diagnostic Residual Gas Analyzer (DRGA) will measure the distribution of gas species, i.e., deuterium (D), tritium (T), and impurities, in the divertor exhaust stream and in the plasma periphery, with time resolution relevant to fusion plasma–wall particle dynamics. The uniqueness of the DRGA, over previous implementations of plasma dynamics residual gas analysis, is an integrated approach, combining mass and low-temperature plasma-activated optical spectroscopy, in a differentially pumped analysis station. A further unique feature of the ITER divertor-specific DRGA is an ~8-m separation of the analysis station from the sampled pumping duct, while still maintaining a ~1-s response time for hydrogen isotopic concentrations. ITER DRGA final design activities are strongly benefiting from testing of prototypical DRGA components and methods on present fusion devices, most currently on JET and W7-X. DRGA systems are implemented on both these devices and include sensors (and pumping methods) that are directly relevant to the ITER DRGA design. The recent JET-DTE2 campaign has provided the first experience on operating the combined ITER DRGA sensors with D-T plasmas. While enhancing system design for ITER, this experience on operating devices has also revealed additional engineering challenges, which further guide the continuing final design project. Meanwhile, the recent determination that the ITER DRGA, with slight optimization, will resolve the helium isotopes well enough to support an ITER pre-DT, He-3-based heating scheme, has greatly increased ITER Research Program interest in the DRGA and its implementation well ahead of the DT phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rapid Airplane Parametric Input Design (RAPID)

RAPID is a methodology and software system to define a class of airplane configurations and directly evaluate surface grids, volume grids, and grid sensitivity on and about the configurations. A distinguishing characteristic which separates RAPID from other airplane surface modellers is that the output grids and grid sensitivity are directly applicable in CFD analysis. A small set of design parameters and grid control parameters govern the process which is incorporated into interactive software for 'real time' visual analysis and into batch software for the application of optimization technology. The computed surface grids and volume grids are suitable for a wide range of Computational Fluid Dynamics (CFD) simulation. The general airplane configuration has wing, fuselage, horizontal tail, and vertical tail components. The double-delta wing and tail components are manifested by solving a fourth order partial differential equation (PDE) subject to Dirichlet and Neumann boundary conditions. The design parameters are incorporated into the boundary conditions and therefore govern the shapes of the surfaces. The PDE solution yields a smooth transition between boundaries. Surface grids suitable for CFD calculation are created by establishing an H-type topology about the configuration and incorporating grid spacing functions in the PDE equation for the lifting components and the fuselage definition equations. User specified grid parameters govern the location and degree of grid concentration. A two-block volume grid about a configuration is calculated using the Control Point Form (CPF) technique. The interactive software, which runs on Silicon Graphics IRIS workstations, allows design parameters to be continuously varied and the resulting surface grid to be observed in real time. The batch software computes both the surface and volume grids and also computes the sensitivity of the output grid with respect to the input design parameters by applying the precompiler tool ADIFOR to the grid generation program. The output of ADIFOR is a new source code containing the old code plus expressions for derivatives of specified dependent variables (grid coordinates) with respect to specified independent variables (design parameters). The RAPID methodology and software provide a means of rapidly defining numerical prototypes, grids, and grid sensitivity of a class of airplane configurations. This technology and software is highly useful for CFD research for preliminary design and optimization processes.

Smith, Robert E.↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

An EvoDevo Study of Salmonid Visual Opsin Dynamics and Photopigment Spectral Sensitivity

Salmonids are ideal models as many species follow a distinct developmental program from demersal eggs and a large yolk sac to hatching at an advanced developmental stage. Further, these economically important teleosts inhabit both marine- and freshwaters and experience diverse light environments during their life histories. At a genome level, salmonids have undergone a salmonid-specific fourth whole genome duplication event (Ss4R) compared to other teleosts that are already more genetically diverse compared to many non-teleost vertebrates. Thus, salmonids display phenotypically plastic visual systems that appear to be closely related to their anadromous migration patterns. This is most likely due to a complex interplay between their larger, more gene-rich genomes and broad spectrally enriched habitats; however, the molecular basis and functional consequences for such diversity is not fully understood. This study used advances in genome sequencing to identify the repertoire and genome organization of visual opsin genes (those primarily expressed in retinal photoreceptors) from six different salmonids [Atlantic salmon ( Salmo salar ), brown trout ( Salmo trutta ), Chinook salmon ( Oncorhynchus tshawytcha ), coho salmon ( Oncorhynchus kisutch ), rainbow trout ( Oncorhynchus mykiss ), and sockeye salmon ( Oncorhynchus nerka )] compared to the northern pike ( Esox lucius ), a closely related non-salmonid species. Results identified multiple orthologues for all five visual opsin classes, except for presence of a single short-wavelength-sensitive-2 opsin gene. Several visual opsin genes were not retained after the Ss4R duplication event, which is consistent with the concept of salmonid rediploidization. Developmentally, transcriptomic analyzes of Atlantic salmon revealed differential expression within each opsin class, with two of the long-wavelength-sensitive opsins not being expressed before first feeding. Also, early opsin expression in the retina was located centrally, expanding dorsally and ventrally as eye development progressed, with rod opsin being the dominant visual opsin post-hatching. Modeling by spectral tuning analysis and atomistic molecular simulation, predicted the greatest variation in the spectral peak of absorbance to be within the Rh2 class, with a ∼40 nm difference in λ max values between the four medium-wavelength-sensitive photopigments. Overall, it appears that opsin duplication and expression, and their respective spectral tuning profiles, evolved to maximize specialist color vision throughout an anadromous lifecycle, with some visual opsin genes being lost to tailor marine-based vision.

59 BASIC BIOLOGICAL SCIENCES↗

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.↗