scholarly journals Investigation of the ship–seabed interaction with a high-fidelity CFD approach

Author(s):  
Weichao Shi ◽  
Mingxin Li ◽  
Zhiming Yuan

AbstractShip–seabed interaction is highly critical to ship safety and ship performance when ship operates over shallow water with uneven seabed. However, the prediction method for ship reaction is still yet to be fully developed. In this paper, a high-fidelity transient numerical simulation method with sliding mesh is developed based on computational fluid dynamics to model a vessel passing a step bank, demonstrating to be a computational cost economical solution. The comprehensive numerical model is validated and verified against benchmarked experimental and numerical studies, which is proved to be highly accurate in predicting force characteristics and wave development. Meanwhile, during the research the impulse effect generated by the step bank was found to have striking effects on the wave elevation, wake development and vessel sinkage. Five regions on the behaviour of vessel sinkage are defined in the present work. According to the result, the vessel will encounter the extreme sinkage after a relatively long distance (12 $${L}_{\mathrm{pp}}$$ L pp at $${F}_{h2}=0.519$$ F h 2 = 0.519 ) passing the step bank instead of immediately.

2020 ◽  
Author(s):  
Ali Raza ◽  
Arni Sturluson ◽  
Cory Simon ◽  
Xiaoli Fern

Virtual screenings can accelerate and reduce the cost of discovering metal-organic frameworks (MOFs) for their applications in gas storage, separation, and sensing. In molecular simulations of gas adsorption/diffusion in MOFs, the adsorbate-MOF electrostatic interaction is typically modeled by placing partial point charges on the atoms of the MOF. For the virtual screening of large libraries of MOFs, it is critical to develop computationally inexpensive methods to assign atomic partial charges to MOFs that accurately reproduce the electrostatic potential in their pores. Herein, we design and train a message passing neural network (MPNN) to predict the atomic partial charges on MOFs under a charge neutral constraint. A set of ca. 2,250 MOFs labeled with high-fidelity partial charges, derived from periodic electronic structure calculations, serves as training examples. In an end-to-end manner, from charge-labeled crystal graphs representing MOFs, our MPNN machine-learns features of the local bonding environments of the atoms and learns to predict partial atomic charges from these features. Our trained MPNN assigns high-fidelity partial point charges to MOFs with orders of magnitude lower computational cost than electronic structure calculations. To enhance the accuracy of virtual screenings of large libraries of MOFs for their adsorption-based applications, we make our trained MPNN model and MPNN-charge-assigned computation-ready, experimental MOF structures publicly available.<br>


2021 ◽  
Vol 11 (13) ◽  
pp. 5900
Author(s):  
Yohei Fujinami ◽  
Pongsathorn Raksincharoensak ◽  
Shunsaku Arita ◽  
Rei Kato

Advanced driver assistance systems (ADAS) for crash avoidance, when making a right-turn in left-hand traffic or left-turn in right-hand traffic, are expected to further reduce the number of traffic accidents caused by automobiles. Accurate future trajectory prediction of an ego vehicle for risk prediction is important to activate the assistance system correctly. Our objectives are to propose a trajectory prediction method for ADAS for safe intersection turnings and to evaluate the effectiveness of the proposed prediction method. Our proposed curve generation method is capable of generating a smooth curve without discontinuities in the curvature. By incorporating the curve generation method into the vehicle trajectory prediction, the proposed method could simulate the actual driving path of human drivers at a low computational cost. The curve would be required to define positions, angles, and curvatures at its initial and terminal points. Driving experiments conducted at real city traffic intersections proved that the proposed method could predict the trajectory with a high degree of accuracy for various shapes and sizes of the intersections. This paper also describes a method to determine the terminal conditions of the curve generation method from intersection features. We set a hypothesis where the conditions can be defined individually from intersection geometry. From the hypothesis, a formula to determine the parameter was derived empirically from the driving experiments. Public road driving experiments indicated that the parameters for the trajectory prediction could be appropriately estimated by the obtained empirical formula.


Author(s):  
Wei Zhang ◽  
Saad Ahmed ◽  
Jonathan Hong ◽  
Zoubeida Ounaies ◽  
Mary Frecker

Different types of active materials have been used to actuate origami-inspired self-folding structures. To model the highly nonlinear deformation and material responses, as well as the coupled field equations and boundary conditions of such structures, high-fidelity models such as finite element (FE) models are needed but usually computationally expensive, which makes optimization intractable. In this paper, a computationally efficient two-stage optimization framework is developed as a systematic method for the multi-objective designs of such multifield self-folding structures where the deformations are concentrated in crease-like areas, active and passive materials are assumed to behave linearly, and low- and high-fidelity models of the structures can be developed. In Stage 1, low-fidelity models are used to determine the topology of the structure. At the end of Stage 1, a distance measure [Formula: see text] is applied as the metric to determine the best design, which then serves as the baseline design in Stage 2. In Stage 2, designs are further optimized from the baseline design with greatly reduced computing time compared to a full FEA-based topology optimization. The design framework is first described in a general formulation. To demonstrate its efficacy, this framework is implemented in two case studies, namely, a three-finger soft gripper actuated using a PVDF-based terpolymer, and a 3D multifield example actuated using both the terpolymer and a magneto-active elastomer, where the key steps are elaborated in detail, including the variable filter, metrics to select the best design, determination of design domains, and material conversion methods from low- to high-fidelity models. In this paper, analytical models and rigid body dynamic models are developed as the low-fidelity models for the terpolymer- and MAE-based actuations, respectively, and the FE model of the MAE-based actuation is generalized from previous work. Additional generalizable techniques to further reduce the computational cost are elaborated. As a result, designs with better overall performance than the baseline design were achieved at the end of Stage 2 with computing times of 15 days for the gripper and 9 days for the multifield example, which would rather be over 3 and 2 months for full FEA-based optimizations, respectively. Tradeoffs between the competing design objectives were achieved. In both case studies, the efficacy and computational efficiency of the two-stage optimization framework are successfully demonstrated.


Author(s):  
Benjamin D. Youngman ◽  
David B. Stephenson

We develop a statistical framework for simulating natural hazard events that combines extreme value theory and geostatistics. Robust generalized additive model forms represent generalized Pareto marginal distribution parameters while a Student’s t -process captures spatial dependence and gives a continuous-space framework for natural hazard event simulations. Efficiency of the simulation method allows many years of data (typically over 10 000) to be obtained at relatively little computational cost. This makes the model viable for forming the hazard module of a catastrophe model. We illustrate the framework by simulating maximum wind gusts for European windstorms, which are found to have realistic marginal and spatial properties, and validate well against wind gust measurements.


2013 ◽  
Vol 353-356 ◽  
pp. 436-439
Author(s):  
De Sen Kong ◽  
Yong Po Chen

In order to forecast the stability of deep roadway and optimize the parameters of bolts, the complex stress environment and the multivariate surrounding rocks characteristics of deep roadway were analyzed. Then the classification prediction method and the numerical simulation method were simultaneously used to analysis the stability of surrounding rocks. Furthermore, the supporting parameters of bolts were also designed optimally. It was shown that the characteristics of stress distribution, deformation and failure zone of surrounding rocks are not ideal. So it is necessary to optimize the supporting parameters of deep roadway. All these research findings will provide the theory basis for bolts of deep roadway and will ensure the optimization of bolts and the stability of deep roadway in the long run.


2021 ◽  
Author(s):  
Ruizi Zhang ◽  
Ian Frigaard

Abstract Many numerical studies have been conducted regarding laminar miscible displacement flow in narrow, vertical, eccentric annuli. For the next decade it is likely that primary cementing flows on the scale of the well will continue to be simulated predominantly with 2D gap-averaged (2DGA) models. However, 3D simulations are less common due to the computational cost. The comparison between 2D and 3D models needs further attention, to understand the main discrepancies and thus help to understand primary cementing flows better. In this paper, comparisons of 3D against 2DGA model results show a range of interesting different phenomena, e.g. static layers, dispersive spikes, and instabilities. The predictions of the 2DGA model are the same as the 3D results to a degree. In particular, they are consistent with each other regarding the evolving process, interface shape, etc. However, the main difference with the 2DGA concentration arises from dispersion on the scale of the annular gap. From the recent research of Renteria and Frigaard (J. Fluid Mech., vol. 905, 2020) [1], a variety of dispersive effects are the main discrepancy between experiments and 2DGA results as well. We give representative examples of these flows in surface casing geometries and suggest methods for improvement of the 2DGA model.


2020 ◽  
Author(s):  
Shine Win Naung ◽  
Mohammad Rahmati ◽  
Hamed Farokhi

Abstract The high-fidelity computational fluid dynamics (CFD) simulations of a complete wind turbine model usually require significant computational resources. It will require much more resources if the fluid-structure interactions between the blade and the flow are considered, and it has been the major challenge in the industry. The aeromechanical analysis of a complete wind turbine model using a high-fidelity CFD method is discussed in this paper. The distinctiveness of this paper is the application of the nonlinear frequency domain solution method to analyse the forced response and flutter instability of the blade as well as to investigate the unsteady flow field across the wind turbine rotor and the tower. This method also enables the aeromechanical simulations of wind turbines for various inter blade phase angles in a combination with a phase shift solution method. Extensive validations of the nonlinear frequency domain solution method against the conventional time domain solution method reveal that the proposed frequency domain solution method can reduce the computational cost by one to two orders of magnitude.


2021 ◽  
Author(s):  
Robert Elian Feteanu

Experimental and numerical studies have been undertaken to examine various aspects pertaining to the interaction of an incident travelling shock wave with a solid rocket motor's head end (forward section), in order to identify any potential gasdynamic mechanism of wave reinforcement pertinent to combustion instability behaviour in these motors. A cold-flow experiment, based on a shock tube scheme tailored to the present application, has proved to be useful in providing information surrounding the interaction process. Both experimental and numerical results (CFD simulations) confirm the existence of substantial transient radial wave development superimposed on the base reflected axial shock wave. These results illustrate the potential weakness of one-dimensional flow models for certain engineering applications, where important multidimensional phenomena, such as those observed in this work, may not be captured. By analogy to actual propulsion system combustion chambers, the transverse wave activity is potentially a factor in supporting an augmentation of the local combustion rate in the head-end region of a rocket motor combustor.


2021 ◽  
Author(s):  
Anuj Dhoj Thapa

Gillespie's algorithm, also known as the Stochastic Simulation Algorithm (SSA), is an exact simulation method for the Chemical Master Equation model of well-stirred biochemical systems. However, this method is computationally intensive when some fast reactions are present in the system. The tau-leap scheme developed by Gillespie can speed up the stochastic simulation of these biochemically reacting systems with negligible loss in accuracy. A number of tau-leaping methods were proposed, including the explicit tau-leaping and the implicit tau-leaping strategies. Nonetheless, these schemes have low order of accuracy. In this thesis, we investigate tau-leap strategies which achieve high accuracy at reduced computational cost. These strategies are tested on several biochemical systems of practical interest.


Author(s):  
Sriram Shankaran ◽  
Brian Barr

The objective of this study is to develop and assess a gradient-based algorithm that efficiently traverses the Pareto front for multi-objective problems. We use high-fidelity, computationally intensive simulation tools (for eg: Computational Fluid Dynamics (CFD) and Finite Element (FE) structural analysis) for function and gradient evaluations. The use of evolutionary algorithms with these high-fidelity simulation tools results in prohibitive computational costs. Hence, in this study we use an alternate gradient-based approach. We first outline an algorithm that can be proven to recover Pareto fronts. The performance of this algorithm is then tested on three academic problems: a convex front with uniform spacing of Pareto points, a convex front with non-uniform spacing and a concave front. The algorithm is shown to be able to retrieve the Pareto front in all three cases hence overcoming a common deficiency in gradient-based methods that use the idea of scalarization. Then the algorithm is applied to a practical problem in concurrent design for aerodynamic and structural performance of an axial turbine blade. For this problem, with 5 design variables, and for 10 points to approximate the front, the computational cost of the gradient-based method was roughly the same as that of a method that builds the front from a sampling approach. However, as the sampling approach involves building a surrogate model to identify the Pareto front, there is the possibility that validation of this predicted front with CFD and FE analysis results in a different location of the “Pareto” points. This can be avoided with the gradient-based method. Additionally, as the number of design variables increases and/or the number of required points on the Pareto front is reduced, the computational cost favors the gradient-based approach.


Sign in / Sign up

Export Citation Format

Share Document