Dynamic flight stability of hovering model insects:theory versus simulation using equations of motion coupled with Navier-Stokes equations

In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model （which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion, so that the flapping wings＇ degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments）; the simulation solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theory and the simulation provides a test to the validity of the assumptions in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency （164 Hz） and the other is a model hawkmoth which has relatively low wingbeat frequency （26 Hz）. The results show that the averaged model is valid for the hawkmoth as well as for the dronefly. Since the wingbeat frequency of the hawkmoth is relatively low （the characteristic times of the natural modes of motion of the body divided by wingbeat period are relatively large） compared with many other insects, that the theory based on the averaged model is valid for the hawkmoth means that it could be valid for many insects.

Special Lie symmetry and Hojman conserved quantity of Appell equations in a dynamical system of relative motion

Special Lie symmetry and the Hojman conserved quantity for Appell equations in a dynamical system of relative motion are investigated. The definition and the criterion of the special Lie symmetry of Appell equations in a dynamical system of relative motion under infinitesimal group transformation are presented. The expression of the equation for the special Lie symmetry of Appell equations and the Hojman conserved quantity, deduced directly from the special Lie symmetry in a dynamical system of relative motion, are obtained. An example is given to illustrate the application of the results.

A model for predicting SST in limited region-I. The dynamical equations

-Starting from physical oceanology characteristics of the China seas and for the short-term operational prediction of SST in the region, a two-dimensional (vertically integrated) primitive equation model, physically reasonable and operationally feasible,on the upper mixed layer is constructed and given here, which consists of three parts, the nondivergent residual current (the monthly mean field of the Kuroshio and its branches) equations, the dynamic forecasting equations, and the equation of model's physics consisting of surface heat flux, coolings of the upper mixed layer due to the Ekman pumping and the entrainment by gale. This model may be used primarily to forecast the sea surface temperature, and to give estimations of the mean wind-driven current and the sea level, for a period of 3-5 d. In part 1 of this series, the physical conditions for establishing model equations are discussed first, that is, 1. the existence of the upper well mixed layer in the region; 2. the distinguishability of currents of all kinds; 3. the splitting of thermodynamical equation. The equations of nondivergent residual current, and the dynamic forecasting equations with initial values and boundary conditions are also discussed.

Monolithic Coupling of the Pressure and Rigid Body Motion Equations in Computational Marine Hydrodynamics

In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume(FV) framework,where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach(i.e. "strongly coupled" approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.

Modeling Dynamic Systems by Using the Nonlinear Difference Equations Based on Genetic Programming

When acquaintances of a model are little or the model is too complicate to build by using traditional time series methods, it is convenient for us to take advantage of genetic programming (GP) to build the model. Considering the complexity of nonlinear dynamic systems, this paper proposes modeling dynamic systems by using the nonlinear difference e-quation based on GP technique. First it gives the method, criteria and evaluation of modeling. Then it describes the modeling algorithm using GP. Finally two typical examples of time series are used to perform the numerical experiments. The result shows that this algorithm can successfully establish the difference equation model of dynamic systems and its predictive result is also satisfactory.

MODELLING OF MECHANICAL MECHANISM OF CHROMATOGRAPHIC SYSTEM AND THEORETICAL EQUATIONS SHOWING DYNAMICAL CHARACTERISTICS OF CHROMATOGRAPHY

A simple model of chromatographic mechanical mechanism is present, and then a scrics of theoretical chromatographic equations and fundamental Formulae are derived. These theoretical equations and formulae not only reserve thermodynamic characteristics in the current fundamental chromatographic formulae, but also introduce one or more kinetic parameter, so it is possible to make the macroscopic-control on the effect of kinetic characteristics on chromatographic system.

A Comparative Study of Adomian Decomposition Method with Variational Iteration Method for Solving Linear and Nonlinear Differential Equations

This study compares the Adomian Decomposition Method (ADM) and the Variational Iteration Method (VIM) for solving nonlinear differential equations in engineering. Differential equations are essential for modeling dynamic systems in various disciplines, including biological processes, heat transfer, and control systems. This study addresses first, second, and third-order nonlinear differential equations using Mathematica for data generation and graphing. The ADM, developed by George Adomian, uses Adomian polynomials to handle nonlinear terms, which can be computationally intensive. In contrast, VIM, developed by He, directly iterates the correction functional, providing a more straightforward and efficient approach. This study highlights VIM’s rapid convergence and effectiveness of VIM, particularly for nonlinear problems, where it simplifies calculations and offers direct solutions without polynomial derivation. The results demonstrate VIM’s superior efficiency and rapid convergence of VIM compared with ADM. The VIM’s minimal computational requirements make it practical for real-time applications and complex system modeling. Our findings align with those of previous research, confirming VIM’s efficiency of VIM in various engineering applications. This study emphasizes the importance of selecting appropriate methods based on specific problem requirements. While ADM is valuable for certain nonlinearities, VIM’s approach is ideal for many engineering scenarios. Future research should explore broader applications and hybrid methods to enhance the solution’s accuracy and efficiency. This comprehensive comparison provides valuable guidance for selecting effective numerical methods for differential equations in engineering.

A stochastic two-dimensional intelligent driver car-following model with vehicular dynamics

The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic flow analysis,simulation,autonomous vehicle development,etc.Two-dimensional(2D)vehicular movement is basically stochastic and is the result of interactions between a driver's behavior and a vehicle's characteristics.Current microscopic models either neglect 2D noise,or overlook vehicle dynamics.The modeling capabilities,thus,are limited,so that stochastic lateral movement cannot be reproduced.The present research extends an intelligent driver model(IDM)by explicitly considering both vehicle dynamics and 2D noises to formulate a stochastic 2D IDM model,with vehicle dynamics based on the stochastic differential equation(SDE)theory.Control inputs from the vehicle include the steer rate and longitudinal acceleration,both of which are developed based on an idea from a traditional intelligent driver model.The stochastic stability condition is analyzed on the basis of Lyapunov theory.Numerical analysis is used to assess the two cases:(i)when a vehicle accelerates from a standstill and(ii)when a platoon of vehicles follow a leader with a stop-and-go speed profile,the formation of congestion and subsequent dispersion are simulated.The results show that the model can reproduce the stochastic 2D trajectories of the vehicle and the marginal distribution of lateral movement.The proposed model can be used in both a simulation platform and a behavioral analysis of a human driver in traffic flow.

The Application of Thermomechanical Dynamics (TMD) to the Analysis of Nonequilibrium Irreversible Motion and a Low-Temperature Stirling Engine

We applied the method of Thermomechanical Dynamics (TMD) to a low-temperature Stirling engine, and the dissipative equation of motion and time-evolving physical quantities are self-consistently calculated for the first time in this field. The thermomechanical states of the heat engine are in Nonequilibrium Irreversible States (NISs), and time-dependent thermodynamic work W(t), internal energy E(t), energy dissipation or entropy Qd(t), and temperature T(t), are precisely studied and computed in TMD. We also introduced the new formalism, Q(t)-picture of thermodynamic heat-energy flows, for consistent analyses of NISs. Thermal flows in a long-time uniform heat flow and in a short-time heat flow are numerically studied as examples. In addition to the analysis of time-dependent physical quantities, the TMD analysis suggests that the concept of force and acceleration in Newtonian mechanics should be modified. The acceleration is defined as a continuously differentiable function of Class C2 in Newtonian mechanics, but the thermomechanical dynamics demands piecewise continuity for acceleration and thermal force, required from physical reasons caused by frictional variations and thermal fluctuations. The acceleration has no direct physical meaning associated with force in TMD. The physical implications are fundamental for the concept of the macroscopic phenomena in NISs composed of systems in thermal and mechanical motion.

Dynamic Modeling and Three-Dimensional Motion Analysis of Underwater Gliders

For consideration of both the eccentric rotatable rigid body and the translational rigid body, the dynamic model of the underwater glider is derived. Dynamical behaviors are also studied based on the model and can be used as the guidance to underwater gliders design. Gibbs function of the underwater glider system is derived first, and then the nonlinear dynamic model is obtained by use of Appell equations. The relationships between dynamic behaviors and design parameters are studied by solving the dynamic model. The spiral motion, swerving motion in three dimensions and the saw-tooth motion of the underwater glider in vertical plane are studied. Lake trials are carried out to validate the dynamic model.

Dynamics model of underwater robot motion control in 6 degrees of freedom

In order to analyze underwater robot control system dynamics features, a system 6-DOF dynamics model was founded. Underwater robot linear and nonlinear hydrodynamics were analyzed by Taylor series, based on general motion equation. Special control system motion equation was deduced by cluster of inertial items and non-inertial items. For program convenience, motion equation matrix format was presented. Experimental principles of screw propellers, rudders and wings were discussed. Experimental data least-square curve fitting, interpolation and their corresponding traditional equation helped us to obtain the whole system dynamic response procedure. A series of simulation experiments show that the dynamics model is correct and reliable. The model can provide theory proof for analyzing underwater robot motion control system physics characters and provide a mathematic model for traditional control method.

An Integrated Dynamic Model of Ocean Mining System and Fast Simulation of Its Longitudinal Reciprocating Motion

An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.

Lie symmetry and its generation of conserved quantity of Appell equation in a dynamical system of the relative motion with Chetaev-type nonholonomic constraints

Lie symmetry and conserved quantity deduced from Lie symmetry of Appell equations in a dynamical system of relative motion with Chetaev-type nonholonomic constraints are studied.The differential equations of motion of the Appell equation for the system,the definition and criterion of Lie symmetry,the condition and the expression of generalized Hojman conserved quantity deduced from Lie symmetry for the system are obtained.The condition and the expression of Hojman conserved quantity deduced from special Lie symmetry for the system under invariable time are further obtained.An example is given to illustrate the application of the results.

Mei symmetry and Mei conserved quantity of the Appell equation in a dynamical system of relative motion with non-Chetaev nonholonomic constraints

The Mei symmetry and the Mei conserved quantity of Appell equations in a dynamical system of relative motion with non-Chetaev nonholonomic constraints are studied.The differential equations of motion of the Appell equation for the system,the definition and the criterion of the Mei symmetry,and the expression of the Mei conserved quantity deduced directly from the Mei symmetry for the system are obtained.An example is given to illustrate the application of the results.

Noether Symmetry and Noether Conserved Quantity of Nielsen Equation for Dynamical Systems of Relative Motion

Noether symmetry of Nielsen equation and Noether conserved quantity deduced directly from Noether symmetry for dynamical systems of the relative motion are studied. The definition and criteria of Noether symmetry of a Nielsen equation under the infinitesimal transformations of groups are given. Expression of Noether conserved quantity deduced directly from Noether symmetry of Nielsen equation for the system are obtained. Finally, an example is given to illustrate the application of the results.

CANONICAL FORMS OF NIELSEN’S AND CENOV’S DYNAMICAL EQUATIONS

In this paper,with Poincare's formalism,and an indirect method,the canonical forms of the generalized equations of motion due to Nielsen and Cenov of a holonomic dynamical system in the velocity-phase space and the acenleration-phase space are obtained in terms of the Poincare parameters.

COORDINATION OF THE MOTION PARAMETERS WITHIN STEP-BY-STEP INTEGRATION FOR DYNAMIC EQUATION

A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.

Linking Structural Equation Modeling with Bayesian Network and Its Application to Coastal Phytoplankton Dynamics in the Bohai Bay

Bayesian networks （BN） have many advantages over other methods in ecological modeling, and have become an increasingly popular modeling tool. However, BN are flawed in regard to building models based on inadequate existing knowledge. To overcome this limitation, we propose a new method that links BN with structural equation modeling （SEM）. In this method, SEM is used to improve the model structure for BN. This method was used to simulate coastal phytoplankton dynamics in the Bohai Bay. We demonstrate that this hybrid approach minimizes the need for expert elicitation, generates more reasonable structures for BN models, and increases the BN model＇s accuracy and reliability. These results suggest that the inclusion of SEM for testing and verifying the theoretical structure during the initial construction stage improves the effectiveness of BN models, especially for complex eco-environment systems. The results also demonstrate that in the Bohai Bay, while phytoplankton biomass has the greatest influence on phytoplankton dynamics, the impact of nutrients on phytoplankton dynamics is larger than the influence of the physical environment in summer. Furthermore, although the Redfield ratio indicates that phosphorus should be the primary nutrient limiting factor, our results show that silicate plays the most important role in regulating phytoplankton dynamics in the Bohai Bay.

Prediction of Aircraft＇s Longitudinal Motion Based on Aerodynamic Coefficients and Derivatives by Surrogate Model Approach

The accuracy of a flight simulation is highly dependent on the quality of the aerodynamic database and prediction accuracies of the aerodynamic coefficients and derivatives. A surrogate model is an approximation method that is used to predict unknown functions based on the sampling data obtained by the design of experiments. This model can also be used to predict aerodynamic coefficients/derivatives using several measured points. The objective of this paper is to develop an efficient digital flight simulation by solving the equation of motion to predict the aerodynamics data using a surrogate model. Accordingly, there is a need to construct and investigate aerodynamic databases and compare the accuracy of the surrogate model with the exact solution, and hence solve the equation of motion for the flight simulation analysis. In this study, sample datas for models are acquired from the USAF Stability and Control DATCOM, and a database is constructed for two input variables （the angle of attack and Mach number）, along with two derivatives of the X-force axis and three derivatives for the Z-force axis and pitching moment. Furthermore, a comparison of the value predicted by the Kriging model and the exact solution shows that its flight analysis prediction ability makes it possible to use the surrogate model in future analyses.

DYNAMIC MODELING FOR AIRSHIP EQUIPPEDWITH BALLONETS AND BALLAST

Total dynamics of an airship is modeled. The body of an airship is taken as a submerged rigid body with neutral buoyancy, i.e., buoyancy with value equal to that of gravity, and the coupled dynamics between the body with ballonets and ballast is considered. The total dynamics of the airship is firstly derived by Newton-Euler laws and Kirchhoff’s equations. Furthermore, by using Hamiltonian and Lagrangian semi-direct product reduction theories, the dynamics is formulated as a Lie-Poisson system, or also an Euler-Poincaré system. These two formulations can be exploited for the control design using energy-based methods for Hamiltonian or Lagrangian system.