4th Order Runge-Kutta Research Articles

On the dynamic characteristics of using track nonlinear energy sinks for structural vibration control

Nonlinear energy sinks (NESs) have emerged as a promising solution to overcome the narrow effective bandwidth of linear dynamic vibration absorbers such as the tuned mass damper (TMD). In particular, track NESs are increasingly attracting attention as they are capable of generating nonlinear restoring forces by the movement of an additional mass on the curved tracks. The present study aims to investigate and explore the versatile characteristics of track NESs in structural vibration control. To this end, three different track profiles are considered, namely the profiles designed by a quartic term only, by combined positive quadratic and quartic terms, and by combined negative quadratic and positive quartic terms. Their profiles and damping coefficients are analytically derived in closed form. The frequency response functions of standalone track NESs are constructed using the Harmonic Balance Method and compared with those obtained by the 4th-order Runge-Kutta method, and the dynamic responses, including the displacements, the phase planes, the nonlinear stiffness forces, and the vibrational frequencies, of the track NESs are presented and discussed under different excitation amplitudes and frequencies. The effectiveness and robustness of using track NESs for structural vibration control are also systematically investigated and compared with TMD under harmonic and white-noise ground excitations. Results show that the closed-form solutions of optimal profile and damping coefficients are sufficiently accurate for lightly damped structures, and the frequency response functions of track NESs can be accurately estimated using the Harmonic Balance Method for low and moderate excitations, although their dynamic responses are sensitive to the excitation amplitudes and frequencies. Furthermore, in comparison to TMD, the track NESs are superior in terms of decreased stiffness and smaller working strokes. In general, the present study demonstrates the potential of using track NESs as an effective and robust alternative for structural vibration control.

The Power Factor Improvement Modellingwith the Runge-Kutta 4thOrder, Analytical, and Experimental Method

The Power Factor Improvement Modellingwith the Runge-Kutta 4thOrder, Analytical, and Experimental Method

Mathematical Model on Distributed Denial of Service Attack in the Computer Network

In this paper, an electronic- epidemic two-folded mathematical model is formulated with help of non-linear ordinary differential equations. Distributed Denial of Service (DDoS) attacks in the computer network are studied. The modeling of both attacking nodes and targeting nodes is performed. Botnet based malicious devices and their threats on computer networks are addressed using appropriate parameters. The basic reproduction numbers for both the attacking and the targeting population are calculated and interpreted. Local and global stability analysis is carried out for the infection-free and endemic equilibrium points. Differential equations are solved with the help of the Runge-Kutta 4th order numerical method and graphs are analyzed using MATLAB software. Simulation shows that the success or failure depends on the number of initially infected computers in the attacking group. The proposed model exhibits the phenomenon of backward bifurcation for different values of transmission parameters. This model gives the theoretical base for controlling and predicting the DDoS attack. This shows the way to minimize the attack in the network. This study will be helpful to identify the botnet devices and run the latest version of antivirus in the network to protect against DDoS attacks from attacking sources. The application of this study is to ascertain online crime and locate the attacking nodes in the field of online transactions of real-life problems that involve the internet and computer networking systems. Moreover, our model can play an important role in policy-making against the distributed attack.

Entropy generation in free convective micropolar couple stress regime in vertical channel

The study reports an application of homotopy perturbation method to investigate entropy generation in the steady micropolar couple stress fluidics inside a vertical porous medium channel. The channel is composed of two isothermal parallel walls maintained at different temperatures. The governing equations are framed in the Cartesian system. The boundary value problem depicting the thermofluidic configuration is treated by the homotopy perturbation method (HPM) and for its validation, the 4th-order Runge-Kutta scheme is also invoked. The quantities of interest, viz. skin friction, wall couple stress, Nusselt number computed by HPM, and 4th order R-K method are found to be in good agreement. After getting the velocity, temperature, and microrotation regime, the entropy involved with the thermodynamic process was modeled to facilitate future thermodynamic optimizations.

Modeling of forced-vibration systems using continuous-time state-space neural network

Dynamic analysis of forced-vibration systems in civil engineering could be computationally inefficient or even hard to converge if the systems are stiff or highly complicated. Rapid advances in machine learning make it possible to formulate surrogate models for forced-vibration systems using neural networks. The widely used neural networks such as the convolutional neural network (CNN), recurrent neural network (RNN), etc., usually require a constant sampling rate and data length, thus they are difficult to be implemented for real-time calculation of the dynamic system with varying sampling rates. Recently, the continuous-time state-space neural network (CSNN) has shown the capability to lift these restrictions and has been drawing growing attention from the community. In this paper, we propose a generalized CSNN model for various forced-vibration systems (linear and nonlinear). The CSNN model comprises two sets of independent neural networks aiming to compute the state derivative and system response, respectively. Both neural networks adopt linear and nonlinear layers in parallel, instead of only fully connected nonlinear layers as adopted in the literature. This configuration is aimed to enhance the CSNN model with its capability to recognize the linear and nonlinear behaviors of systems. Additionally, the bias options in the CSNN model are all turned off to improve the stability of the model in the long-term time-series forecast, premised on the assumption that the forced-vibration systems are dissipative systems without drift, which is the most common case in civil engineering. Integration on the state derivative at the current time step is executed to obtain the state at the next time step using the explicit 4th-order Runge–Kutta method. Both numerical and experimental illustrative examples are provided, demonstrating that the CSNN model can achieve high performance and training efficiency with a few hyper-parameters, and thus is highly promising for engineering applications.

Practical considerations for HIL simulations of power converters using different numerical methods

The use of Hardware-In-the-Loop systems implemented in FPGAs is constantly growing due to their performance. However, hardware implementation of numerical methods to solve differential equations presents some challenges when applied to power converters. This paper shows a comparison of several numerical methods: Euler, Heun, Midpoint and 4th order Runge-Kutta, taking into account the accuracy of the methods and how they can be applied to power converters, where the equations change depending on the status of the switches. Results show that the speed in the solution of the numerical method is the main variable that affects the accuracy of the simulation, so in order to keep a fast resolution speed numerical method with low resources should be used.

Komparasi Skema Numerik Euler, Runge-Kutta dan Adam-Basforth-Moulton untuk Menyelesaikan Solusi Persamaan Osilator Harmonik

This article discusses the comparison of different numerical schemes to visualize the solution of 2nd-order differential equations. One-step methods such as the Euler method and the 4th-order Runge-Kutta method are combined with the 3rd-order Adam-Bashforth-Moulton method to solve the solution of 2nd-order differential equations. This combination of methods solves the Harmonic Oscillator equation, an 2nd-order differential equation widely applied in various oscillation contexts. The order of accuracy and order of approximation error are determined analytically. Finally, simulations are given with different steps for the three methods to confirm the behavior of the solution to the Harmonic Oscillator equation. The results show that the Euler method with the lowest order of accuracy has good accuracy at the beginning of the oscillation but not when time t is increased. The Runge-Kutta method, with the highest order of accuracy, shows excellent and consistent accuracy and solution stability, while the Adam-Bashforth-Moulton method, although it has a lower accuracy than the Runge-Kutta method of order 4, can be improved by choosing a one-step method with a high order of accuracy to approximate some of the required initial solutions. All three methods can provide approximation values with excellent accuracy and stability if a small step, h, is chosen, but this step can increase the time duration to display the solution. Thus, it is necessary to choose the right h according to the context of the equation and the method used to obtain accurate solutions with optimal time duration.

Propagation of magnetogasdynamics spherical shock wave in a gravitating gas with radiation heat flux

This article presents a mathematical model for characterizing the dynamic behavior of spherical shock waves in a self-gravitating, radiating ideal gas with the presence of an azimuthal magnetic field, emphasizing adiabatic conditions. The analysis assumes that the environment just ahead of the shock wave is stationary, and it accounts for variations in density, magnetic field, and fluid velocity within the disturbed medium just behind the shock front. Furthermore, the impact of thermal radiation within the context of an optically thin limit is incorporated into the energy equation of the governing system. Employing the Lie invariance method, the set of partial differential equations governing the flow within this medium is transformed into a system of nonlinear ordinary differential equations through the use of similarity variables. Four distinct cases of similarity solutions are derived by selecting different values for the arbitrary constants associated with the generators. Among these four cases, only two yield similarity solutions, one assuming a power-law shock path and the other an exponential-law shock path. In the case of a power-law shock path, the resulting set of nonlinear ordinary differential equations is numerically solved using the 4th-order Runge-Kutta method in MATLAB software. The article thoroughly discusses the influence of various parameters, including γ (adiabatic index of the gas), Ma−2 (Alfvén–Mach number), ϕ (ambient density exponent), and G0 (gravitational parameter), on the flow properties. The findings are presented graphically to provide a comprehensive understanding of the effects of these parameters.

Effects of Electrification and Transverse Force on Dusty Flow over a Linear Stretching Sheet

The numerical investigation of the flow and heat transfer of steady dusty flow over a linear stretching sheet has been carried out. The effects of Transverse force and electrification has been incorporated in this problem. Modelling of the problem comprises of highly nonlinear partial differential equations that have been transferred to systems of ordinary differential equations by implementing suitable transformations. Since the equations are of boundary value problems in nature, have been transferred to initial value problem by using shooting method and then solved by RungeKutta 4th order technique. All the above numerical methods incorporated in the BVP4C tool of MATLAB and has been solved the systems of differential equations by using MATLAB software. The effects of different parameter like Prandtl number, Eckert number, transverse force, and electrification parameter on flow and heat transfer profile has been presented through graphs and tables. The results have been validated with previous authors. From the graph it is observed that the transverse force reduces the velocity of the particle phase in the flow where as electrification of particles raises the temperature of particle phase in the flow. The computations carried out with tolerance of order less than O (10^(-6) ).

Growth kinetics of Bacillus cytotoxicus in liquid Egg yolk during treatment with phospholipase A2 – A one-step global dynamic analysis

During commercial production of liquid egg yolk (LEY), phospholipase A2 (PLA2) is used to improve its emulsification capacity and thermal stability. The enzymatic treatment may occur at elevated temperatures such as 50 °C, potentially allowing foodborne pathogens, such as Bacillus cereus, to grow. Little knowledge is available concerning growth of B. cereus in LEY during PLA2 treatment. Therefore, the objective of this study was to investigate the growth kinetics of B. cereus during PLA2 treatment using pathogenic B. cytotoxicus NVH391-98, the most thermotolerant member in the B. cereus group, as a surrogate.Inoculated LEY samples were placed in precision programmable incubators to observe the growth of B. cytotoxicus NVH391-98 under multiple isothermal and dynamic temperature conditions between 20 and 53 °C. The bacterial growth was described using the differential Baranyi model coupled with two different secondary models. The kinetic parameters were determined using one-step dynamic inverse analysis of multiple growth curves. The least square method was used in combination with the 4th order Runge-Kutta method to solve the differential Baranyi model using multiple growth curves to determine the cardinal kinetic parameters.The results showed that B. cytotoxicus NVH391-98 can grow prolifically at 50 °C. The estimated minimum, optimum and maximum temperatures were 16.7 or 18.5, 47.8 or 48.1, and 52.1 or 52.4 °C, respectively, depending on the secondary models, with an optimum growth rate of 2.1 log colony-forming-unit (CFU)/g per hour. The dynamic model is validated using isothermal curves with reasonable accuracy. B. cytotoxicus died off slowly at 15 °C. At 55 °C, thermal inactivation was observed, with a D value of approximately 2.7 h. Holding at 55 °C or below 15 °C can effectively prevent the growth of B. cytotoxicus in egg yolk.

Heat transfer rate and thermal energy analysis of MHD powell-eyring fluid in a permeable medium

This study is of significant importance as it focuses on exploring the intricate dynamics of viscous dissipation over a layered stretching surface in the context of MHD free and forced convective flow influenced by Powell-Eyring fluid. The primary aim of this research is to analyze the impact of heat emission and immersion phenomena on the flow, and to achieve this, conformable transformations of similarity have been pragmatic to adapt the governing equations from a set of non-linear PDE's into ODE's. The research methodology employed in this study involves numerical analysis using a 4th order Runge-Kutta formulation, implemented through MATLAB. The study precisely examines the stimulus of different flow constraints on velocity and temperature profiles. The outcomes of this investigation are presented through graphical representations, providing a comprehensive view of the flow characteristics. Furthermore, this study goes beyond the basic analysis by presenting important engineering metrics, the Nusselt number and skin-friction coefficient, explicitly across dissimilar variables. In conclusion, the findings of this research reveal as the Eyring-Powell fluid velocity profile rises, the temperature profile decreases, while the temperature profile rises with rising radiation.

Wave energy conversion using heaving oscillator inside ship: Conceptual design, mathematical model and parametric study

A concept of utilizing heaving oscillators inside a ship as a means of capturing wave energy is proposed, drawing inspiration from the two-body self-contained wave energy point absorber. This concept involves the ship acting as a small-amplitude-motion platform and the heaving oscillator as a large-stroke body, where the kinetic energy of their relative motion can be converted into electrical energy through a power take-off (PTO) system. The dynamic equations of the ship and the inboard heaving oscillator under regular waves are derived using the boundary element method and solved using a Python code with 4th-order Runge-Kutta integration. The numerical results exhibit favorable agreement with publicly available experimental data. Furthermore, the effects of forward speed, encounter angle, and PTO mechanical parameters on various aspects, including ship motion, oscillator response, wave energy capture performance, and inner loads, are investigated. The results show that the concept of heaving oscillator inside ship has no significantly adverse effects on the ship seakeeping performance. Both forward speed and encounter angle significantly influence the encounter period and thereby impact the sensitivity of capture width (CW) at a broader band of incident wave periods. Conversely, PTO mechanical parameters exerting a more localized impact primarily under long incident waves.

Inclined magnetic force impact on cross nanoliquid flowing with widening shallow and heat generating by using artificial neural network (ANN)

Artificial neural networks (ANNs) have a wide range of applications in science and technology. Artificial neural networks (ANNs) widely utilized in image and speech recognition, neural language processing, drug discovery, genomics and bioinformatics, Robotics and control systems, material science and energy and power system. Due to the above applications, the main focus of this article is to scrutinize the impact of inclined magnetic field on a cross nanofluid flow with slip velocity and convective boundary conditions over a variable porosity. Gold (Au) nanoparticles are suspended in base liquid blood. Thermal stratification and heat generation with joule heating impact is taken in the energy equation in order to see the heat transfer. To transform the non-linear partial differential equations (PDE's) into non-linear ordinary differential equations (ODE's) utilized the von-Karman similarity transformation parameters. For the solutions of the non-linear ODE's, the 4th- order Runge-Kutta numerical method is employed. A dataset for the employed neural network back propagated Levenberg-Marquard scheme (NN-BLMS) is generated for various estimations of the embedded parameters like Weissenberg number, magnetic parameter, angle of inclination, porosity and permeability parameters, thermal stratification parameter by utilizing the 4th-order Runge-Kutta numerical scheme. The testing, validation and training methods of NN-BLMS are utilized to scrutinize the approximate solution of cross nanofluid with slip velocity and convective boundary conditions. The performance of the proposed NN-BLMS to successfully solved the cross nanofluid model is endorsed via mean squared error, error histogram and regression analysis. Streamlines and heatlines are plotted for different parameters including in velocity and temperature equations.

Application of homotopy perturbation method to solve a nonlinear mathematical model of depletion of forest resources

Reduction in forest resources due to increasing global warming and population growth is a critical situation the World faces today. As these reserves decrease, it alarms new challenges that require urgent attention. In this paper, we provide a semi-analytical solution to a nonlinear mathematical model that studies the depletion of forest resources due to population growth and its pressure. With the help of the homotopy perturbation method (HPM), we determine an approximate series solution with few perturbation terms, which is one of the essential power of the HPM method. We compare our semi-analytical results with numerical solutions obtained using the Runge-Kutta 4th-order (RK-4) method. Furthermore, we analyze the model’s behaviour and dynamics by changing the parametric coefficients that represent the depletion rate of forest resources and the growth rate of population pressure and present these findings using various graphs.

MHD rotating flow over a stretching surface: The role of viscosity and aggregation of nanoparticles

The magnetohydrodynamic (MHD) rotating flow that occurs across a stretching surface has numerous practical applications in a variety of domains. These fields include astronomy, engineering, the material sciences, and space exploration. The combined examination of magnetohydrodynamics rotating flow across a stretching surface, taking into consideration fluctuating viscosity and nanoparticle aggregation, has significant ramifications across several different domains. It is essential for both the growth of technology and the attainment of deeper insights into the complicated fluid dynamics to maintain research in this field. Given the aforementioned motivation, the principal aim of this study is to examine the effects of variable viscosity on the bidirectional rotating magnetohydrodynamic flow over a stretching surface. Aggregation effects on nanoparticles are used in the analysis. Titania (TiO2) is taken nanoparticle and ethylene glycol as base fluid. The nonlinear ordinary differential equations and the boundary conditions that correspond to them can be transformed into a dimensionless form by using a technique called similarity transformation. To get a numerical solution to the transformed equation, the Runge-Kutta 4th order (RK-4) method is utilized, and this is done in conjunction with the shooting method. The impact of various leading variables on dimensionless velocity, the coefficients of temperature, skin friction and local Nusselt number are graphically represented. Velocity profiles in both direction increases with increasing values of φ. The Nusselt number increases with increasing values of the radiation and temperature ratio parameters. When a 1 % volume fraction of nanoparticles is introduced, the Nusselt number exhibits a 0.174 % increase for the aggregation model compared to the regular fluid in the absence of radiation effects. When the aggregation model is used with a 1 % volume fraction of nanoparticles, the skin friction increases by 0.1153 % in the x direction and by 0.1165 % in the y direction compared to the regular fluid. Tables show the variation in Nusselt numbers, as well as a comparison of the effects of nanoparticle's aggregation model without and with radiation. Moreover, the numerical results obtained were compared with previously published data, demonstrating a satisfactory agreement. We firmly believe that this finding will have extensive implications for engineering and various industries.

Analysis of Nonlinear Flow and Heat Transfer Issues in Nanofluid by Using RK-4th Order Method

Nanofluids, composed of a base fluid with dispersed nanoparticles, have gained significant attention due to their unique thermal and flow characteristics. Understanding the nonlinear flow and heat transfer phenomena in nanofluids is crucial for optimizing their performance in engineering applications. In this study, we propose using the Runge-Kutta 4th order (RK-4) numerical method to analyze these issues. The RK-4 method is known for its accuracy and efficiency in solving ordinary differential equations. By applying this method to the governing equations of fluid flow and heat transfer, we can accurately capture the nonlinear behavior of nanofluids. We investigate various nonlinear phenomena, such as boundary layer separation, flow instability, and heat transfer enhancement, considering factors like nanoparticle concentration, size, shape, fluid properties, Reynolds number, and temperature gradient. The RK-4 method provides accurate solutions, revealing critical parameters that affect flow and heat transfer behavior. Comparisons with experimental data and other numerical methods validate the approach, demonstrating its reliability. Parametric studies are conducted to optimize nanofluid characteristics for improved heat transfer and fluid flow. This research enhances our understanding of nonlinear flow and heat transfer in nanofluids and enables efficient numerical analysis, opening new avenues for improving nanofluid-based systems.

Mathematical study of viscoelastic polymer during roll-over-web coating

Roll-coating, also known as roll-to-roll coating or web coating is a continuous and high-speed method used in various industrial applications to apply a thin and uniform layer of material onto a substrate. This process is commonly employed in industries such as electronics, photovoltaic, printing, packaging, and other fields where precise and consistent coating is essential. Here we present an asymptotic analysis of forward roll coating over a rigid substrate for low Weissenberg number viscoelastic fluids. The analysis uses lubrication theory coupled with regular perturbation theory assuming a simplified Phan-Thien-Tanner (SPPT) constitutive equation and small values of the extensibility parameter ε. The analytical expressions of velocity fields, pressure gradient and pressure fields are obtained while the load carrying force, and power input are estimated by numerical solutions (4th order Runge-Kutta). Increases in both Wi and ε increases the pressure gradient in the nip region and increases in viscoelasticity reduce the overall force on the roll. Further, SPTT model predicts 14% lower pressure for [Formula: see text] in the nip region when compared to Newtonian model.

Bio-Convective Hartmann Flow of Couple Stress Hybrid Nanofluid Between Two Dilating Parallel Walls with Heat and Mass Transfer

The study of bio-convective flow of hybrid nanofluid attracted many researchers because of tremendous applications in the fields of biofuel biotechnology, enzyme-based biosensors and biomedical science. The present work addresses a comparative study of CuO/Al2O3-water and CuO-water nanoparticles on heat and mass transfer characteristics of the squeezing flow of MHD couple stress fluid between two parallel plates by suspending motile micro-organisms. An approximated numerical technique (Shooting method along with Runge-Kutta 4th order scheme) have been employed to analyse the system of coupled nonlinear ordinary differential equations. The above numerical investigations were carried out for various governing parameters such as couple stress parameter, Hartmann number, bioconvection Peclet number, squeezing parameter etc. The effects of these physical parameters are illustrated graphically over velocity components, temperature distribution, diffusion of concentration and density of motile microorganisms. In addition to this the numerical values of skin friction, the local Nusselt number and local Sherwood number are tabulated at the upper plate for CuO-water and CuO–Al2O3-water at the expanding and squeezing cases. The numerical results for temperature profiles are in good consistency with earlier research.

Microrotating chemically reactive hybrid nanomaterial in a high porous medium influenced by Cattaneo-Christov double diffusion and non-linear slip

ABSTRACT A simultaneous heat and mass transfer due to microrotating Darcy-Forchheimer flow of hybrid nanofluid over a moving thin needle is investigated. Darcy-Forchheimer medium accommodating hybrid nanofluid flow yields greater heat transfer rate, thereby leading to greater mass transfer rate over thin needle in industrial applications such as blood flow problems, aerodynamics, transportation, coating of wires, lubrication, and geothermal power generation. The thermophoresis and Brownian motion phenomena are introduced to enrich thermal treatment. Heat and mass transfer are accompanied by Cattaneo-Christov heat and mass flux. The hybrid nanofluid is radiative and dissipative in nature. Arrhenius pre-exponential factor law is introduced. Entropy generation analysis is carried out. The 4th order Runge-Kutta method along with shooting technique is devised to get requisite numerical solution of the transformed non-dimensional system of equations. Darcy-Forchheimer effect to simultaneous heat and mass transfer of microrotating hybrid nanofluid flow over thin needle subject to non-linear slip is the novelty of present study which is beyond of previous investigations. Rise in Forchheimer number (strengthening Darcy Forchheimer medium) leads to surface viscous drag decreases by 11.11% for hybrid nanofluid and 10.78% for pure nanofluid indicating the control of momentum transfer, thereby regulating heat transfer rate effectively.

LibFastMesh: An optimized finite-volume framework for computational aeroacoustics

This paper presents an optimized framework for simulating aeroacoustic flow on unstructured meshes. The libFastMesh (LFM) framework is based on a low-dissipation finite–volume discretization, together with high-order explicit time integration, for direct computation of aeroacoustic fields. These methods were previously implemented in caafoam, an OpenFOAM-based solver developed by D'Alessandro et al. [1] to resolve aeroacoustic flows. The new framework is specifically developed from scratch to alleviate two of the main bottlenecks currently limiting the performance of OpenFOAM-based solvers at large scale: memory bandwidth and inter-node communication. Compact data structures and compute kernel fusion are employed to enhance cache utilization and re-use, respectively. Separate treatment of inter-process boundary cells, together with non-blocking MPI communication enable effective overlap between communication and flux computations. The accuracy and robustness of libFastMesh are established via simulations of well-established aeroacoustic flow benchmarks. Comparison to previous simulation results obtained with caafoam, and to DNS data serve to validate the code. Finally, the effectiveness of the aforementioned optimizations is demonstrated through rigorous analysis of the proposed solver performance in single-node and multi-node operation modes. The obtained results show that libFastMesh offers a speed-up of up to 20x with respect to the previously developed OpenFOAM-based solver, caafoam, in large-scale computations. Moreover, LFM is shown to scale extremely for cell-per-core ratios of less than 500. These results are particularly appealing given the highly resource-demanding nature of direct computational aeroacoustics (CAA). Consequently, LFM could be an extremely attractive tool for scientists who wish to conduct large-scale CAA simulations on modern, exascale architectures. Program summaryProgram Title: libFastMeshCPC Library link to program files:https://doi.org/10.17632/7w9fy2xtcf.1Developer's repository link:https://github.com/TRC-HPC/LFM_PublicLicensing provisions: GPLv3Programming language: C++Nature of problem: This software solves compressible Navier–Stokes equations. The adopted numerics and flow models allow to capture with high accuracy aerodinamically generated sound. The solver is also designed to take advantage of the massively parallel computing systems which are strictly required for facing real–life aeroacoustic problems.Solution method: The libFastMesh (LFM) framework is based on a low-dissipation finite–volume discretization, together with high-order explicit time integration, for direct computation of aeroacoustic fields on unstructured meshes. Convective terms can be approximated using either the Kurganov–Noelle–Petrova (KNP) scheme or Pirozzoli's energy conserving approach. Time integration is fully–explicit and relies on a low-storage, 4th order Runge–Kutta scheme. The solver is specifically developed to run efficiently on modern many-core CPUs, and at large-scale where it exhibits excellent strong scaling features. This is made possible thanks to several optimizations that alleviate memory bandwidth and inter-node communication bottlenecks.Additional comments: OpenFOAM-v2112+ is an installation prerequisite. Standard OpenFOAM pre-processing and post-processing tools, as well as mesh decomposition and reordering methods may be used with the solver.