Assistance #185
Amortissement visqueux en RD
Ajouté par Frank Petitjean il y a environ 6 ans.
Mis à jour il y a presque 6 ans.
Description
Gérard,
J'ai mis en place un amortissement mixte cinétique-visqueux dans mes calculs de ballon pour la prochaine version de Omher. En prenant la précaution de rester en visqueux dans les incréments qui suivent le 1er basculement la convergence est fortement améliorée (facteur 2 à 3 !).
Problème : les contraintes sont toujours nettement plus faibles qu'avec amortissement cinétique seul.
J'ai l'impression que l'amortissement est trop rapide (sur-critique ?). J'utilise les paramètres standards :
type_calcul_visqu_critique= 1 opt_cal_C_critique= 1 f= 0.9_
Question : quel autre jeu de paramètres utiliser pour relâcher un peu l'amortissement (sous-critique) ?
Merci
Frank
Fichiers
|
modele.CVisu (7,88 ko)
modele.CVisu |
|
Frank Petitjean, 16/10/2018 09:38
|
|
|
maillage.her (93,5 ko)
maillage.her |
|
Frank Petitjean, 16/10/2018 09:38
|
|
|
modele.info (6,04 ko)
modele.info |
|
Frank Petitjean, 16/10/2018 09:38
|
|
|
modele_mixte.info (6,04 ko)
modele_mixte.info |
|
Frank Petitjean, 16/10/2018 09:39
|
|
|
modele.info (6,04 ko)
modele.info |
|
Frank Petitjean, 16/10/2018 11:14
|
|
- Statut changé de Nouveau à En cours
Je ne comprend pas ce qui se passe, car normalement à convergence la partie visqueuse n'intervient pas dans l'équilibre.
Donc est-ce qu'il s'agit des contraintes obtenues pour un équilibre assez précis (du genre < 0.001) ?
Si oui, c'est anormal et je suis intéressé pour regarder plus en détail un cas test représentatif
Difficile d'extraire un cas test élémentaire significatif. Les écarts se produisent avec un maillage déjà un peu fin.
Voici le cas du BSO100Z mono-fuseau avec un maillage 100x4 éléments. Les 2 fichiers .info ne diffèrent que par la fonction f_propCinetique qui gère le basculement en mode visqueux au 1er incrément en fonction de la norme_de_convergence.
Le calcul avec amortissement mixte dure quand même 20min...
Si j'augmente la précision de convergence du calcul mixte je converge vers le cas RD_cinématique mais les temps de calcul sont plus long que le calcul de réf. !
Remarque : Une fois le basculement fait le calcul se poursuit avec amortissement visqueux. Comme l'a également observé Hugo, si l'on revient en mode mixte la convergence devient très très lente. Je ne comprends pas pourquoi !
Le lambda du cas non visqueux ne convient pas. Si tu fais tourner le calcul prends plutôt ce fichier avec lambda=2 pour le 2e incrément.
Temps de calcul 55min !
- % réalisé changé de 0 à 50
suite aux dernières simulations: en particulier celles effectués par Hugo, tout dépend de la précision utilisée. Si cette précision est suffisamment faible on obtient les mêmes résultats quelque soit les algo utilisés. Par contre pour une faible précision, l'amortissement visqueux conduit à une erreur de traînée : peut-être caractéristique d'un amortissement critique ??
Mais a priori rien d'illogique (??)
Formats disponibles : Atom
PDF