Search NASASearch

SEARCH · Search NASA

Results for “Molecular 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.

558 records · Page 25

Probing the Environment around GW170817 with DESI: Insights on Galaxy Group Peculiar Velocities for Standard Siren Measurements

We present a new measurement of the Hubble constant, H 0 , following the gravitational-wave event GW170817 and Dark Energy Spectroscopic Instrument (DESI) observations. A standard siren measurement with a nearby (luminosity distance ∼40 Mpc) event such as GW170817 is typically sensitive to the peculiar motion of the host galaxy owing to local dynamics. Previous measurements from this event have taken advantage of peculiar velocity measurements of nearby galaxies, including a handful of objects in the galaxy group that the host of the event, NGC 4993, has been associated with. Still, the group’s properties and NGC 4993’s membership were debated. We present DESI observations of thousands of galaxies in the vicinity of NGC 4993, resulting in 39 group galaxies and a fivefold increase in galaxies compared to previous observations, with many contributing to a peculiar velocity measurement. Examining the local dynamics, our observations support the presence of a galaxy group of which NGC 4993 is a part with a halo mass of order ∼10 13 M ⊙ . Using peculiar velocity measurements from our fundamental plane galaxy observations, we find $H_0 = 70.9^{+6.4}_{-8.5}$ km s −1 Mpc −1 . In addition, using a peculiar velocity measurement for NGC 4993 from surface brightness fluctuations in Cosmicflows-4, we find $H_0 = 73.4^{+3.3}_{-3.9}$ km s −1 Mpc −1 . We study the impact of different galaxy selection criteria on the determination of the peculiar velocity and, in turn, on the H 0 measurement. Our results demonstrate the value of multiplexed spectroscopic observations for probing the local environments of gravitational-wave events used in standard siren measurements.

Amsellem, A. J. [Carnegie Mellon University, Pitts

Development of global/chemistry model for jet-fuel thermal stability based on observations from static and flowing experiments

Two global-chemistry models for oxidative deposition of jet fuels are evaluated by integrating them into a Computational Fluid Dynamics with Chemistry (CFDC) code. A previously developed two-step global-chemistry model was found to be insufficient to describe the thermal-oxidation and -deposition rates associated with a Jet-A fuel. A new global-chemistry model has been developed systematically based on observations from flowing and static experiments. The global-autoxidation reaction is modified such that the reaction rate becomes zeroth-order with respect to the dissolved oxygen concentration. The generation of deposit-forming precursor is coupled with the autoxidation reaction by introducing a radical species ROO. A formulation for the sticking probability has also been developed. Deposition profiles are well represented by this new model under a variety of temperature and flow conditions. The model correctly predicts the changes in magnitude and spatial location of the deposition peak due to changes in flow. The CFDC model, which is designed for flowing systems, has been extended to static experiments. The model incorporates a non-depleting species F(sub s) representing all non-oxygen compounds responsible for deposition. Static experiments were found to provide a useful and inexpensive method for estimating the concentration of F(sub s) in the fuel.

V. R. Katta

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science

Interface Evolutionand Long-Term Performance of NegativeCarbon Fiber Structural Electrodes

Abstract Laminated structural batteries present a transformative solution to reducing weight constraints in electric vehicles. These structural batteries are based on a multifunctional material that incorporates an energy storage function within a carbon fiber-reinforced polymer. Despite the potential of this technology, the intricate morphology of fiber–matrix or electrode–electrolyte interfaces and the impact of long-term cycling at low current rates (C-rates) on these interfaces remain insufficiently understood. This study addresses these critical knowledge gaps by examining the influence of matrix composition on the long-term electrochemical performance of structural battery electrodes and exploring advanced techniques to investigate carbon fiber–matrix interfaces. Localized imaging and X-ray scattering techniques were used to characterize morphological changes at the electrode–electrolyte interfaces by analyzing negative structural electrodes. The findings revealed that the matrix composition influences long-term electrochemical behavior and fiber–matrix interface formation. While the intrinsic properties of carbon fibers largely remain unaffected by long-term cycling, cycling promotes debonding at fiber–matrix interfaces. Nonetheless, residual regions of adhesion persist, underscoring the potential for preserving multifunctionality even under prolonged cycling conditions. These insights advance the understanding of interface dynamics, which is critical for optimizing structural battery technologies.

Chemistry

Changes in Arctic and Antarctic Sea-Ice Properties and Processes

Although sea ice covers ~9% of the world’s oceans, its presence, and absence, vastly alter the energetic balance of Earth’s climate. Because of geographic differences, Arctic and Antarctic sea ice have different dominant properties and processes, which produce different responses to natural variability and anthropogenic change. In this Review, we survey changes in sea-ice properties and processes of the Arctic and Antarctic environments. In the Arctic, the melt season has lengthened by ~7.4 days decade −1 (1979-2024). There has been a corresponding decrease in winter sea-ice thickness (1.6 m mean total over 1980-2023), summer surface albedo (0.03 decade −1 over 1979-2020), and spring snow depth (2.5 cm decade −1 over 1954-2024). In comparison, corresponding changes in the Antarctic are small overall due to contrasting regional and temporal trends. From a dynamic perspective, sea-ice motion in both hemispheres has increased, with 0.63 cm s −1 decade −1 in the Arctic (1978-2024) and 0.69 cm s −1 decade −1 in the Antarctic (1982-2024), whereas trends in polynya occurrence differ by region. Future research should prioritize sea-ice observations that can inform how the climate system responds to natural variability and anthropogenic change, be standardized to increase uptake by the broader scientific community, and target process-oriented interactions to support model development.

Melinda A Webster

The impact of chemistry on anion migration in bixbyite-structured lanthanide oxides

This manuscript describes atomistic calculations of oxygen vacancy and interstitial migration in bixbyite structured lanthanide oxides. We examine two types of compounds, one in which only one type of lanthanide cation is present and a second class in which two lanthanides are present in a 3:1 ratio as dictated by the symmetry of the bixbyite crystal structure. Using temperature accelerated dynamics and the nudged elastic band method, we quantify the role of chemistry on the energy barriers for the most important pathways for both vacancy and interstitial migration. We then analyze the impact of these variations on the overall diffusivity of each defect, quantifying the contribution of each pathway using the theory of kinosons. We find that vacancy mobility can vary by as much as three orders of magnitude through changes in chemistry at 500 K. Changes in interstitial mobility are more modest but can still vary by an order of magnitude. This points to the ability to tune the mass transport characteristics of these compounds through appropriate choices in chemistry. We have also included supplementary information containing the atomic structures of the relevant pathways.

36 MATERIALS SCIENCE

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING

Interaction of magnetic fields with spinons in a fractionalized state

The 4d-electron trimer lattice Ba₄Nb₁₋ₓRu₃₊ₓO₁₂ exhibits either a quantum spin liquid (QSL) or a heavy-fermion strange metal (HFSM) phase, depending on Nb content. In the QSL state, itinerant spinons act as effective heat carriers, enhancing thermal conductivity. Strikingly, applying a magnetic field up to 14 T causes an abrupt, up-to-5000% increase in heat capacity below 150 mK, disrupting the linear temperature dependence typical of both phases. Meanwhile, AC susceptibility and electrical resistivity remain nearly unchanged, while thermal conductivity drops by up to 40% below 4 K. These results suggest spinons, despite being charge-neutral, are highly sensitive to magnetic fields at low temperatures. We propose that the magnetic field could induce Anderson localization of spinons, creating emergent non-magnetic two-level systems responsible for the Schottky-like anomaly in heat capacity. These findings point to a previously unexplored regime of spinon dynamics, potentially governed by field-induced localization and distinct from conventional magnetic or transport signatures.

Condensed-matter physics

Generalizable Web User Interface for Scalable and Streamlined Deployment of Building Energy Management Systems in Small and Medium-Sized Commercial Buildings

Small and medium-sized commercial buildings (SMCBs) comprise 94% of US commercial buildings yet face significant barriers to implementing building energy management systems despite advances in smart device technology. Existing solutions present critical limitations: cloud-based API solutions simplify deployment but create vendor lock-in constraints; commercial integrated software solutions ensure compatibility via standardized protocols but require substantial cost and technical expertise; open-source IoT platforms offer cost-effective vendor independence but provide insufficient standardized protocol support for commercial building automation. This research presents a generalizable web user interface framework that bridges the gap between evolving smart device capabilities and lagging software infrastructure for SMCBs. The proposed system integrates VOLTTRON open-source middleware with an automated configuration converter that transforms unified specifications written in YAML, a human-readable data-serialization format, into system-specific files, streamlining manual setup processes. The vendor-agnostic architecture supports industry-standard protocols (BACnet and Modbus) and semantic building models while providing adaptive web interfaces that dynamically adjust to various building configurations. Demonstrations through simulation-based testing and a field deployment show automatic interface adaptation across heterogeneous HVAC systems and multizone monitoring. The automated configuration converter also substantially reduces labor-intensive setup.

Chung, Jihoon [ORNL] (ORCID:0000000184880815)

UCB-GLOBES: An open-access mass spectral database of identified and unidentified atmospheric organic compounds

Chemical characterization of atmospheric organic aerosols using gas chromatography with 70 eV electron ionization mass spectrometry (GC/EI-MS) has been used for decades in advancing molecular marker detection and identification, though primarily through suspect screening and/or targeted analyses. To advance non-targeted analyses of environmental samples, we have catalogued approximately 27 000 mass spectra (MS) of the trimethylsilyl derivatives of semi-volatile organic aerosol (OA) analytes in the open-access University of California Berkeley Goldstein Library of Organic Biogenic Environmental Spectra (UCB-GLOBES). Analytes were observed in ambient samples from the U.S. and the Central Amazon and/or laboratory simulations of secondary OA (SOA) formation. These samples are representative of OA under urban and biomass burning influences as well as SOA derived from biogenic precursors (e.g., isoprene, monoterpenes, sesquiterpenes) and biomass burning intermediates. MS are documented in UCB-GLOBES without regard to known chemical identity, annotated with extensive metadata such as sample source/experimental conditions, any structural information gained from MS analyses, and predicted chemical properties such as average carbon oxidation state and carbon number. UCB-GLOBES MS are compatible for importing into the NIST MS Search program, and we have also provided a Jupyter Notebook for MS visualization and comparisons. We demonstrate the utility of UCB-GLOBES through MS reanalyses of prior analytes observed in ambient data, finding a 20 % reduction in the number of analytes assigned to OA source categories reliant solely on time series correlation and an overall 11 % increase in new MS-based OA source categorization for the Southeast U.S. For 1513 analytes observed previously in the Central Amazon, we found 375 MS matches using UCB-GLOBES vs. 136 MS matches during prior analyses, representing a 14 % gain in newly confirmed or newly categorized OA species. While OA from laboratory oxidation experiments in UCB-GLOBES are highly diverse chemically, on average only 29 % of UCB-GLOBES MS have a mass spectral match to another MS entry in UCB-GLOBES and/or in databases of known compounds (i.e. NIST MS Database, Adams Essential Oil, MANE Flavor and Fragrance Company). This indicates that roughly 70 % of UCB-GLOBES MS are unique thus far, not observed more than once among the laboratory oxidation samples and ambient data in UCB-GLOBES MS. Further, only 18 % can be positively identified using these databases or known authentic standards. This points to a large gap between these laboratory simulations and ambient OA. Overall, the UCB-GLOBES database can be utilized for improving confidence in OA source categorization and/or identification, novel chemical marker discovery, tracking chemical diversity, de novo structure and properties prediction, and improving MS search and matching algorithms. This can ultimately inform future research priorities for the chemical characterization of atmospheric organic samples.

Mass spectrometry

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown

Prime Time for Model-Predictive Control? Assessing the Technical and Market Readiness of Advanced Controls in Buildings

Despite three decades of extensive research and field testing that have consistently validated the benefits of Model Predictive Control (MPC) in building applications, the technology has seen limited market adoption. This paper evaluates the readiness of MPC for widespread deployment, showcases recent demonstrations and field tests across diverse building types, including residential, small commercial, large commercial, and campus settings. Our results demonstrate that MPC can optimize system operations to achieve load shifting, minimize curtailment of on-site generation, and reduce energy costs by up to 80 %, while maintaining or improving occupant comfort. We also show that MPC can effectively control large assets, such as MW-sized thermal storage systems, and respond to dynamic pricing signals. However, achieving scale remains difficult due to labor-intensive workflows, reliance on a “PhD-in-the-loop” for MPC design and maintenance, susceptibility to fragile data infrastructure, and persistent workforce education and acceptance barriers. To bridge this gap, we outline a transition from bespoke, labor intensive prototypes toward streamlined, segment-targeted deployment strategies that leverage model templates, semantic tools, and generative AI. By automating control configuration and reducing engineering effort, these recommendations provide a pathway for transforming successful research demonstrations into scalable, market ready solutions for MPC-based controls.

Pritoni, Marco

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)

Nanoscale Compositional and Strain Gradients Enable High‐Speed and Amplitude‐Resolved Pyroelectric Sensing

The frequency response of pyroelectric sensors is fundamentally governed by thermal time constant (τth, determined by thermal mass and thermal conductance) and electrical impedance arising from film capacitance and readout circuit. Conventional bulk LiTaO3 detectors are optimized for high responsivity at low modulation frequencies (0.1-10 Hz), possessing a large τth that thermally averages rapid temperature oscillations at elevated modulation frequencies, limiting fidelity in resolving dynamic varying thermal signals. Here, compositional and strain gradients are introduced into 100-nm-thick relaxor-ferroelectric films reducing τth to ≈2 µs and producing built-in potentials (≈1.45 V or 145 kV cm-1) that enhance the pyroelectric coefficient and suppress the dielectric constant. This enables complementary dual-mode operation by enhancing current-mode electrical responsivity and improving the voltage-mode figure of merit - advantageous for superior temperature resolution (ΔTmin ≈ 30 µK). The responsivity peak shifts to near 1 kHz (>2500-times higher than conventional bulk sensors), with measurable responsivity extending to a carrier frequency of 100 kHz and amplitude-resolved detection at modulation frequencies up to 15 kHz. These results establish nanoscale internal-field engineering can reshape electro-thermal trade-off in pyroelectric thin films toward zero-bias, high-thermal-sensitivity, and amplitude-resolved thermal sensing across a wide frequency bandwidth.

Lin, Ching‐Che

Optimizing Renewable Ammonia Production for a Sustainable Fertilizer Supply Chain Transition

Local renewable ammonia production using electrolytic hydrogen is an emerging approach to alleviate emissions attributed to synthetic nitrogen fertilizer production while also insulating against fluctuations in fertilizer prices and mitigating transportation costs and emissions. However, replacing ammonia currently produced using fossil fuels will not be immediate. To this end, we develop a supply chain transition model, which first optimizes the design and hourly operation of new renewable ammonia facilities to minimize production costs and then optimizes the annual installation timing, production scale, and location of these new renewable facilities along with ammonia transportation to meet county resolution demands. The objective is to augment and eventually replace conventional ammonia market imports in an economically competitive manner. We performed a case study for Minnesota's ammonia supply chain and found that a full transition to in-state renewable production by 2032 is optimal. This is incentivized by the U.S. federal government's clean hydrogen production credits. This transition results in 99 % reduction in carbon intensity along with stable supply costs below $475 per metric tonne. New renewable production facilities are an order of magnitude smaller than existing conventional plants. They use both wind and solar resources and operate dynamically to minimize expensive battery and hydrogen storage capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue

From PINNs to PIKANs: recent advances in physics-informed machine learning

Physics-Informed Neural Networks (PINNs) have emerged as a key tool in Scientific Machine Learning since their introduction in 2017, enabling the efficient solution of ordinary and partial differential equations using sparse measurements. Over the past few years, significant advancements have been made in the training and optimization of PINNs, covering aspects such as network architectures, adaptive refinement, domain decomposition, and the use of adaptive weights and activation functions. A notable recent development is the Physics-Informed Kolmogorov-Arnold Networks (PIKANS), which leverage a representation model originally proposed by Kolmogorov in 1957, offering a promising alternative to traditional PINNs. In this review, we provide a comprehensive overview of the latest advancements in PINNs, focusing on improvements in network design, feature expansion, optimization techniques, uncertainty quantification, and theoretical insights. We also survey key applications across a range of fields, including biomedicine, fluid and solid mechanics, geophysics, dynamical systems, heat transfer, chemical engineering, and beyond. Lastly, we review computational frameworks and software tools developed by both academia and industry to support PINN research and applications.

Kolmogorov-Arnold networks