time step size
Recently Published Documents

Tank sloshing is widely present in many engineering fields, especially in the field of marine. Due to the trend of large-scale liquid cargo ships, it is of great significance to study the coupled motion response of ships with tanks in beam waves. In this study, the CFD (Computational Fluid Dynamics) method and experiments are used to study the response of a ship with/without a tank in beam waves. All the computations are performed by an in-house CFD solver, which is used to solve RANS (Reynold Average Navier-Stokes) equations coupled with six degrees-of-freedom solid-body motion equations. The Level Set Method is used to solve the free surface. Verification work on the grid number and time step size has been conducted. The simulation results agree with the experimental results well, which shows that the numerical method is accurate enough. In this paper, several different working conditions are set up, and the effects of the liquid height in the tank, the size of the tank and the wavelength ratio of the incident wave on the ship’s motion are studied. The results show the effect of tank sloshing on the ship’s motion in different working conditions.

AbstractThe sensitivity of moving particle semi-implicit (MPS) simulations to numerical parameters is investigated in this study. Although the verification and validation (V&V) are important to ensure accurate numerical results, the MPS has poor performance in convergences with a time step size. Therefore, users of the MPS need to tune numerical parameters to fit results into benchmarks. However, such tuning parameters are not always valid for other simulations. We propose a practical numerical condition for the MPS simulation of a two-dimensional wedge slamming problem (i.e., an MPS-slamming condition). The MPS-slamming condition is represented by an MPS-slamming number, which provides the optimum time step size once the MPS-slamming number, slamming velocity, deadrise angle of the wedge, and particle size are decided. The simulation study shows that the MPS results can be characterized by the proposed MPS-slamming condition, and the use of the same MPS-slamming number provides a similar flow.

AbstractIn this paper, we develop an optimized hybrid block method which is combined with a modified cubic B-spline method, for solving non-linear partial differential equations. In particular, it will be applied for solving three well-known problems, namely, the Burgers equation, Buckmaster equation and FitzHugh–Nagumo equation. Most of the developed methods in the literature for non-linear partial differential equations have not focused on optimizing the time step-size and a very small value must be considered to get accurate approximations. The motivation behind the development of this work is to overcome this trade-off up to much extent using a larger time step-size without compromising accuracy. The optimized hybrid block method considered is proved to be A-stable and convergent. Furthermore, the obtained numerical approximations have been compared with exact and numerical solutions available in the literature and found to be adequate. In particular, without using quasilinearization or filtering techniques, the results for small viscosity coefficient for Burgers equation are found to be accurate. We have found that the combination of the two considered methods is computationally efficient for solving non-linear PDEs.

This paper presents computational fluid dynamics (CFD) simulations of the flow around a horizontal axis hydrokinetic turbine (HAHT) found in the literature. The volume of fluid (VOF) model implemented in a commercial CFD package (ANSYS-Fluent) is used to track the air-water interface. The URANS SST k-ω and the four-equation Transition SST turbulence models are employed to compute the unsteady three-dimensional flow field. The sliding mesh technique is used to rotate the subdomain that includes the turbine rotor. The effect of grid resolution, time-step size, and turbulence model on the computed performance coefficients is analyzed in detail, and the results are compared against experimental data at various tip speed ratios (TSRs). Simulation results at the analyzed rotor immersions confirm that the power and thrust coefficients decrease when the rotor is closer to the free surface. The combined effect of rotor and support structure on the free surface evolution and downstream velocities is also studied. The results show that a maximum velocity deficit is found in the near wake region above the rotor centerline. A slow wake recovery is also observed at the shallow rotor immersion due to the free-surface proximity, which in turn reduces the power extraction.

We consider the long-time behavior of an explicit tamed exponential Euler scheme applied to a class of parabolic semilinear stochastic partial differential equations driven by additive noise, under a one-sided Lipschitz continuity condition. The setting encompasses nonlinearities with polynomial growth. First, we prove that moment bounds for the numerical scheme hold, with at most polynomial dependence with respect to the time horizon. Second, we apply this result to obtain error estimates, in the weak sense, in terms of the time-step size and of the time horizon, to quantify the error to approximate averages with respect to the invariant distribution of the continuous-time process. We justify the efficiency of using the explicit tamed exponential Euler scheme to approximate the invariant distribution, since the computational cost does not suffer from the at most polynomial growth of the moment bounds. To the best of our knowledge, this is the first result in the literature concerning the approximation of the invariant distribution for SPDEs with non-globally Lipschitz coefficients using an explicit tamed scheme.

A numerical formulation based on the precise-integration time-domain (PITD) method for simulating periodic media is extended for overcoming the Courant-Friedrich-Levy (CFL) limit on the time-step size in a finite-difference time-domain (FDTD) simulation. In this new method, the periodic boundary conditions are implemented, permitting the simulation of a wide range of periodic optical media, i.e., gratings, or thin-film filters. Furthermore, the complete tensorial derivation for the permittivity also allows simulating anisotropic periodic media. Numerical results demonstrate that PITD is reliable and even considering anisotropic media can be competitive compared to traditional FDTD solutions. Furthermore, the maximum allowable time-step size has been demonstrated to be much larger than that of the CFL limit of the FDTD method, being a valuable tool in cases in which the steady-state requires a large number of time-steps.

AbstractIn this paper, we will consider a geometrically exact Cosserat beam model taking into account the industrial challenges. The beam is represented by a framed curve, which we parametrize in the configuration space $\mathbb{S}^{3}\ltimes \mathbb{R}^{3}$ S 3 ⋉ R 3 with semi-direct product Lie group structure, where $\mathbb{S}^{3}$ S 3 is the set of unit quaternions. Velocities and angular velocities with respect to the body-fixed frame are given as the velocity vector of the configuration. We introduce internal constraints, where the rigid cross sections have to remain perpendicular to the center line to reduce the full Cosserat beam model to a Kirchhoff beam model. We derive the equations of motion by Hamilton’s principle with an augmented Lagrangian. In order to fully discretize the beam model in space and time, we only consider piecewise interpolated configurations in the variational principle. This leads, after approximating the action integral with second order, to the discrete equations of motion. Here, it is notable that we allow the Lagrange multipliers to be discontinuous in time in order to respect the derivatives of the constraint equations, also known as hidden constraints. In the last part, we will test our numerical scheme on two benchmark problems that show that there is no shear locking observable in the discretized beam model and that the errors are observed to decrease with second order with the spatial step size and the time step size.

Fuel starvation is a major cause of anode corrosion in low temperature polymer electrolyte fuel cells. The fuel cell start-up is a critical step, as hydrogen may not yet be evenly distributed in the active area, leading to local starvation. The present work investigates the hydrogen distribution and risk for starvation during start-up and after nitrogen purge by extending an existing computational fluid dynamic model to capture transient behavior. The results of the numerical model are compared with detailed experimental analysis on a 25 cm2 triple serpentine flow field with good agreement in all aspects and a required time step size of 1 s. This is two to three orders of magnitude larger than the time steps used by other works, resulting in reasonably quick calculation times (e.g., 3 min calculation time for 1 s of experimental testing time using a 2 million element mesh).

Purpose A transient magneto-quasistatic vector potential formulation involving nonlinear material is spatially discretized using the finite element method of first and second polynomial order. By applying a generalized Schur complement the resulting system of differential algebraic equations is reformulated into a system of ordinary differential equations (ODE). The ODE system is integrated in time by using explicit time integration schemes. The purpose of this paper is to investigate explicit time integration for eddy current problems with respect to the performance of the first-order explicit Euler scheme and the Runge-Kutta-Chebyshev (RKC) method of higher order. Design/methodology/approach The ODE system is integrated in time using the explicit Euler scheme, which is conditionally stable by a maximum time step size. To overcome this limit, an explicit multistage RKC time integration method of higher order is used to enlarge the maximum stable time step size. Both time integration methods are compared regarding the overall computational effort. Findings The numerical simulations show that a finer spatial discretization forces smaller time step sizes. In comparison to the explicit Euler time integration scheme, the multistage RKC method provides larger stable time step sizes to diminish the overall computation time. Originality/value The explicit time integration of the Schur complement vector potential formulation of eddy current problems is accelerated by a multistage RKC method.