Explicit Dynamics Research Articles

Numerical analysis of initial cloud of powder dynamic dispersion

Numerical simulations were performed to investigate the explosion dispersion characteristics of powder materials with respect to the implicated velocity, and a fluid–solid coupling simulation model of the explosion dispersion process of the initial cloud was established. Explicit dynamics models of Finite-element method and smooth particle hydrodynamics for the detonation driving stage of central charge and the cloud isentropic expansion stages were developed, the pressure loads of the explosive-driven were generated using the user-defined function program, and computational fluid dynamics models for the isentropic expansion and free expansion stages of the initial cloud were developed. Additionally, the maximum expansion radius of the cloud in explosion dispersion tests and theoretical calculations were compared with the radius in the numerical simulation results to verify the accuracy of the simulation model. Analyses of the morphology, concentration, and turbulence field in the numerical simulation revealed that vortex stratification occurred in the “two wings” of the initial cloud in its development process and that a particle-free “hollow” area is subsequently formed. Analyses also revealed that the vortex stratification occurs earlier as the implicated velocity of powder material increases. Additionally, an increase in the implicated velocity can effectively weaken the concentration of cloud and make the distribution of cloud particles more uniform.

Mapping function model representing the relationship between explosion shock wave pressure and wavefront temperature

The evaluation of ammunition damage power and guidance for ammunition design heavily relies on the shock wave pressure and wavefront temperature produced by an ammunition explosion. However, temperature test results are often inaccurate and unreliable. Therefore, this study utilized Autodyn explicit dynamics simulation software to conduct finite element numerical simulations of explosion shock wave pressure, wavefront temperature propagation, and distribution for trinitrotoluene explosives weighing 10, 20, 50, and 100 kg. The shock wave pressure and wavefront temperature were measured at different measuring points. The cloud maps of wavefront propagation evolution obtained at different explosion times were analyzed to determine the attenuation laws of pressure and temperature propagation in the near field and far field. Based on the similarity law of explosions and the dimensional analysis method, a mapping function model was established to represent the relationship between a shock wave’s peak pressure and peak temperature. The parameters of the model include explosive mass (w), measuring point radius (r), specific heat capacity in the air (c), and the peak pressure of an explosion shock wave (p). The model’s accuracy in calculating the explosion shock wavefront temperature exceeds 89.75%, effectively resolving the issue of low accuracy in the results of explosion field temperature tests and calculations. Therefore, this study provides data and theoretical support for testing and evaluating the damage power of ammunition explosives, and the proposed model has high applicability in the engineering field.

Evaluation of Accelerator Pedal Strength under Critical Loads Using the Finite Element Method

The core idea of the research consists in a formulation of boundary conditions of a mechanical accelerator pedal’s strength in an Ansys environment, whose conditions are equivalent to full-scale tests under the critical loads defined by the UNECE’s Regulation No. 13. The lack of regulatory requirements for the strength of pedal types other than brake pedals is a major gap in vehicle certification, especially when it comes to agricultural machinery. In such cases, the authors suggest being guided by UNECE R 13 regarding the strength of the accelerator and other types of pedals and checking their behavior under loads of at least 1000 N. The real value of the yield strength of the material (Silumin 4000) is very important, both in the physical real-life experiments and in FEA simulation. The critical case of a short-term shock loading of the pedal in its extreme position has been considered separately. With the help of the Ansys Explicit Dynamics module, results of a pedal’s behavior were obtained; it lost its integrity and suffered destruction. It is also necessary to check the intermediate stress values depending on the loads for direct and hybrid tasks using the Transient Structural module in order to estimate other critical cases of the pedal behavior.

Shock absorption capability of corrugated ring yield mount subjected to high impact loading

Various types of mechanical energy-absorbing devices are known that operate by plastic deformation. The corrugated ring mount that is used in this study relates to a device that absorbs energy by plastic deformation. This energy-absorbing device has reduced volumetric proportions, simple in design, and therefore has small overall dimensions and can be mass-produced at low cost. This study aims to determine the shock absorption capability and efficiency of this mount against impact loading. For this, Finite Element Method Analysis (FEA) and experimentation are done. The FEA is done using the Explicit Dynamics (AutoDyn) module of ANSYS Workbench and for experimentation Drop Test Machine (DTM) is used. In this study impact load from low g up to 85 g is applied and a very close agreement is found between FEA and experimental results. There is just a 5–10% deviation between the findings. The results show that this mount is plastically deformed to absorb the impact energy with a maximum efficiency of 70%. It concludes that it is a reliable and safer shock energy device.

Sparse identification of Lagrangian for nonlinear dynamical systems via proximal gradient method

The autonomous distillation of physical laws only from data is of great interest in many scientific fields. Data-driven modeling frameworks that adopt sparse regression techniques, such as sparse identification of nonlinear dynamics (SINDy) and its modifications, are developed to resolve difficulties in extracting underlying dynamics from experimental data. However, SINDy faces certain difficulties when the dynamics contain rational functions. The Lagrangian is substantially more concise than the actual equations of motion, especially for complex systems, and it does not usually contain rational functions for mechanical systems. Few proposed methods proposed to date, such as Lagrangian-SINDy we have proposed recently, can extract the true form of the Lagrangian of dynamical systems from data; however, these methods are easily affected by noise as a fact. In this study, we developed an extended version of Lagrangian-SINDy (xL-SINDy) to obtain the Lagrangian of dynamical systems from noisy measurement data. We incorporated the concept of SINDy and used the proximal gradient method to obtain sparse Lagrangian expressions. Further, we demonstrated the effectiveness of xL-SINDy against different noise levels using four mechanical systems. In addition, we compared its performance with SINDy-PI (parallel, implicit) which is a latest robust variant of SINDy that can handle implicit dynamics and rational nonlinearities. The experimental results reveal that xL-SINDy is much more robust than the existing methods for extracting the governing equations of nonlinear mechanical systems from data with noise. We believe this contribution is significant toward noise-tolerant computational method for explicit dynamics law extraction from data.

Biomechanics of calcaneus impacted by talus: a dynamic finite element analysis

This paper aimed to investigate the biomechanical changes during the talus impact with the calcaneus at varying velocities. Various three-dimensional reconstruction software was utilized to construct a finite element model that consisted of the talus, calcaneus, and ligaments. The explicit dynamics method was used to explore the process of the talus impacting on the calcaneus. The velocity of impact was altered from 5 m/s to 10 m/s with a 1 m/s interval. Stress readings were collected from the posterior, intermediate, and anterior subtalar articular (PSA, ISA, ASA), calcaneocubic articular (CA), Gissane Angle (GA), calcaneal base (BC), medial wall (MW), and lateral wall (LW) of the calcaneus. The changes in the amount and distribution of stress in the different regions of the calcaneus that varied with velocity were analysed. The model was validated through comparison with findings from the existing literature. During the process of impact between the talus and calcaneus, the stress in the PSA reached its peak first. Notably, stress was concentrated mainly in the PSA, ASA, MW, and LW of the calcaneus. At varying impact velocities of the talus, the mean maximum stress of the PSA, LW, CA, BA, and MW exhibited statistically significant differences (P values were 0.024, 0.004, <0.001, <0.001, and 0.001, respectively). However, the mean maximum stress of the ISA, ASA, and GA was not statistically significant (P values were 0.289, 0.213, and 0.087, respectively). In comparison with the velocity at 5 m/s, the mean maximum stress increases in each region of the calcaneus at a velocity of 10 m/s were as follows: PSA 73.81%, ISA 7.11%, ASA 63.57%, GA 89.10%, LW 140.16%, CA 140.58%, BC 137.67%, MW 135.99%. The regions of stress concentration were altered, and the magnitude and sequence of peak stress in the calcaneus also varied according to the velocity of the talus during impact. In conclusion, the velocity of the talus during impact had a significant influence on the magnitude and distribution of stress within the calcaneus, which was a crucial factor in the development of calcaneal fractures. It was possible that the magnitude and sequence of stress peaks played a vital role in determining the emergence of fracture patterns.

Critical time-step size analysis and mass scaling by ghost-penalty for immersogeometric explicit dynamics

In this article, we study the effect of small-cut elements on the critical time-step size in an immersogeometric explicit dynamics context. We analyze different formulations for second-order (membrane) and fourth-order (shell-type) equations, and derive scaling relations between the critical time-step size and the cut-element size for various types of cuts. In particular, we focus on different approaches for the weak imposition of Dirichlet conditions: by penalty enforcement and with Nitsche’s method. The conventional stability requirement for Nitsche’s method necessitates either a cut-size dependent penalty parameter, or an additional ghost-penalty stabilization term. Our findings show that both techniques suffer from cut-size dependent critical time-step sizes, but the addition of a ghost-penalty term to the mass matrix serves to mitigate this issue. We confirm that this form of ‘mass-scaling’ does not adversely affect error and convergence characteristics for a transient membrane example, and has the potential to increase the critical time-step size by orders of magnitude. Finally, for a prototypical simulation of a Kirchhoff–Love shell, our stabilized Nitsche formulation reduces the solution error by well over an order of magnitude compared to a penalty formulation at equal time-step size.

Multiscale analysis of wheel-rail rolling contact wear and damage mechanisms using molecular dynamics and explicit finite elements

As traffic density and axle loads increase, the macroscopic wheel-rail contact can lead to severe rail wear. This paper constructs a multi-scale analysis framework from macroscopic force characteristics to microscopic wear mechanisms based on explicit finite element methods and molecular dynamics models. The study shows that the mutual proximity and annihilation of reverse dislocations in the rail material under the action of wheel-rail forces is an important cause of grain boundary fragmentation, while the proliferation of dislocations on the Subsurface leads to work hardening of the rail. The analytical framework proposed in this paper achieves a multiscale correlation between wheel-rail macroscopic contact and material microscopic response, and the conclusions obtained reveal the microscopic wear mechanism of rail materials.

Effect of debonding damage on transient dynamic response on stern ramp door of passenger ship 6311 GT

The sandwich, as an advanced material, increases innovation in the field of ship construction. One of the applications is on the ramp door structure. The characteristics of sandwich material which are light and strong, are the main reasons for the application of this material. However, sandwich material also has drawbacks, one of which is debonding damage, where the bond between the face and core is released. This damage can affect the structure’s stiffness, which impacts reducing the strength of the structure on the ramp door. The effect of debonding is identified by performing numerical tests using finite element software. This test aims to identify the ramp door’s dynamic response due to the debonding damage’s influence. This research uses finite element software with an explicit dynamics method to determine the dynamic response of stress and displacement due to loading on the ramp door. The results of testing with finite element software show that an increase in the ratio of debonding damage causes the amplitude of the dynamic response to increase when subjected to transient loading.

English

Recently, there has been a high interest in using lightweight aluminum foams for automotive, railway and aerospace operations. Because of its high ductility and deformability, Aluminum foam is generally used for energy absorption for crashworthiness applications. To keep safe and avoid occupant injuries, it is necessary to absorb the kinetic energy generated during impact. Thus, to absorb high kinetic energy, the crash box material needs a special material microstructure, which is light in weight and can absorb further energy than the being one like CaCo3, CBC, or SiC. B4C etc. In particular, the analysis of energy absorption of aluminum foam in automotive for energy absorption applications is limited. The main objective of this exploration is to analyze and optimize the porosity size and voids percentage on impact energy absorption of aluminum foam using a numerical approach. For this purpose, first, fifteen CAD Al foam samples were developed by using Digimat multi-scale material modeling software. Second, cubic elements with circular bubble shape at 5, 10 and 15 void percentage and at 1.5 mm, 2 mm,2.5 mm, 3 mm and 3.5 mm bubble sizes were modeled. Finally, the numerical analysis of impact energy by using ANSYS workbench19.2 Explicit dynamics by applying initial low velocity was performed. The parameters were compared to each other to optimize the proper percentage composition and cell size for the best energy absorption applications. The effects of bubble shape, foaming agent and percentage composition on energy immersion were discussed. In this study, the analysis was fulfilled by determining and comparing the energy absorptions of all the models and, finally, comparing them with existing foaming agents.

Gross-Neveu model with O(2)L×O(2)R chiral symmetry: Duality with Zakharov-Mikhailov model and large N solution

The two-flavor Gross-Neveu model with ${\mathrm{U}(2)}_{L}\ifmmode\times\else\texttimes\fi{}\mathrm{U}(2{)}_{R}$ chiral symmetry in $1+1$ dimensions is used to construct a novel variant of four-fermion theories with ${\mathrm{O}(2)}_{L}\ifmmode\times\else\texttimes\fi{}\mathrm{O}(2{)}_{R}$ chiral symmetry. The spontaneous breaking of the group O(2), a continuous group with two connected components (rotations and reflections), gives rise to new phenomena. It is ideally suited to describe a situation where two distinct kinds of condensation compete, in particular chiral symmetry breaking (particle-hole condensation) and Cooper pairing (particle-particle condensation). After solving the O(2) chiral Gross-Neveu model in detail, we demonstrate that it is dual to another classically integrable model due to Zakharov and Mikhailov. The duality enables us to solve the quantum version of this model in the large $N$ limit with semiclassical methods, supporting its integrability at the quantum level. The resulting model is the unique four-fermion theory sharing the full Pauli-G\"ursey symmetry with free, massless fermions (``perfect Gross-Neveu model'') and provides us with a solvable model for competing chiral and Cooper pair condensates, including explicit soliton dynamics and the phase diagram.

Design and Simulation of Banana Pseudo-stem Fibre Extracting Raspador (Cylinder) Using FEA Technique

ABSTRACT Development of mechanical fiber extractors without considering engineering knowledge can cause failures of the machine’s components during operations. Therefore, the present study was conducted to examine the structural strength of the banana pseudo-stem fiber-extracting raspador using finite element analysis. For measuring the force required to extract the fibers, a measuring setup along with an embedded system was developed using the ATmega328P microcontroller, DYMH-103 load cell, and 24-Bit ATD converter, and it was found that an average of 99.09 N force was required to scratch the fibers from single pseudo-stem sheath. Three computer-aided design (CAD) models of raspador having tip angles of blades of 25°, 35°, and 45° were developed in the Creo-Parametric software and simulated in the Ansys® software using three material properties, i.e. high-carbon steel, medium-carbon steel, and low-carbon steel. The results of both static structural analysis and explicit dynamics simulation analysis suggested that the CAD model with 45° tip angle of blade having low-carbon steel has a higher value of total deformation and maximum principal stress and a comparatively lower value of equivalent elastic strain than other models and under the yield limit of the parent material. Hence, it may be used for the development of a prototype of a banana pseudo-stem fiber-extracting machine.

3D hydrodynamics simulations of internal gravity waves in red giant branch stars

ABSTRACT We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the ppmstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a $1.2\, {\rm M}_{\odot }$ upper RGB star. We perform simulations for different grid resolutions (7683, 15363, and 28803), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial set-up increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\, H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers.

Gamma-Hemolysin Components: Computational Strategies for LukF-Hlg2 Dimer Reconstruction on a Model Membrane.

The gamma-hemolysin protein is one of the most common pore-forming toxins expressed by the pathogenic bacterium Staphylococcus aureus. The toxin is used by the pathogen to escape the immune system of the host organism, by assembling into octameric transmembrane pores on the surface of the target immune cell and leading to its death by leakage or apoptosis. Despite the high potential risks associated with Staphylococcus aureus infections and the urgent need for new treatments, several aspects of the pore-formation process from gamma-hemolysin are still unclear. These include the identification of the interactions between the individual monomers that lead to the formation of a dimer on the cell membrane, which represents the unit for further oligomerization. Here, we employed a combination of all-atom explicit solvent molecular dynamics simulations and protein-protein docking to determine the stabilizing contacts that guide the formation of a functional dimer. The simulations and the molecular modeling reveal the importance of the flexibility of specific protein domains, in particular the N-terminus, to drive the formation of the correct dimerization interface through functional contacts between the monomers. The results obtained are compared with the experimental data available in the literature.

SlicerCBM: automatic framework for biomechanical analysis of the brain

PurposeBrain shift that occurs during neurosurgery disturbs the brain’s anatomy. Prediction of the brain shift is essential for accurate localisation of the surgical target. Biomechanical models have been envisaged as a possible tool for such predictions. In this study, we created a framework to automate the workflow for predicting intra-operative brain deformations.MethodsWe created our framework by uniquely combining our meshless total Lagrangian explicit dynamics (MTLED) algorithm for computing soft tissue deformations, open-source software libraries and built-in functions within 3D Slicer, an open-source software package widely used for medical research. Our framework generates the biomechanical brain model from the pre-operative MRI, computes brain deformation using MTLED and outputs results in the form of predicted warped intra-operative MRI.ResultsOur framework is used to solve three different neurosurgical brain shift scenarios: craniotomy, tumour resection and electrode placement. We evaluated our framework using nine patients. The average time to construct a patient-specific brain biomechanical model was 3 min, and that to compute deformations ranged from 13 to 23 min. We performed a qualitative evaluation by comparing our predicted intra-operative MRI with the actual intra-operative MRI. For quantitative evaluation, we computed Hausdorff distances between predicted and actual intra-operative ventricle surfaces. For patients with craniotomy and tumour resection, approximately 95% of the nodes on the ventricle surfaces are within two times the original in-plane resolution of the actual surface determined from the intra-operative MRI.ConclusionOur framework provides a broader application of existing solution methods not only in research but also in clinics. We successfully demonstrated the application of our framework by predicting intra-operative deformations in nine patients undergoing neurosurgical procedures.

The influence of explicit local dynamics on range expansions driven by long-range dispersal.

Range expansions are common in natural populations. They can take such forms as an invasive species spreading into a new habitat or a virus spreading from host to host during a pandemic. When the expanding species is capable of dispersing offspring over long distances, population growth is driven by rare but consequential long-range dispersal events that seed satellite colonies far from the densely occupied core of the population. These satellites accelerate growth by accessing unoccupied territory, and also act as reservoirs for maintaining neutral genetic variation present in the originating population, which would ordinarily be lost to drift. Prior theoretical studies of dispersal-driven expansions have shown that the sequential establishment of satellites causes initial genetic diversity to be either lost or maintained to a level determined by the breadth of the distribution of dispersal distances. If the tail of the distribution falls off faster than a critical threshold, diversity is steadily eroded over time; by contrast, broader distributions with a slower falloff allow some initial diversity to be maintained for arbitrarily long times. However, these studies used lattice-based models and assumed an instantaneous saturation of the local carrying capacity after the arrival of a founder. Real-world populations expand in continuous space with complex local dynamics, which potentially allow multiple pioneers to arrive and establish within the same local region. Here, we evaluate the impact of local dynamics on the population growth and the evolution of neutral diversity using a computational model of range expansions with long-range dispersal in continuous space, with explicit local dynamics that can be controlled by altering the mix of local and long-range dispersal events. We found that many qualitative features of population growth and neutral genetic diversity observed in lattice-based models are preserved under more complex local dynamics, but quantitative aspects such as the rate of population growth, the level of maintained diversity, and the rate of decay of diversity all depend strongly on the local dynamics. Besides identifying situations in which modeling the explicit local population dynamics becomes necessary to understand the population structure of jump-driven range expansions, our results show that local dynamics affects different features of the population in distinct ways, and can be more or less consequential depending on the degree and form of long-range dispersal as well as the scale at which the population structure is measured.

Modeling of wound ballistics in biological tissues simulators

Modern weapons cause severe damage, accompanied by high rates of complications and mortality. The investigation of such kinds of weapons is in high demand considering the ongoing active phase of russia's war against Ukraine since February 2022. In order to understand the pathological processes that occur in and outside the gunshot wound, we conducted an experimental study using mathematical simulation. The results presented in the article will help to choose the appropriate surgical management and improve the results of treatment. The aim of this study was to investigate and evaluate the damaging effect of a 5.45 mm 7N6M bullet and a 5.45 mm V-max expansive bullet using numerical modeling of wound canals in ballistic plasticine. The Ansys Explicit Dynamics engineering complex was used to simulate the dynamics of the bullet’s motion. The basic equations, solved by the explicit dynamic analysis, express the conservation of mass, momentum, and energy in Lagrange coordinates. Together with the material model and the set of initial and boundary conditions, they determine the complete solution to the problem. Taking into account that the initial velocity of the bullet is 1185 m/s at a mass of 3.9 g, we obtain energy 2740 J. All this energy acts at the area of the wound canal with a depth of 150 mm. Injury with a conventional 7N6M bullet is characterized by the fact that it passes through the block and loses only part of the kinetic energy. The simulation results showed that the velocity of the bullet at the outlet is 220 m/s. Taking into account the initial velocity of the bullet 918 m/s with a mass of 3.4 g, we obtain the kinetic energy acting on the walls of the wound canal with a depth of 200 mm of about 830 J. Mathematic analyses showed that the expansive bullet has a soft core that deforms and transfers all the kinetic energy to the tissues immediately after penetration into the tissues. The loss of kinetic energy of the bullet (ΔE, J) is defined as the difference between the kinetic energy at the time of injury (Ec, J) and the residual energy of the bullet when leaving the material (Er, J). Numerical modeling of wound ballistics in biological tissue simulators allows us to determine with high accuracy the features of wound canal formation and tissue response to damage of bullets having different kinetic energy, which contributes to the choice of adequate surgical management during surgery for gunshot wounds.

Effective Molecular Dynamics from Neural Network-Based Structure Prediction Models.

Recent breakthroughs in neural network-based structure prediction methods, such as AlphaFold2 and RoseTTAFold, have dramatically improved the quality of computational protein structure prediction. These models also provide statistical confidence scores that can estimate uncertainties in the predicted structures, but it remains unclear to what extent these scores are related to the intrinsic conformational dynamics of proteins. Here, we compare AlphaFold2 prediction scores with explicit large-scale molecular dynamics simulations of 28 one- and two-domain proteins with varying degrees of flexibility. We demonstrate a strong correlation between the statistical prediction scores and the explicit motion derived from extensive atomistic molecular dynamics simulations and further derive an elastic network model based on the statistical scores of AlphFold2 (AF-ENM), which we benchmark in combination with coarse-grained molecular dynamics simulations. We show that our AF-ENM method reproduces the global protein dynamics with improved accuracy, providing a powerful way to derive effective molecular dynamics using neural network-based structure prediction models.

Crashworthiness Study of a Newly Developed Civil Aircraft Fuselage Section with Auxiliary Fuel Tank Reinforced with Composite Foam

Over the past two decades, aircraft crashworthiness has seen major developments, mainly with modern computing systems and commercial finite element (FE) codes. The structure and the material have been designed to absorb more kinetic energy to ensure enough safety during a controlled crash condition. However, the fuselage section with an onboard auxiliary fuel tank requires special arrangements, since the inclined strut system with an efficient energy absorber is difficult to install under the cabin floor due to the space occupied by the fuel tank. To solve this shortcoming, a PVC composite foam along with an aluminum plate is introduced beneath the fuel tank to improve the crashworthiness metrics of the fuselage. Drop tests for both the conventional design and the proposed model are investigated by adopting the nonlinear explicit dynamics code Ansys Autodyn, with an impact velocity of 9.14 m/s. It was found that the kinetic energy absorption of the original fuselage section can be improved by 3.54% by reinforcing the foam and the plate. Moreover, they contribute to 20% of total internal energy dissipation. Numerical outcomes also suggest that the cabin floor surface experiences a 41% lower maximum stress, in addition to the mitigation of the maximum peak acceleration responses of the cabin floor at different measured locations from 6% to 36%.

3D hydrodynamics simulations of core convection in supermassive main-sequence stars

ABSTRACT Supermassive stars are Population III stars with masses exceeding $10^4\, {\rm M}_{\odot }$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code ppmstar to simulate the hydrogen-burning core of a $10^4\, {\rm M}_{\odot }$ supermassive main-sequence star. These are the first three-dimensional hydrodynamics simulations of core convection in supermassive stars. We perform a series of 10 simulations at different heating rates and on Cartesian grids with resolutions of 7683, 11523, and 17283. We examine different properties of the convective flow, including its large-scale morphology, its velocity spectrum, and its mixing properties. We conclude that the radiation pressure-dominated nature of the interior does not noticeably affect the behaviour of convection compared to the case of core convection in a massive main-sequence star where gas pressure dominates. Our simulations also offer support for the use of mixing-length theory in one-dimensional models of supermassive stars.